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NASA AMES THREE-DIMENSIONAL POTENTIAL FLOW ANALYSIS SYSTEM (POTFAN) 


EQUATION SOLVER CODE (SOLN) VERSION 1 


J. E. Davis 

Computer Sciences Corporation 

W. S. Bonnet t 
Acurex Corporation 

and 

R. T. Medan 
Ames Research Center 


SUMMARY 


This document describes a computer program known as SOLN which has been 
developed as an independent segment of the NASA-Ames three-dimensional poten- 
tial flow analysis system (POTFAN) , and which is used to solve small to 'large 
systems of linear algebraic equations by any of several methods including LU 
decomposition, Householder’s method, a partitioning scheme, and a block suc- 
cessive relaxation method. Due to the independent modular nature of the pro- 
gram, it may be used by itself and not necessarily in conjunction with other 
segments of the POTFAN system. 



1 INTRO DU CTIO N 


This document describes version 1 of an equation solver 
computer code (SOLS) which is a segment of the NASA^Ames 
three-dimensional potential flow analysis system (POTFAN) . 
This segment or the system selves the set of linear 
algebraic equations that are generated by the collocation 
method of satisfying the boundary condition cf a specified 
flow at various locations oh components in the flow field 
under consideration. 

The sets of constant coefficients and right hand sides 
associated with these simultaneous eguations are calculated 
tv other programs in the POTFAN system and transmitted to 
the S01N program as files through auxiliiary storage 
devices. The 30LN code reads in these files, determines the 
solutions, and then writes them cut as files to be read in 
by the next program in the POTFAN system. See Fig. 1-1. 

The SOLb! code provides a variety of solution techniques 
inclu'ding LU decomposition. Householder’s method, a 
partitioning scheme, and block successive over- relaxation. 
These various solution techniques are provided so that each 
.different type of system of equations encountered in 
potential flow aerodynamics may be handled in the most 
efficient manner availaole. The iU decomposition method is 
the fastest technique available for solving a small, closely 
coupxea, square system of equations that car. be placed in 
core all at once; the Householder procedure is the best 
method for solving an overdetermined .system of equations; 
the partitioning scheme is best for large closely coupled 
systems of equations that cannot be placed. in core ail at 
one time; and the block successive over-relaxation procedure 
is the potentially fastest method available for large, 
diagonally dominant systems of equations such as those 
generated in solving the problem cf an aircraft component xn 
a wind tunnel. 

The number of eguatiens that can be simultaneously 
solved by the code is mainly limited only by the amount of 
tape, disc, or drum storage available tc the user. 

The SOLN code i's constructed in a modular fashion so 
that any modifications or improvements to particular 
portions of the code do not affect the rest of the code. 
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2 _PBO B LEB T ASK DES CR IPTI ON 


This computer program was developed under Task. 3 of 
NASA-Ames Contract NA52-7571, and Task 26 of HASA-Ames 
Contract NAS2-6912. The purpose of these tasks was to 
develop a computer code which will, efficiently solve the 
systems of linear algebraic equations arising from 
attempting to satisfy the boundary conditions of a speciried 
flow on the various components in the flow field by the 
collocation method. These systems of equations may or may 
not exceed- the size of available computer storage, aiay be 
partitioned in an arbitrary manner with each partition 
representing the influence of one component on itself or 
another component, and may be overdetermined in some cases. 
In addition, a block iteration procedure was to be provided 
•where each block represents the aerodynamic influence of a 
single aircraft component on itself or on another aircraft 
component. This iteration procedure must allow an add-on 
solution capability without unnecessarily repeating any 
calculations (for example, obtaining the solution for a wing 
elone and then adding in wind tunnel walls tc determine the 
wall correction) and must allow an assemblage of components 
(for example, wing t body * tail) to be solved iteratively. 
Therefore, a variety of solution procedures were to be 
developed to best handle each type of situation likely to be 
encountered. — 

Also the program was required to be modular in nature 
so that it could be used independently and so' that any 
modifications or improvements to the code would not affect 
the other segments of the POTFAN system. 

Furthermore, the code was tc be constructed so that it 
is versatile, yet easy to use and easy to modify. 

Finally the program was required tc he able to handle 
large problems consisting of many separate components and/or 
large influence matrices. 
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i METH OD OF SO LUTIO N 


This section describes how the problems posed in tie 
previous section were solved. 


1.1 EQUATION SOLVING OPTIONS 


The SOLN program was designed to compute the solutions 
:f linear algebraic systems. It can compute these solutions 
:y one of several different procedures depending on the 
articular circumstances of the problems. In general, the 
ystem of equations to be solved may be overspecified (more 
guations than unknowns) or properly specified. These 
ypes of systems produce a rectangular or square matrix, 
espectively. The SOLN code uses a different solution 
rocedure for each. In addition, differert procedures are 
sed depending on the size of the system of eguations to be 
olved. one type of procedure is used for systems of 
guations that are small enough to contain all necessary 
ata within allowable core storage at one time, while otner 
rocedures are used when the system of equations is toe 
arge to reside in core all at once. The latter procedures 
ake extensive use of temporary storage devices such as 
ape, disc, or drum files. A special block iteration 
rocedure is also available which best handles problems 
here certain blocks of a matrix are only weakly coupled 
ith - other blocks. This iteration procedure also operates 
ifferently depending on whether the system to be solved is 
ectangular or sguare in nature, and whether each block fits 
n core at one time or not. 

A detailed discussion of each of the solution 
rocedures available in the SOLN code is contained in tne 
allowing subsections. In the following discussion, it will 
e assumed that capital letters stand for matrices of at 
sast two columns and rows, and that lower case letters 
tand fer column vectors. In that case, the problem to be 
alved can be stated as: given A and B, find X such that 


AX = B or Ax = b. 


(3'- 1-1) 
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The only assumptions that are made about A are that either A 
cr the transpose cf A times A is nonsingular. It is also 
assumed that A, X and B have real elements. No structure 
assumptions have teen made about A: it is' not handed in any 
systematic manner, it is not symmetric or positive definite 
except by chance, and it is not even necessarily square, 
although squareness does cut down on execution time. 

There are four main procedures used in this code, the 
XU Decomposition method (or the Crout-Doolittle method) , the 
Householder method, a matrix partitioning scheme for large 
matrices that was developed especially for this program, and 
a block Gauss-Siedel iteration procedure. The references 
provided in the reference secticn are very helpful in the 
matrix theory involved. Forsythe (1967) ' has a Doolittle 
method written in Algol 60 and Fortran, and provides a 
glimpse of error considerations. Nering (1963) gives a 
complete development of matrix theory in a theoretical 
manner, while Ralston { 1963) provided the basis for the 
presently included decomposition subroutine. Westlake ( 1968) 
is perhaps the most valuable for a user interested in 
selecting and implementing his own procedure. 


3.1.1 The Decomposition flethod 


This is one of the best all around methods ror a 
general, nonbanded, nonsymmetric real matrix. It is a 
method of solving a set of equations without calculating 
either the inverse of the matrix or any sequence of large 
matrices. The procedure is approximately as fast and just 
as accurate as any eguaticn solving technique - so far found. 
This procedure also gives a clear-cut singularity indicator 
before the entire procedure has teen executed. The double 
precision version of this routine is very readily arrived 
at, if desired, as well. 

The method Is based on decomposing a given matrix A 
into the matrix product of a unit lower triangular matrix L 
and an upper triangular matrix U, so that 

A = XU. (3.1. 1-1) 

Cnee this. has been accomplished, the matrix problem becomes 
a set of two easily solved eguatiens 


Ly = b (3.1. 1-2) 

and 
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Ux = y 


(3.1. 1-3) 


The decomposition is almost as easily done as said. To 
produce the first row of U, the first row ci A is taken, 
since the components of the first row of 1 (which is, l f! = 
1, l )a = 0, 1|3 = 0, . ..l (f1 = 0) times the columns of U 

produce the first row of U identically. And, by the 
eguality statement, the' first row of U equals the first row 
of A. Now, with the first row cf U available, the first 
column of L can be calculated. This proceeds to the nth row 
of [J, alternately calculating rows of 0 and columns of 1 
with data provided l>y the previous rows and columns. The 
equations at the rth step are 


(here 


ti 


TI 


r-| 

Z 1 rfe u ki 

R=1 


(3.1. 1-4) 


mr 


( a r)T 


u kr ) / u rr 
k=i 


(3.1. 1-5) 


This decomposition procedure breaks down at any point 
it which the diagonal element u rr vanishes. To prevent 
:his, each time the procedure starts a new row of U , it 
learches the column under consideration for the largest 
;lement. Having, identified the largest element on the 
liagonal, the row about to he reduced and the row containing 
.he largest element are interchanged. Thus, if a column is 
ver searched, and nothing but zeros are found, the 
rocedure stops, since a singular matrix has been found. 

The decomposition of the matrix takes roughly n 3 /3 
ultiplications, contrasted tc the usual n 3 operations of 
auss-Jordan reduction, and is therefore much faster. Since 
he inverse has not been explicitly determined, that extra 
torage has not been needed since the intermediate results 
f the decomposition are stored where the matrix wasl 

The auxiiliary equations are solved very quickly 
ecause they are triangular. First the vector y is 
etermined, then the solution vector x. The first component 
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The second component 


of y is b* , the first component of b. 
of y is determined from 


1 21 I| + ?2 ^ z (3* 1_» 1~ 6) 


and so on. Once y has been determined, the solution vector 
is J *unzipped M from 0 in the same manner, only starting from 
the bottom. 

In summary, the equations to be solved are 


A = LU 


(3.1. 1-7} 


Ly = b 


(3.1. 1-8) 


Ox = y 


(3.1. 1-9) 


3.1.2 Householder’s Method 


This method involves triangularization of the matrix 
and the use of elementary Hermitean unitary transformations. 
At present, the operations done to the matrix to 
triangularize it must also be done to the constant vectors 
simultaneously. In other words, unlike the decomposition 
method, there are no provisions for storing the' operations 
done to the matrix so that they may be done to constant 
vectors independently ' at a later time. Every constant 
vector must be present at the time of triangularization, or 
the original matrix must be reread and xetriacgularized . 

The method involves premultiplying a matrix by a 
sequence of other matrices, all of which are unitary 
matrices, so that the condition number of the matrix is 
unchanged at each step. The condition number is a function 
or the size of the matrix inverse, ‘and determines how badly 
roundoff effects confuse the answer. It, can' be shown that 
the condition number does not change as the matrix is 
manipulated and therefore gets no worse. 

Basically the Householder method seeks to triaiigularize 
the matrix A, by multiplying it by a series of the unitary 
matrices, called P| . for the first step we would liKe 
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(3. 1.2-1) 





here the first cclunin of A (call it a) has a nonzero 
leaent in the first row, but zeros in all subsequent rows 
f this cclunn. 


Pa = ke, where e, = (1 , 0, 0, G. . . 0) 


(3. 1.2-2) 


he size of ’’a” can be detersined in the usual method; 


I j a | 1 


^ 2 

5Z < a k, > 



s 


(3. 1.2-3) 


=> I tall 2 = s 2 


(3. 1.2-4) 


: can be shown that the size of a unitary transformation of 
vector is the same as the vector itself; 


l!Pa|| * 


n_ 

> 

kli 


but 


H ) 


= s (3. 1.2-5) 


I I Pal I 2 - - like, || 2 - k 2 


(3. 1.2-fc) 


and therefore: 


k 


= +s 


(J. 1.2-7) 


It can be shown that for any vector w such tnat J|w|l = 
and w is real', then a matrix constructed bj 


* 


~-5 



a 


= I - 2ww T 


(3. 1.2-8) 


is unitary. The choice of w with which to construct a 
matrix of this form so that it has the properties of f is 
the key to the procedure. 


Pa = (I - 2vw r )a * a - 2w (w T a) ' (3. 1.2-9) 


but the matrix product w^a = K, is a constant. All elements 
other ,than the first of this resultant vector should vanrsh, 
since the first column is being triangularized- That 
implies the following set of equations: 

a; - 2Kwj = 0 i?1 (3.1.2-10) 


a, - 2Kv, - k = ±s (3.1.2-11) 


or 


wj = a /2K i ^ 1 and w, - (a { + s)/2K (3.1.2-12) 


but, J| w jl =1, by the requirements for constructing the 
matrix, and therefore 




1 



(a 


n 



fe =2 



(3, 1.2-13) 




2a, s • *■ s 


2 


+ 




(3. 1.2-14) 
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(3. 1.2-15) 


4K* 


| i a) | 2 + 2aj s + s 2 


4K? 


s 2 + 2a, s + s 2 


13. 1.2- 16) 


4K* 2 


s * a ( s 


(3. 1.2-17) 


Therefore; 


2 . 5 

2K = s + a ( s 


(3. 1.2- 18) 


f the vector u is defined as; 


u = ( (a ( + s) ,a z ,a g .. .a n ) 


{3. 1.2- 19) 


ad w as w = u/2K, then the following relations can be shown 
3 be true; 


a) w T w — 1 {3. 1.2-20) 


b) P = I - 2ww T is unitary (3. 1.2-21) 

and 


c) Pa = +se, . (3. 1.2-22) 

At this point, the first column of A has been 
:iangularized. The other columns of the matrix are then 
aerated on by the unitary matrix (or its vector form) , 
.though none of them are triangularized if the matrix is 
.ngular. The next step in the procedure is to develop the 
•ctor w for producing the same triangularization of the 
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second column, from its diagonal element cn down. The third 
column will be treated from the third row down, and so on. 
In this way, a series of matrices are generated which 
premultiply the’ original matrix. The (n - 1) +h matrix 
produces an upper triangular matrix. If the same 
transformations have been done to the constant vectors, then 
all that need be done is a back-substitution to produce the 
solution. ’ 

Another important feature of this procedure is the fact 
that no division by diagonal elements, and therefore no 
pivoting of rows, need be done. .also, -this procedure can be 
used to produce "least-squared” solutions if the matrix has 
more rows than columns, without explicitly multiplying the 
matrices. In ether words, the system 


(A T A) x = a T b (J. 1. 2 - 23 ) 


can he solved without actually multiplying the matrices A . 
and A together. 


•3.1.3 The Partitioning Schemes 


There are several possible methods for partitioning a 
very large matrix into submatrices in order to create the 
matrix inverse. Most of these methods could be termed 
"block” methods, in that they apply well defined methods 
for dealing with individual elements cf a matrix to 
submatrices.’ for instance, a very popular method- for 
solving a predictably sparse, well behaved matrix system is 
the method of successive over- relaxation. The block analog 
cf this is the method of block successive cver-relaxation 
described in section 3.1.h. 

In this section, a partitioning scheme designed to take 
maximum advantage of available cere storage is described. 
In this manual it is known as the large decomposition 
solution procedure. The partitions developed in this 
algorithm are not related to the natural partitions 
associated with the separate ' components of the 
configuration. 

Basically, the method starts in the upper left corner 
of the large, dense matrix and inverts directly as large a 
block as is possible. From that point on, pieces of the 
original matrix are processed as the matrix inverse "grows” 
toward the lower right corner. Ihe number of words in each 
partition is as large as possible. The sequence of the- 
inversion is shown in Figure 3. 1.3-1, and the sequence of 
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quations at each step are gives by 


B is replaced by AB 
D = B - CB 


(3. 1.3-1) 


D = D ' 

E = — BD 
C = CA 
A = A - EC 
C = -DC 


(3. 1.3-2) 
■ (3. 1.3-3) 
(3. 1.3-4) 
(3. 1.3-5) 
(3. 1.3-6) 


Let the matrix be partitioned as is shown in Figure 
.1.3-2. The inverse of that matrix is assumed to be ot the 
ollowing form 


LM 

UP 


{3. 1 . 3-7) 


Since the product of the matrix and its inverse is, the 
dentity matrix, the following set cf eguations can be 
evelo pea 


P - (D - CA f B) 


(3. 1 . 3-8) 



(3. 1.3-9) 


N = -PC A 


(3. 1.3-10) 


L 




(3. 1.3-11) 


le program eguations are identical except that they are 
^signed to take as little extra storage as possible. 
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3.1.4- Block Successive Over-Relaxation 


The block successive over-relaxation procedure provides 
an alternate method for. the solution cf large systems of 
equations- It is particularly advantageous for matrices 
that are block diagonally dominant. This situation occurs, 
for example, in the case of a wing in a wind tunnel. The 
method consists of an iteration sequence which operates on 
an initial guess to the solution cf the equation 


Ax = b 


(3. 1.4-1} 


where the matrix A consists of several submatrices 'A 
submatrix represents the influence cf an aircraft componen 
cn itself or on another aircraft component. The figure 
ielow illustrates this matrix structure. 


*tz 


A zi 



(3. 1.4-2) 


A, 


m 


^nn 



Ihe subscript n refers to the number of subvectors Xj , each 
cf which represents the singularity strengths of a single 
component. ■ 

In order to begin the iteration procedure, an initial 
guess for each- of the subvectcrs xj must be either given as 
input data or computed or both. If ' the initial guess is to 
be entirely computed by the program, the jirocedure begins by 
approximating the solution for the first subvector x“ as the 
solution to the system of equations 
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ff s* 



( 3 . 1 . 4 - 3 ) 



where the superscript refers to the iteration number ana the 
superscript "o n refers to the initial guess. The initial 
guesses to the succeeding subvectors x» are computed by 
sequentially solving each of the following systems of 
equations; 


A 22 X 2 




( 3 . 1 . 4 . 4 ) 




Vl Vl 


At any step in the above procedure, the initial guess 
for the sub-vector under consideration may be entered as 
input data rather than being computed. Once an initial 
guess has been provided for each solution sub-vector, the 
succeeding iteration proceeds by sequentially solving each 
cf the following systems of equations: 


\ x i * B i 


. — ft v-M 
X 2 



( 3 . 1 . 4 - 5 ) 


A x. = b. 
22 2 l 



\ 


jtr 


23 





A x = n 

mm m u m 


mi 






A mm+i' 


x k '' - 
X m-H * * * 
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After solving each system of equations, the solution 
sobvectors may be further adjusted by the method of 
successive o ver-relaxaticn or under- relaxation. That is, 
let xf be an arbitrary solution sub vector of the kth 
iteration step. Then define, the improved estimate of that 
solution subvector as 


xj* * xH + w(xj k - xH ) (3. 1.4-6) 


That is to say, the improved estimate at step k is 
extrapolated from the Gauss-Siedel estimate and the previous 
improved estimate. If w=1, the method reduces to that of 
Gauss-Siedel. The quantity w is called a relaxation 
parameter, the choice of which determines the rapidity of 
convergence. for problems where the successive iterations 
indicate an oscillatory convergence behavior in- a 
predominant number of the elements of the solution 
subvectors, the rate of convergence will probably be 

increased by setting w to some value between 0 and 1.0. If 

the convergence behavior is predominantly cverdaraped, the 
rate of convergence may be increased by setting w to some 
value between 1.0 and 2.0. For a given problem there exists 
an optimum w such that the number of iterations to 

convergence is minimized. Several procedures have been 
suggested for determining this, optimum w. Ames(1969) and 

Eratkovich (1975) provide further information on determining 
an optimum relaxation parameter. In the SOLN code, the 
value of w is input by the user as part of the input data. 

Convergence of the iteration procedure is defined such 

that 


I xf - xf-' I 


< S 


|x!"i 


(3. 1.4-7) 


for each element in each subvector x ( > . If the specified 
maximum number of iterations have been performed and 
convergence has not been achieved, the last estimates of the 
solution subvectors may be saved and used as initial guesses 
for a subsequent series of iterations. 
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In solving each system of equations in' the series of 
systems described by equations 3.1. 4-4 and 3. 1.4-5, any one 
of the direct solution techniques available in the SOLN code 
may be used. To eliminate redundant calculations in each 
iteration step, the inverse or decomposition of the diagonal 
submatrices A*j is computed or read in and stored before the 
initial guess or iteration procedure begins. Then the 
matrix multiplications or forward elimination and backward 
substitution are the only computations made in solving each 
system or equations involved in the initial guess and 
subsequent iterations. The decomposition is used with 
square systems of equations for each diagonal sub-matrix 
that fits in core at one time. The inverse is used in all 
ether circumstances. This feature significantly reduces the 
amount of computet time required. 

This iteration procedure may also be used to solve a 
rectangular system of equations. In this case, _ the 
Householder procedure must be used to compute the inverses . 
of the diagonal submatrices. 


3.2 SYSTEM MODULARITY 


The required degree of modularity, which implies that 
the program should operate independently from and without 
interferring with other portions of the POTPAN system, has 
been guaranteed by the use of auxilliary storage devices 
(tapes, disks, or drums) as the only method by which the 
various segments of POTPAN can communicate with one another. 
This has the disadvantage that communicating through 
auxilliary devices is relatively slow and on some computer 
systems will require many job control cards to manage the 
devices. However the d isadvantages of this approach are 
more than offset by the advantages. The principal advantage 
is that it strictly guarantees segment independence. Also, 
this approach is necessary to maximize the size of problem 
that can be solved with a fixed amount of core memory. 


3.3 VERSATILITY AND EASE Op USE AND MODIFICATION 


These factors are all-important to the usefulness of a 
program. Although there is no unigue formula for 
guaranteeing them, two important elements of a sufficient 
formula relevant to the SOLN code have been identified. 
These are command format programming and the liberal and 
thoughtful use of comment cards in the source deck. 

Command format programming is a phrase coined to 
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describe a programming method which uses words or acronyms 
(called commands) to control the actions taken by the 
program. Although this technique is not new, the particular 
style employed in POTFAN programs originated with the 
computer program reported in Medan (1 973) . This technique 
contributes significantly to the above important factors. 

The manner in which command format programs work is the 
hollowing. First of all, there is a . program known as the 
control program, which may be either the main program or the 
principal subroutine called by the main program. For the 
SOLN code the control program is the subroutine SOLNIO. The 
first two actions performed by the control program are the 
initialization of default values and the establishment of 
whether or not the program is being run in a batch or 
conversational mode. In the conversational mode the program 
prompts the user for commands and data, and in the event of 
a recoverable error, pauses to allow the user to perform a 
fixup. In the batch mode the program echoes each command 
read in and generally always stops upon detecting an error. 
Other than these differences, the two modes of operation are 
identical. Following these actions the program enters the 
command phase. In this phase the program reads in various 
four character commands, takes the action associated with 
the command, and then (provided the command is not STOP) 
reads in the next command and so on. The specific commands 
available are given in section 5. 2 . 

Regarding the use of comment cards, there are no 
specific' rules that can guarantee that the inner workings of 
a program are adequately documented. However, some 
guidelines have been developed. 

The first of these concerns the quantity of comment 
cards. Some programs that were thought to be-adeguately 
commented were examined and it was found that approximately 
25 per cent of the source lines were significant comment 
lines.' Most programs are not considered to be well 
documented unless they contain approximately this ratio of 
comment -lines.' Therefore, the SGLN code is documented to 
this extent. 

The second guideline concerns subprogram documentation. 
The comment cards should state as a minimum the basic 
purpose of the subprogram and all inputs and outputs except 
printed output. Each subprogram in the SOLN code was made 
to comply with this guideline. 

In summary, the use of command format programming and 
careful commenting have made the SOLN code versatile, easy 
to use, and easy to modify. 
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3.4 MAXIMIZING PROGRAM SIZE 


Two major techniques were used for treating problems 
involving a large number of components and/or components 
with large influence matrices. These are the -use of 
auxilliary storage units and dynamic memory ■ allocation. 


3.4.1 Auxilliary Storage 


The degree to which auxilliary storage units are used 
in the SOLN code depends on the size and type of problem for 
which a solution is desired. for - systems of equations that 
can be placed in core all at once, a minimum of auxilliary 
storage is required. In fact, the influence matrix and 
right hand sides are placed in auxilliary storage only when 
they are first read in. • When the code automatically 
determines that the system can be solved entirely in core, 
the influence matrix and right hand side are read back into 
core and solved. 

For larger systems of eguations that cannot fit within 
core simultaneously, the influence matrix must be 
partitioned and solved -in a piecemeal fashion which requires 
extensive usage of auxilliary storage units. For problems 
involving the iteration procedure, auxilliary storage units 
are used to store each block of the matrix as well as the 
right hand sides, the current estimate of the solution 
suivectors, and the previous estimate of the solution 
subvectors. 

This usage of auxilliary storage units allows for the 
solution of systems of equations whose size is only limited 
mainly by the amount of auxilliary storage available rather 
than the amount of core storage available. 


3.4.2 Dynamic Memory Allocation 


Dynamic memory allocation as used herein refers to 
packing several small arrays into one large single array. 
In this way all of the available core space can be utilized. 

A necessary procedure in the use of dynamic memory 
allocation is that dynamically allocated arrays must be 
passed as variably-dimensioned arrays through subroutine 
argument lists. This in turn makes the program easy to 
modify since the storage array into which the smaller arrays 
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are packed needs to be dimensioned only once, e.g., in a 
small main program. 

The storage array in SOLN is denoted as (A) . This 
array is packed differently depending upon the type of 
calculations being performed. When bulk data is read in, 
storage is only allocated for that particular set of data. 
Then the data is placed in temporary storage. The 
allocation of storage is more complicated at the point where 
a system of equations is to be solved, storage in this case 
is allocated depending upon the method of solution, whether 
the iteration procedure or the direct procedure is used, and 
whether the system of equations can bo solved in core all at 
once. In general, however, the arrays that are the length 
of the number of fragments and the number of rows in the 
matrix are allocated first. Then the remainder of the 
available storage is allocated for the influence matrix'. 
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4 PROGRAM D ESCRIPTION’ 


To a great extent the description of the inner workings 
of the program has been relegated to comment cards in the 
FORTRAN source decks. This includes descriptions of the 
functions of the subroutines and- their input and output. 
The remainder of this section presents relevant descriptive 
data which could not effectively be placed on comment cards. 


4.1 CALLING STRUCTURE 


Figure 4.1-1 shows the subroutine calling structure. 
Table 4.1-1 shows the calling structure in a different 
format. 


4.2 FLOW CHART 


Figure 4.2-1 is ' a flow chart of the control program 
SOLNIO . 


4.3 COMMON BLOCKS 


Table 4.3—1 shows the common blocks used, their sizes, 
and the subprograms which they appear in. 


4.4 LOGICAL UNITS 


Table 4.4-1 summarizes the logical units (tape, disks, 
cr drums) which the program uses. Note that not all units 
would be used for each program run. ,,For the worst case the 
number of units used would be ten (10) . The specific data 
input from or output onto each of these units except for the 
line printer unit is discussed in detail in sections 6 and 
7.- 


4-1 


REFROPUCSBILrrY OF THE 
0KH8NAL PAGE IS POOR 



4.5 SEHGBY REQUIREMENTS 


Without the working storage array, the SOLN code 
requires approximately 24000 decimal words of core storage. 
This requirement includes all system subroutines and 
internal symbol' dictionaries and was determined on the 
INFONET Univac 1108 operating system without using overlays. 
The size of the working storage array must be added to this 
number to determine the total amount of storage required by 
the program. The memory requirements can be' reduced by 
usiny overlays. 


4.6 RESTRICTIONS AND IIM1TATI0NS 


The mcst important restrictions regarding the SGLN code 
are those associated with the limitation in size of the 
system of equations that can be solved. There are three 
different types of limitations depending cn the manner in 
which the system is solved. First, there is a limitation to 
the number of equations that can fce solved in core at one 
time. For a given computer system, this limitation can be 
determined by the following formula 

NR = J 4 + <HXS-3*NIR) ’ -2 (sguare matrix) 


NR - MS ~ 3*NIB 
4 + NC 


(non-square matrix) 


where NR - Number of equations that can be solved in core 
at one time. 

MIS = Size of the working storage array (A). 

NC = Number of unknowns to be determined. 

NIR = Number of records of input containing the 
influence matrix. This is also referred 
to as the number of influence matrix fragments. 

Second, there is a limitation to the number of 
equations that can be solved even using auxilliary storage 
devices. This limitation depends upon whether the iteration 
procedure is used or whether a direct procedure is used. 
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for a direct procedure the limitation is computed as 
follows: 


NB = HXS-3*NIB+1 {sguare matrix) 

7 


BB = HXS-3*NIB 
4 


(non-sguare matrix) 


Por the iteration procedure the limitation is computed as 
fellows: 


NBB = «XS-3 *«IB +1 1 

? ) j 

] whichever is smaller 
j (sguare matrix) 

NBB = M XS -3*N IB - BAX V NF j 

4 J 


NRB = MXS— 3*NXB— BAXVN P 
4 


(non-square matrix) 


where NBE = The order of the largest diagonal nloclc 

of the matrix. 

MAXVNF = Number of influence coefficierts in the 
largest influence matrix input record. 

The third limitation concerns the amount of auxiliiary 
storage available. The temporary storage units INA, INAT, 
1NCBD, KST, INB, INBT, and INVB must each be able to contain 
the entire influence matrix. Hence, the largest system of 
equations that can be solved may be limited by the size of 
the auxiliiary storage devices. 
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5 OPER ATING IN STRUCTION S 


The purpose of this section is to provide the user with 
information necessary to execute the program. 


£.1 GENERAL DATA INPOT CONSIDERATIONS 


The input data for the equation solver program consists 
cf punched cards and tape, disc, cr drum files when the code 
is utilized in the hatch mode, and on-line data with tape, 
disc, or drum files when it is executed in the- 
conversational mode. All punched card and conversational 
terminal input is prescrifced in NAMELIST cr regular formats, 
whereas all bulk data input through tape', disc, or drum 
files is unformatted. The latter data is described in 
sections 6 and 7. 

The program is designed to use commands as the basic 
fcrm of input to control the program flew. These commands 
consist of four letter words placed in the first four 
characters of an input record (first four columns of an 
input card) and are recognized as keys that cause the 
program to perform particular operations. After the 
operations are performed, the program flew returns to the 
beginning of the program and reads the next command. This 
continues until a STOP command is encountered, whereby the 
program terminates. Any command input record (card) whose 
first four characters cr columns are left blank is 
considered a "comment” command. Any command that is not 
recognized by the program is printed and program flow is 
either returned to the next command without any operations 
being performed or terminated depending on the value of the 
variable CONTIN (see section 5.2.5). The particular 
commands recognized by the equation solver code are 
discussed in more detail in section 5.2. 

All of the commands available in the equation solver 
code are related in one way or another to the solution of a 
system of equations, although several commands in sequence 
are required to complete all the steps necessary to obtain 
this solution. Since there are several different solution 
techniques available in the code, the best method of 
solution and the associated commands for a given problem 
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will depend upon whether the system of equations (matrix) 
under consideration is over- deter mined (rectangular) or not, 
and whether the system of equations under consideration can 
te subdivided or blocked effectively so that the iteration 
procedure can be implemented. 

For solving larger systems cf eguations, the 
potentially fastest, solution technique is 'the block 
successive over-relaxation method- This is true whether the 
system of equations is over deter mined or not. However, the 
iteration procedure sill only converge tc a solution if the 
matrix produced by the system of equations has certain 
properties. In particular, the matrix must possess 
submatrix blocks such that the diagonal blocks are more 
dominant than the off diagonal blocks. The more dominant 
the diagonal blocks, the faster the iteration procedure will 
converge to a solution. 

Matrices that dc not possess this property are best 
solved by the direct technigues. For sguare matrices,, the 
IU decomposition method is the fastest cf the direct 
methods, and for rectangular matrices the Householder 
procedure should be used. In the event that the data 
defining the system of eguations under consideration 
reguires more computer cere storage than is available, the 
equation solver code automatically places the data in 
temporary storage and determines the solution to the system 
cf equations by the procedures for large systems as 
discussed in Section 3- 3. 


5-2 INPUT DESCHIPTION 


The first line of input tc the program consists of a 
single logical variable. Its value depends on whether the 
command inputs are to be made utilizing the conversational 
cr batch modes, -TRUE, fer conversational, .FALSE, tor batch 
mode. The format for this input record (card) is Li. 

The remaining input for the equation solver code 
consists mainly of operation commands, input file names, and 
option flags- A command statement must begin in the first 
column or with the first character of an input card or 
record. • Only the first four characters of any command are 
recognized by the program. If the first four characters or 
columns are left blank, the command is assumed to be for 
comment purposes, and the entire card or line is printed. 
In the conversational mode, an unrecognizable command 
returns control to the next command after the unrecognizable 
characters are printed.; otherwise, the program terminates 
unless the variable CONTJN has been set to .THUE. (see DATA 
command- ' A detailed description of the commands recognized 
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by the code, as well as the input associated with them, is 
given below. All logical unit numbers may range from 1 to 
59, except 5 and 6. File identification numbers may range 
from 1 to 9999. See Appendix A for an explanation of file 
identification numbers. 


5-2.1 (9 -Blanks ) Command 
Ex plan at ion 

Any command with blanks in the first four positions is 
considered a comment card. All information cn this card or 
input record is printed, then ccntrcl proceeds to the next 
command. 


5.2.2 B CB-S Command 
Ex planati on 

This command causes the boundary condition constant 
vector files to be read in and placed in temporary storage. 
This boundary condition constant vector file is ordinarily 
created by the boundary condition program. This command 
should not be given before the CNTL command (when required) 
or the VNRE command, and must precede the SGLV, RIVNS, or 
KBIT commands. 

In put 

( NTBCR , NBCR (I) , 1= 1 , NCMP) — The identification number and 
logical unit number of each boundary condition constant 
vector file. 

Fo rma t 

215 

C omment s 

Enter one identification number and one logical unit 
number per line or card of input. The input sequence must 
be consistent with the input sequence of the input influence 
sufamatrices (see VNRE command). 


5.2.3 C NIL Command 
Explanat ion 

This command defines the number of aircraft components 
and the overall type of solution procedure to be used 'in 
solving the system of eguations under ‘-consideration, it is 
the first command given if the internal logical units do not 
require initializing with the DATA command. This command 
must precede all other commands except blank and DATA 
commands except when there is only a single component whose 
influence matrix can fit in core. In the latter case, the 
CNTL command is not required. 
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Input 

NC HP" Number- of aircraft components associated with the 
system of. eguations under consideration. This determines 
the ‘number of boundary condition and influence matrix files 
that must be read in. 

ITEBAT — Flag which designates whether an iteration 
procedure or a direct solution procedure is to be used to 
solve this system of- equations. 

HE BAT = 

0 — Direct solution 

1 — Iteration solution 

F ormat 

215 


5.2.4 CO MB Command 
Ex planation 

This command is used only when the system of eguations 
under consideration is toe large for available core storage 
and the direct soluticn procedures are being used to solve 
the system. It combines all submatrices read in by the VNR'E 
command into one matrix which is stored in column length 
records on a temporary storage device. The command need not 
ie used unless both circumstances mentioned above exist. No 
■further data input is reguired with this command. If .used, 
it must immediately follow the VN BE command. 

5.2.5 DATA Command 
Ex planati on 

The DATA command is used to initialize any or all 
internal logical units used by the program. This includes 
any temporary storage devices as well as the input and 
output logical units. The DATA command is' also used to set 
the variable CONTIN. This command need not be used if the 

« default values of the variables are satisfactory. This 
command , if needed, roust precede all of the commands except 
the blank command. 

I nput 

NTCP — Logical unit number of the output device used for 
conversational mode operations. The default value is 6. 

NTP — logical unit number of the output device ■ used for 
batch mode operations. The default Value is 6. 

NTCB — Logical unit number of the input device used for 
conversational mode operations. The default value is. 5. 
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NTSE — logical unit number of the input device used for 
batch mode operations. The default value is 5. The value 
of NT B is set egual to NTRR immediately following the read 
statement. 

INA — Logical unit number for a temporary storage- 
device. The default value is 11. However, during a COME 
command, it may be changed to 12. See section 7.4 for an 
explanation of the temporary storage devices. 

INAT — logical unit number for a temporary storage 
device. The default value is 12. However, during a COME 
command, it may be changed to 11. 


INCBD — Logical unit number 
device. The default value is 14. 

for 

a 

t era porary 

storage 

KST — Logical unit number 

device. The default value is 13. 

for 

a 

temporary 

storage 

INE — Logical unit number 
device. The default value is 15. 

for 

a 

tem porary 

storage 

INBT — Logical unit number 
device. The default value is 17. 

for 

a 

temporary 

storage 

INVE — Logical unit number 
device. The default value is 18. 

for 

a 

tem porary 

storage 

NTIR — Logical unit number 

for 

a 

temporary 

storage 


device. The default value is 16. 

CONTXN — Logical variable governing program action upon 
encountering an unrecognizable command. Unless CONTXN is 
.TEUE. , execution will terminate with an unrecognizable 
command. The default value for CONTIN is .TRUE. under 
conversational operation and .FALSE, under batch operation. 

F ormat 

Namelist MDATA. 


5.2.6 IN VE Command 
E xplana tion 

This command causes the inverse of a matrix to be 
computed and stored on a tape, disc, or drum file. The user 
must input the identification number and logical unit number 
of the tape, disc, or drum file where this inverse is to be 
stored. In addition, the user must provide a title for the 
file. A special option with this command allows the user to 
compute and store the decomposition of a matrix instead of 
the inverse (see Section 3.1). The only restrictions for 
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the decomposition option are that the matrix to be 
decomposed must fit in core all at once and that it must be 
a square matrix. The INVE command must follow the VNRE and 
(if needed) COMB commands. 


METHOD — Option flag which designates which solution 
technique is to be used to create the inverse or 
decomposition. 

METHOD * 

0 — LU decomposition method (square matrix only) 

1- -Householders* technique (rectangular or square 

matrix) 

N1INVW — Identification number of tape, disc, or drum 
file where the computed inverse is to be stored. The 
logical unit number associated with this file is INVR which 
may be input with the DATA command at the beginning of each 
test case. 

I DC — Option flag which designates whether the inverse 
cr the decomposition of the matrix is to be computed. 

IDC = 

0 — Inverse 

1- - Decomposition 

NTITL — Number of words in the title describing the 
inverse or decomposition file. NTITL must be greater than 
2 ero but not exceed 100. 

(TITL (I) , 1=1, NTITL) — A one dimensional alphanumeric 

array that may be used to. describe the contents of the 
inverse file. 

f ormat 

3X5/15/ (20A4) 


5.2.7 IVNfi Command 
E xplana t ion 

This command causes an inverse or decomposed matrix 
file to be read in and placed in temporary storage and is a 
sequel to the INVE command. The intent of this command is 
to allow the user to read in an inverse without computing 
it, then multiply the inverse by any number of constant 
vectors to obtain the solution vectors. In the case of a 
decomposed matrix, a forward elimination and backward 
substitution are performed on the decomposed matrix and 
constant vectors to obtain the solutions. The IVNR command, 
however, only reads the inverse or decomposed matrix. Other 
commands read the constant vectors (BCRE command) and 
perform the multiplication or forward elimination and 
backward substitution (MULT command). The IVNH command must 


r* 
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f ollow a CNTL command and a BCBE ccmntarid and precede a EUL1 
command. 

Input , , 

NTINVR — Identification number of the trie where the 

inverse or decomposed matrix is stored. 

HTIHV — Logical unit number of the file where the 
inverse cr decomposed matrix is stored. If entered as zero, 
the program uses 18 instead. 

Fo rmat 

215 


£.2.8 MULT Command 

Explanat io n 

This command is 
command. It causes an 
any number of constant 
elimination and backwar 
decomposed matrix and a 
information stored on 
command indicates wheth 
existed on the file and 
appropriate computation 
Ihe IVNB and BCRE comman 


used in conjunction with the IVNK 
inverse matrix tc be multiplied by 
vectors or it causes a forward 
d substitution to be performed on a 
ny number of constant vectors. The 
the file that was read by the IVNIi 
er an inverse or a decomposed matrix 
the code automatically performs the 
s to obtain the solution vectors, 
ds must precede the MOLT command. 


Input . . 

(OUTA (I) , OUT AMT (I) , 1= 1 , NCM P) — One dimensional _ integer 

arrays defining the identification numbers and logical unit 
numbers of the files where the computed solution -vectors are 
to be stored. 


Format 

215 


£.2.9 PART Command 

Ex plana t ion . . 

This command causes the program to partition the matrix 

that was read in by the VNRE command and combined, if 
necessary, by the COMB command. It is a utility command 
that need net ordinarily be used unless- it is desiranle tc 
partition a matrix into submatrices in a manner that is 
different from the way the matrix is already partitioned. 
The PART command will cause the matrix to be subdivided into 
equal sized square blocks, if possible. If this is not 
possible, then the blocks will be made as equal m size a 
is possible. This command must follow the VNRE and (i 
used) COMB commands and precede tbe SOLV and INVE commands. 
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Input 

NP.fi — The number of partitions that a rcw or column of 
the matrix - is -to he subdivided into, Ihe number of 
partitions in the matrix will be NPA**2. 

(NP (I) ,1=1 ,NPA**2) — fi one dimensional array defining 
the logical unit numbers of the files where each subaatrix 
partition is to be stored. A total of . NPfl**2 entries must 
be made. Place the values for one row of partitions on each 
line or card. 

F ormat 

15/(1615) 


5.2.10 SGI V Command 

This command causes the program to compute the solution 
to the system of eguations as defined by the v NRE and BCfiE 
commands. The user need not define a solution technique for 
solving the system of eguations; the code will automatically 
determine this for him. However, if he chooses to override 
the automatically chosen technique, he may dc so. If the 
user chooses an inappropriate technique for a given system, 
the code will automatically choose the best technique and 
override his decision. This may occur if the user attempts 
to solve a rectangular matrix using a technique for square 
matrices. The various solution techniques available in the 
code were discussed in detail in Section 3. If the 
iteration procedure option is in effect, this command 
defines the method for solving the system of equations 
associated with each diagonal submatrix cf the blocked 
matrix. The SOLV command must follow the V N RE , COMB (if 
used), and BCRE commands. If a user input solution 
technique is used, the calling sequence for this procedure 
should be placed in the subroutine INVRS at the designated 
location. 

If the iteration procedure is in effect and is 
n cn-con verges t , the last estimate of the solution subvectors 
will be printed and stored on files. The iteration 

procedure may be restarted by reading these final 

approximations into the program as initial guesses at the 
solution subvectors instead of have them computed. In fact, 
any predetermined solution subvectors may he read in as an 
initial guess for the iteration procedure instead of having 
them computed when the iteration option is in effect. In 
ether words this is the number components for which an 
initial guess of the solution will be supplied. 

Input 

METHOD — Option flag defining which solution technique 
will be used in solving the system of eguations. 

METHOD = 
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0 — LU decomposition method 

1— -Householders * raethcd 

2 — User input method 

IPS — The number of predetermined solutions bo be read 
when the iteration option is in effect. In ether words this 
is the number components for which an initial guess of the 
solution will be supplied. 

(OUTA(I) ,OUTANT (I) ,1=1, NOME) — Two one dimensional 

arrays defining the identification number and logical unit 
number where the solution subvectors cf the system of 
equations under consideration are to be stored. One pair of 
numbers is required for each aircraft component. 

IPBINT — A print parameter that prescribes whether the 
solution vector is printed out. 

IPBINT=0 — Solution vector is printed. 

=1 — Solution vector is not printed. 

(I,INPS(I) ,IHPSHT(I) ,J=1,IPS) — A Simple variable and 
two one dimensional arrays that indicate the aircraft- 
component sequence number, the identification number, and 
logical unit number of any predetermined solution subvector 
files to be used as an initial guess to the solution for the 
iteration procedure. This input is only required ix the 
iteration procedure option is in effect. If IPS is zero, 
this information must not be input. The solution subvectors 
cn these files will be used as the Initial guess for the 
iteration procedure instead of having the initial guesses 
computed. One set of these numbers is required for each 
aircraft component for which an initial guess is to be 
input.- 

MAXIT1 — The maximum number of iterations allowed for 
the iteration procedure tc converge. This is only required' 
if the iteration procedure option is in effect. If 
KAXTT1=0, only the initial guess will be computed, printed, 
and stored. 

OMEGA — Relaxation factor that alters the convergence 
characteristics of the iteration solution. 


0<OMEG A< 1 • (under-relaxation) 

KCMEG A<2 (over-relaxaticn) 

The default value is 1.0. This input is only required when 
the iteration solution option is in effect. 

ITEEP — Print parameter that defines the amount of 
intermediate results tc be printed during the iteration 
solution procedure. 
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ITEBP = 

0 — No intermediate results will be printed. 

Not-zero — Intermediate results will be printed every 
JTERP iterations. 

The intermediate results consist of the influence of 
each submatrix block on every ether sufcmatrix block, the 
constant vectors, and the most recent estimate of the 
solution. If XTERP 0, these quantities are also printed 
for the initial guess at the solution. This input is only 
required if the iteration procedure option is in effect. 

f orma t 

(215) — lor METHOD, IPS, OUT A (I) , OUTANT (I) , IPRIHT. 

(315) — For I, INPS(I), INPSNT(I). 

(I5/F5.2/I5) — For MAX IT 1 , OMEGA, ITEBP. 

C omment s 

The inverse or decomposition of the diagonal 
submatrices associated with the iteration procedure may also 
te read in instead of being computed. This is discussed in 
the VNRE command. 


5.2.11 STOP Command 
Ex plan ation 

This command terminates the program. No further 
operations can follow the STOP command. 


5.2.12 VNRE Command 


Ex planation 

This .command reads the submatrices making up the 
influence coefficients for the system of equations under 
consideration ^and places them in temporary storage. The 
matrix is usually subdivided into submatrices that are 
placed on separate tape, disc, or drum files. For potential 
flow problems, these submatrices usually represent the 
influence of a single aircraft component cn another single * 
aircraft component. - If the iteration procedure is in 
effect, each submatrix also represents a block of the block 
iteration procedure. A special option is available with 
this command when the iteration option is in effect. It 
allows the user to avoid calculating the inverse or 
decomposition of any one of the diagonal submatrix blocks 
that would ordinarily be computed during the iteration 
procedure. Instead, the user may read in the inverse or 
decomposition itself. For those diagonal submatrices that 
are read in the ordinary manner, the inverse or 
decomposition is computed during the iteration procedure 
calculations. The VNRE command must precede the BCRE 
command and the COMB command (if used) and must also follow 
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the CNT1 and DATA commands. 


Input 

NT V NR , NT VN — The identification number and logical unit 
number of a single submatrix cf the influence matrix file. 
Cue pair of numbers is entered per line. The submatrices 
must be entered in rowwise order. That is, referring to the 
figure in Section 3.1.4, the submatrices must be entered in 
the following order: A (| , ...,A ,A^ ) ,...,A 2n , 

...,A n| ,... - 

INVF(I) — Option flag that designates whether the 
associated file contains a sutmatrix or whether it contains 
the inverse or decomposition of that submatrix. This is 
cnly input with the diagonal submatrices and is entered on 
the same line as the logical unit number or identixication 
number fox that sutmatrix. 

IN VP (I) * 

0 — Submatrix. 

-1 — Inverse or decomposition of the submatrix. 

Forma t 

215 {tor off-diagonal submatrices) 

315 (for diagonal suimatrices) 


5.3 SISTEH CONTROL CARDS 


This section describes the control cards that are 
necessary to run the SCIN code on the various computer 
systems that have been cr are being used tc run the— SOLN 
code. 


5.3.1 infonet univac 1108 system 

Since this system allcws automatic file definition 
commands determined from the file identification numbers 
{see Appendix A) , the only control card required is the 
program name, SOIN/PGTF. All files created except scratch 
files will automatically show up in the user's catalog. The 
names of all created files will be identical to the input 
file identification numbers. Also each created file will 
have a version identifier that defines the type of file that 
it is. The version identifier for solution files is SO; for 
inverse files it is IN. All files created will reside in 
the LIB$ library. Thus, fcr example, a solution file 
created with an identification number cf 1023 will be 
assigned the name 1023. SC/LIB$. 


4 
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6 PBOG BA B INPUT BINA R Y FI LES 


The SOLN code uses two types of files created by other 
EOTFAN segments. These files contain the coefficient matrix 
(influence matrix files) and the right hand sides (boundary 
condition files) . These files- are accessed according to the 
procedures in Appendix A and Section 5.3. Each of these 
files conforms to the standardized POTFAN fcrmat, which is 
discussed in Appendix B. The notation used in the remainder 
cf this section will be clear to the reader if he reads 
Appendix B. All of the data available cn these files is not 
used by SOLN. The following subsections describe the files 
briefly and which quantities affect the SOLN code. The 
program may also input an inverse file, a decomposition 
file, or a solution file that was previously computed by the 
SCLN code. These tiles are discussed in Section 7. 


6.1 INFLUENCE HATBIX FILES 


Each influence matrix file contains all the information 
necessary to describe the boundary condition influence of 
cne aircraft component on some other aircraft component. 
The file itself consists of two or more records of 
information. The first record is referred to as the 
introductory record, while all subsequent records contain 
the influence coefficients. In many cases the influence 
matrix is very large and cannot be input or output all at 
cue time on one record of information. Therefore, several 
records are often required to Input or output the entire 
influence matrix file. Each record of influence 
coefficients is referred to as a matrix fragment. Each 
influence matrix file is referred to as a submatrix. , Many 
influence matrix files (submatrices) may be required to 
completely define the entire influence matrix associated 
with a given flow field problem consistiry of several 
aircraft components. 

The following information from the influence matrix 
files is used directly by the SOLN code. The variables 
marked with an asterisk are necessary for proper execution 
of the program- All variables without an asterisk are 
informative only and are not necessary fcr proper execution 
cf the SOLN program but are transferred tc the solution 



files to be used as input for subseguent programs in the 
FOTFAN system. 

1s t Rec ord 

* NCTIHE — Number of words in the data and time array, 

{CTiME (NCTIME) ) — Array containing the date and time of 
creation of the influence matrix file. 

*NI1TL — Number of words in the title of a 
self-influencing submatrix file. 

(Till (NTITL) ) — Title of 'a self- i n fluencing- sunraatrix 

file. 

*NK£CS — one plus the number of submatrix fragments. 

ID ( 1) --Identification number of a seif-inf luencmg 
submatrix geometry file. 

ID (3) — Identification number of unconstrained 

self-influence submatrix. 

ID {5) — Identification number of constrained 

self-inf lue nee submatrix. 

*NiOG 

LOG (2) =TOP — Flag indicating whether an upper surface 
boundary condition has been applied. 

L0G(3)=B0T — Flag indicating whether a lower surface 
boundary condition has been applied. 

LOG (11) =DBLT — Flag indicating whether doublet type 
singularities have been used to represent an aircraft 

component. 

LOG (12) =SOOBC£ — Flag indicating whether source type 
singularities have been used to represent an aircraft 

component. 

LOG (13) =CNSTRN — Flag indicating whether the influence 
matrix is constrained. 

*NINT 

*INT (8) =NSING — Number of unknowns associated with the 
self-influence matrix file. 

INI (12) =ICTYPE — Not used anymore. 
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INT(14)=NFC — Number of constraint function coefficients 
representing the singularity distribution of a 
self-influence submatrix file. 

INI (16) =NCFFX — Number of constraint function fragments 
in the N1 direction. 

XNT (17) =NCFFY — Number of constraint function fragments 
in the N2 direction. 

1NT(18)=ICV — Flag indicating the type of constraint - 
variables used to generate the constraint function. 

NF1T 

*FLT (1 ) =XMACH — Mach number. The SOIN cede checks all 

influence matrix files tc see that they correspond to the 
same Mach number. If there is an inconsistency, the program 
will step and print a diagnostic message. 

2n d Re cord an d Su bse que nt Records 

*J1=NC — Number of columns in the sufcmatrix. 

*J2=NR; — Number of rows of the submatrix cn this record. 


J3= 1 

l 

*NW=NC#NR — Number of influence coefficients on this 
record. 

* (h (I) ,1— 1 , NN) — Influence coefficients on this record. 
These are stored columnwise. Khen multiple records are 
required to store the submatrix, then a certain number of 
complete ro.ws of influence coefficients are contained on 
each record but are still stored columnwise. 


6.2 BOUNDARY CONDITION FILES 


Each boundary condition file contains all the 
information necessary tc describe the influence of the 
uniform and nonuniform freestream and the rotation rates on 
the boundary conditions of a single aircraft component. The 
file itself consists of three cr more records of 
information. The first record is referred to- as the 
introductory record. Each group of two records thereafter 
describes all the How field parameters associated with a 
single set of boundary conditions as well as their influence 
cn this aircraft component. There may be several sets of 
boundary conditions cn a single boundary condition, file. 
The second, fourth, sixth, etc., records on the file- contain 
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the right hand side suhvector associated kith the system of 
eguations describing the specified flow houndary condition. 
The third, fifth, seventh, ninth, etc., records contain the 
parameters defining the freestreas conditions; that is, the 
freestreara velocity vector, the angles of attack and 
sideslip, the rotation rates, and the center of gravity 
location for the flow field. 

The following information from the boundary condition 
files' is used by the SGLN code. The variables marked with 
an asterisk are necessary for proper execution of the 
program. All variables without an asterisk are informative 
only and are not necessary fcr proper execution of the SOLN 
program but are transferred to the solution files to he used 
as input for subseguent programs in the EOIFAN system- 
asterisk. 

1 st Recor d 

♦NCTIBE — Number of words in the date and time array. 

(CUBE (NCTIBE) } — Array containing the date and time of 
creation of the boundary condition file- 

♦NTIT1 — Number of words in the title. 

(TI1L (NTITL) ) — Title of the boundary condition file. 

♦NBHCS — One plus two times the number of sets of 
boundary conditions. 

I D ( 1) — Identification number of the geometry file 
associated with this aircraft component. 

ID (2) — Identification number of the boundary condition 

file. 

• L0G(5)=DBLT — Flag indicating whether doublet type 
singularities have been used to represent an aircraft 
component. 

LOG (6) =SOORCE — Flag indicating whether source type 
singularities have been used to represent an aircraft 
component. 

LOG (7) =TOF — Flag indicating whether an upper surface 
boundary condition has been applied. 

LOG (8) =BOT — Flag indicating whether a lower surface 
boundary condition has been applied. 

*IN1(2)=NB0W — Number cf rows. 

*INT (3) -=NSETS — Number of sets of boundary conditions on 
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this file 


2n d, 4t h, 6 th, . . . Rec ords 

*J1=NB0H — Number of elements in the constant subvector. 

J 2— 1 

J3=1 

NW=l!iBOH — See above. 

* (A (NBOH) ) — Elements of the constant subvector. 

3r d , 5t h , 7th, .. . Eec ords 


J 1= 1 1 


J 2= 1 


J3=1 

NW=11 

A ( N W) --Array containing the freestream velocity vector, 
the angles of attack and sideslip, the rotation rate vector, 
and the center of rotaticn location for this set of boundary 
conditions. 
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7 PROGR AM OUTPUT 


Output from the program consists cf line printer 
output, data left on scratch files, and various POTFAN 
files. The line printer output is meant to be self 
explanatory and will not be discussed further. The POTFAN 
files that are created are created according to the 
procedure in Appendix A and conform to the format in 
Appendix B, Control cards for managing the scratch and 
EOTFAN files are given in Section 5.3. The POTFAN and 

scratch files are discussed in mere detail in the following 
subsections. 


7.1 SOLUTION FILES 


Each solution file contains all the information 
relevant to that portion cf the solution vector associated 
with a single aircraft component. Each such portion ot the 
solution vector is referred to as a solution subvector. The 
file itself consists of three or more records. The first 
record is the introductory record. The second record 
contains the parameters defining the . freestream conditions 
for all sets of boundary conditions for which a solution was 
computed; that is, the freestream velocity vector, the 
angles of attack and sideslip, the rotation rates, and the 
center of rotation location for the flow field. The third 
and subsequent records contain the solution subvectors, one 
record for each boundary condition. 

The following information is written on each solution 

file. 

1st R ecor d 

NCTIME’ — Number of words in the date and time array- 

(CTIHE (NCTIHE) ) — Array containing the date and time 
that the solution file was created. 

NTITL — Number of words in the title arraj- 

(TITL (NTITL) ) — Title of the geometry file for this 
aircraft component. 

NHECS — Two plus the number of sets of solution 

subvectors. 

{IF CRH (NHECS) ) — 0, 1, 1, . . . 

NID=5 
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ID (1) --Identification number of the geometry file 
associated with this aircraft component. 

ID (2) — Identification number of the self-influence 
submatrix file associated with this aircraft component. 

ID (3) — Identification number of the inverse file used 
to generate this solution sufcvector, if such a file was 
used. Otherwise ID (3) is zero. 

ID (4) — Identification number of the boundary condition 
file associated with this aircraft component. 

ID (5) — Identification number of this solution file. 

NLOG=5 

LOG ( 1) =DB1T--F lag indicating whether doublet-type 

singularities have been used to represent an aircraft 
component. 

LOG (2) =SOURCE — flag indicating whether source type 
singularities have been used to represent an aircraft 
component. 

• LOG(3)=TO.P — flag indicating whether an upper surface 
boundary condition has been applied. 

LCG(4)=B0T — Flag indicating whether a lower surface 
boundary condition has been applied. 

LOG (5) =C NSTRN — Flag indicating whether the influence 
matrix is constrained. 

NIN1=8 

INI (1) =ICIYPE — Not used at present. 

INI (2) — number of elements in the solution subvector 
(i-e. , number of unknowns) . 

INI (3) — Number of sets of boundary conditions. 

I NT (4) =HETHOD — Method of solution. 

INT (5) =NSING — Number cf elements in the solution vector 
if it were unconstrained. 

INT (6) =NCFFX — Number of constraint function fragments 
in the N1 direction (if the matrix is constrained) . 

INT (7) =NCFFY--Number of constraint function fragments 
in the N2 direction (if the matrix is constrained) . 
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INT(8)=ICV — Flag indicating the type of constraint 
variables used to generate the constraint function. 

NF1I=1 

FIT (1) =X£4ACH — Mach number. 

2nd_fiecord 

J 1 =NSOL-Number of solutions. 

J 2= 1 1 


J3=1 

NW=11*NS0L 

(UlNF (NSOL) ) , (VINF (NSOL) ) , (WINF (NSOL) ) — Components of 

the unit freestream velocity vector for each set of boundary 
conditions, 

(ALPHA (NSOL) ) — Angle of attack for each set of boundary 
conditions. 

(BETA (NSOL) ) — Angle of sideslip for each set of 
boundary conditions. 

(P (NSOL) ) , (Q (NSOI) ) , (HR (NSCLJ ) — Components of the 
freestream rotation rate vector for each set of boundary 
conditions. 

(BCGX (NSOL) ) f (RCGY (NSOL) ) , (FCGZ (NSOL)) — Components of 

the position vector to the center of rotation of tho flow 
field fen each set of boundary conditions. 

3rd an d S ubs equent Re cor ds 

J 1= NR- -Number of elements in the solution subvector. 


J 2= 1 


J 3= 1 


N W = NS — See above. 

* ♦ 

(X (NS) ) — Elements of the solution subvector for a 
particular set of boundary conditions. 
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7.2 INVIESE MATE IX FILES 


Each inverse matrix file contains all the information 
pertaining to the inverse cf a matrix that was assembled and 
computed from one or more influence matrix files. That is, 
the system of equations from which this inverse matrix was 
generated may have involved cne cr more aircraft components. 
The file itself consists of .three or more records- The 
first record is the introductory record. The second record 
contains all of the information from the influence matrix 
files that must be transferred to the solution files when 
the inverse is multiplied 'by a constant vector to produce a 
solution vector. The third and subsequent records contain 
the inverse matrix coefficients. Since the storage required 
to output the entire inverse matrix on a single record may 
be more than is available, several records may be necessary 
to output the complete matrix. 

The following information is written 
matrix file. 

1st Re c or d 

NCTIME — Number of words in the date and 

(CTIME iNCTIME) ) — Array containing the 
that the inverse matrix file was created. 

NTITL — Number of words in the title array. 

(TJTL (NTITL) ) — Title of the inverse matrix file. 

NfiECS=Tso plus the number cf records required to 
contain the entire inverse matrix. 

(IFGRM (N8ECS) ) =0 ,0 , 1 , 1 , . . . 

NID=1 

ID ( 1) identification number cf the inverse matrix file. 

N10G= 1 

LOG (1) =FIT — Flag indicating whether the matrix fits in 
core all at once. 

NINT— 7 

INT(1)=KM — Number of rows in the inverse matrix. 


on an inverse 

time array, 
date and time 


KSPR0DUC3BIUTY0F 
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1N1(2) ; =NCPF — Maximum number cf columns in each inverse 
matrix output record. 

INI (3) =NCMTS — Number of aircraft components associated 
with, the system of equations from which the inverse matrix 
was created. 

INT (4) — INV — Method used to create the inverse. 

INT(5)=NCFC — Number cf records required to output the 
entire inverse matrix. Ibis is also referred to as the 
number of inverse matrix fragments. 

INT(6)=IDC — Flag indicating whether this matrix is an 
inverse or an 1U decc ^position matrix. A value of 0 
indicates an inverse, while a value of 1 indicates an LU 
decomposition. 

INT (7) =NCGL — Number cf columns in the inverse matrix. 

NFLT-1 

FLT(1)=XhACH — Mach number. 

2n d Rec ord 

J1=NCMFS — See above. 

J 2=7 


J3=1 

NW=7*NCMPS 


(ICKST (NCMPS) ) — Array indicating whether the 
singularity distribution associated with each aircraft 
component is constrained cr net. 

(ICT (NCMPS}) — Not used at present 

(NCFFXA (NCMFS) ) — Array indicating the number of N1 
direction constraint function fragments associated with each 
aircraft component. 


(NCFFYA (NCMPS)) — Array indicating the number of N2 
direction constraint function fragments associated witn each 
aircraft component. 


(ICVA (NCMPS) ) — Array indicating the type cr constraint 
variables used with the constraint function associated with 
each aircraft component. 

(IDA (NCMPS) ) --Identification number of the 
self-inf luence submatrix associated with each aircraft 
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component. 

• 

{NNC (NCMFS) ) — Ending column number oi the 

self-influence submatrix associated with each aircraft 
component as positioned in the entire influence matrix from 
which the inverse was computed. 

3r d an d Subsequent Records 

J1=KM — Number of ' rows in the inverse matrix. 

J2=NCPF — Number of columns in this fragment of the 
inverse matrix. 


J3=1 


NW=J1*J2 — Number of inverse matrix coefficients in this 
record. 

(A (NN) ) — Inverse matrix coefficients for this fragment. 


7.3 LU DECOMPOSITION HATBIX FILES 

Each LU decomposition matrix file contains all the 
information pertaining to the LU decomposition of a matrix 
that was assembled from cne or mere influence matrix files* 
An LU decomposition matrix is analogous to an inverse matrix 
in that it is determined without knowledge cf the constant 
vector. Hence, any number of constant vectors may be solved 
using the same LU decomposition matrix. The system of 
eguations from which this LU decomposition matrix was 
generated may have involved cne or more aircraft components. 
The file itself consists of feur records. The first two 
records are exactly the same as the first two records of the 
inverse files (Section 7.2) . The third record contains the 
permutation vector associated with the LU decomposition 
matrix which tells how the row interchanges tcok place in 
the decomposition. The fourth record contains the LU 
decomposition matrix coefficients. At present, an LU 
decomposition matrix may be computed by the SOLN code only 
if the entire matrix fits in core at one time. 

The following information is written on an LU 
decomposition file. 

1s t Rec ord 

Same as 1st record of inverse file except that IFOfiM(3) 
is 1. See Section 7.2. 

2n d Rec ord 
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See Section 7.2. 


Same as 2nd record of inverse file. 

3r d Re cord 

J1=KM — Number of rows in the LU decomposition matrix. 

J 2= 1 
J3 = 1 

NW=KN — See above. 

(IP {NW) ) "Elements of the pexmutati.cn vector. 

4t h Re cord 

— Number of rows in the LD decomposition matrix. 

J2=NCPF — Number . of columns in the IU decomposition 
matrix. 

J3=1 


NW=J1*J2 — Number of coefficients in the LD 
decomposition matrix. 

{ A {NW) ) — Coefficients of the LU decomposition matrix. 
7.4 SC BATCH PILES 

A brief description c± each scratch file is given below 
as well as an indication of the circumstances under which 
each file is used. 

INA — This file is used in a variety of ways. In 
general, it is needed in those circumstances where the 
influence matrix, inverse matrix, or decomposition matrix 
does not fit in core at one time. In addition, it is used 
to store the solution sub-vectors in the iteration procedure 
and the matrix partitions created during a PART command. 
The default logical unit number is 11. 

INAT — This file is used in tandem with the file INA. 
That is, information on these files are manipulated and 
stored by transferring the information back and forth from 
INA to INAT and vice versa. It is needed in the same 
circumstances that reguire the use of INA. The default 
logical unit number is 12. 

KST — This file is used to store the fragments of the 
influence matrix as read in from the POTFAN influence matrix 
files. This file is reguired for any cases using an RVNS 
command. The default logical unit number is 13. 
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INCED — This file is used to store the partitions of the 
entire influence matrix in these situations where the 
influence matrix does not fit in core at once and a direct 
solution procedure is desired. In addition, this file is 
used to temporarily store the right hand side vectors. 
Bence, the file is required for any cases using an RBCS 
command or where the direct solution cf a lar ger-than-cor.e 
sguare matrix is desired. The default logical unit number 
is 14. 


INE — This file is used during inversion of a matrix 
that does not fit in core at one time, individual blocks of 
the inverse are temporarily stored on INB. This file is 
also used to- manipulate the right hand side vectors and 
place them in the appropriate order. Hence, the INB file is 
required when inversion of large square matrices are 
performed and whenever the RBCS command is used. The 
default logical unit number is 15. 

INBT — This file is used in tandem with the file INB. 
That is, information on these files are uanipulated and 
stored by transferring the information back and forth. It 
is needed in the same circumstances that require the use of 
INB. The default logical unit number is 17. 

IN VR — This file is used as a temporary storage device 
when - the iteration procedure option is in effect. Under 
these circumstances, the fragments of the -influence matrix, 
are temporarily placed on this file. Otherwise, INVR is 
used as the logical unit for any POTFAN inverse files that 
are output from the SOLN code. The default logical unit 
number is 18. 
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8 T EST C ASE S 


This section describes a comprehensive series of test 
cases which were developed to thoroughly test the program. 
These test cases are very useful as an aid in familiarizing 
the user with the manner in which the SC1N code operates. 
They may also he used as a debugging aid when transferring 
the S01N code to a different computer system. 

A total of seven test cases are included in this 
section. Only two different systems of equations were used 
tc perform these test cases; however, a different method of 
solution was used in each case. 

The first system of equation's produces a sguare 
influence matrix of order ten. An illustration of this 
system of equations is shown in figure 8.0-1- A total of 
three aircraft components are associated with this system of 
equations. Therefore, the matrix is subdivided into nine 
submatrices and three constant subvectors. The constant 
subvectors contain 5, 3, and 2 elements, respectively . Each 
submatrix and constant sutvector are placed on a separate 
file. The matrix coefficients associated with each 
submatrix file were stored on two cr acre records. A total 
cf twenty-two records were used to store the entire 
influence matrix. Two sets of boundary conditions — were 
placed on each constant subvector file. The solution to the 
first set of conditions consists of the integers from one 
through ten. The solution to the second set consists of 
every other integer from two through twenty. The formulae 
used to generate the matrix coefficients and constant vector 
elements are the following: 


= (n - |i-j!) 


( 8 . 0 - 1 ) 



(1st set) (8'. 0-2 ) 
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n 

bj = 2a:t 3 (2nd set) {8.0-3) 

J=t 

The matrix produced by the above formulation is 
diagonally dominant and amenable to iteration techniques. 
This system of equations was used in test cases 1, 3, 4 and 
c. 


The second system of equations is overspecified and 
hence produces a rectangular influence matrix. In this case 
the matrix had seventeen (17) rows and ten (10) columns. 
An illustration of this system of equations is shown in 
figure 8.0-2. Three aircraft components are also associated 
with this system of equations. Therefore, the matrix is 
.subdivided into nine submatrices and three constant 
subvectors. The constant subvectors contain 7, 6 , and 4 
elements, respectively while the solution sub vectors contain 
E, 3, and' 2 elements, respectively. As with the first 
system of equations, a total cf tventy-twc ( 22 ) records were 
used to store the entire influence matrix, and two sets of 
boundary conditions were placed on each constant subvector 
file. The solutions to both sets of conditions are the same 
as those for the first system of equations. . The formulae 
used to generate the matrix coefficients and constant vector 
elements are the following: 


7 

a- = (n - Ji-yl) 


(8.0-4) 


m 

b: = a. ; j (1st set) (8.0-5) 

yr IJ 


where 



( 2 nd set) ( 8 . 0 - 6 ) 


J 


C o +C I 1 * c 2 ^ 


(8.0-7) 


n = number of rows in the matrix 
m = number of columns in the matrix 


8-2 



The coefficients c 0 , c, , c 2 are determined by satisfying in 
a least spares sense the rectangular systecr of equations 
that is generated by setting y=j at the location of the 
first and last element of each diagonal submatrix- This 
ncn-simple equation fcr y is required in order to achieve a 
system that can be solved by iteration. This system of 
equations was used for test cases 2, 6 and 7. 

The computer programs used to generate the two systems 
cf equations described above are presented in Figs. 8.0-3 
and 8.0-4. The input used to run these programs and. place 
the results on the various files is shown m Figs. 8.0-5 and 
8 . 0 - 6 . 


8. 1 TEST CASE NO. 1 


The first test case solves the first system of . 
eguatlons by direct LU decomposition. It is assumed that 
the entire influence matrix will fit in core at one time. 
This situation causes the fastest direct solution technique 
to be used to solve the system. The input for this case is 
snown in Fig. 8.1-1 and the output is shown ir. Fig. 8.1-2. 


8.2 TEST CASE NO. 2 


The second test case solves the second system — of 
equations by the direct Householder procedure since it is a 
rectangular matrix. The LU deccmpositicn tecnnigue could 
not solve this system of eguations. It is assumed that the 
entire influence matrix will not fit in core at one time, 
tut that the largest influence matrix fragment will fit in 
core. Hence, the COHB command is used to combine all of the 
influence matrix fragments into a single matrix and place 
the matrix in temporary storage. The extended Householder 
procedure for matrices toe large to fit in cere at one time 
is automatically used to sclve the system. The input for 
this case is shown in Fig. 8.2-1 and the output is shown in 
Fig. 8.2-2. 
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8.3 TEST CASE NO. 3 


The third test case solves the first system of 
equations again by direct LU decomposition. However, it 
illustrates the use of a different series c± commands that 
may be advantageous for situations, where the same influence, 
matrix is used on several different occasions in conjunction 
with different constant vectors. • That is, the matrix 
decomposition is created and stored on a tile for later use. 
This decomposition is the most time consuming portion of the 
LU decomposition procedure and is independent of the 
constant vector. Later, when the solution to a particular 
system of equations N is desired, the decomposition file and 
constant vector files are read in and only the forward 
elimination and backward substitution are performed to 
determine the -solution. It is assumed that the entire 
influence matrix fits in core at one time. In fact, the 
SCLN code at present only possesses the capability for 
creating and storing a decomposed matrix if the influence 
matrix fits in core at one time. The input for this case is 
shown in Pig. 8.3-1 and the output is shown in Fig. 8*3-2. 


8.4 TEST CASE NO. 4 


The fourth test case solves the first system of 
‘equations utilizing the block successive ever- relaxation 
iteration procedure. Some of the additional iteration 
procedure options are also exercised in order to illustrate 
some of the capabilities available with this procedure. The 
diagonal submatrices are solved by LU decomposition. The 
decomposition of one of the diagonal submatrices is assumed 
to have been computed and stored at some previous time. 
Hence, instead of being computed during the course of the 
iteration procedure, it need only be read in. In addition, 
it is assumed that an approximate solution has been 
previously obtained fox one of the aircraft components. 
Hence, at need also only be read in as an initial guess at 
the solution rather than having an initial guess computed 
for that, aircraft component during the course of the 
iteration procedure. A maximum of . ten iterations are 
allowed for convergence. A mild amount of cver-relaxation 
is used to accelerate convergence. Printouts of some 
intermediate iterations are also specified. The input for 
this case is shown in Fig. 8.4-1 and the output is shown in 
fig. 8.4-2. 
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6.5 TEST CASE HO 


The fifth . test case selves the first system of 
equations by a two step combination cx techniques. In the 
xirst step, the solutions for the first two aircraft 
components are determined as though there were no third 
aircraft component. That is, cnly the submatrices and 
constant subvectors associated with the first two aircraft 
components were read in and solved. Ther. the simultaneou 
solutions for all three aircraft components were determine 
by the block successive over-relaxation iteration procedure 
with the solutions from the first step input as initial 
guesses to the solutions for the first two aircraft 
components. The influence cn the third aircraft component 
cf the first two and vice versa is thus determined 

iteratively. This situation is analogous to the problem of 
determining the effect of wind tunnel walls on the 
characteristics of an aircraft. That is, initially a 
free-air solution is obtained fer the aircraft and then a 
solution is obtained with the aircraft in the wind tunnel 
using iteration. The input fox this case is shown m Fig. 
6.5-1 and the output is shGwn in Fig. 8.5-2. 


6.6 TEST CASE NO. 6 


The sixth test case solves the second system of 
equations in a manner similar tc that used tc perform rhe 
fourth test case. That is, the system is solved by the 
• flock successive over-relaxaticu iteraticn technique, 
however, since the second system cf eguations is rectangular 
rather than square, the diagonal submatrices are solved by 
the Householder procedure, and the inverse of one of the 
diagonal submatrices is read in instead of the decomposition 
cf the matrix. The input for this case is shown in Fig. 
6.6-1 and the output is shown in Fig. 8.6-2. 


6.7 TEST CASE NO. 7 


The seventh test case solves the second system of 
equations in a manner similar tc that used to perform the 
fifth test case. That is, the system is solved by a two 
step combination c± techniques. However, since the second 
system of eguations is rectangular rather than 


.EPaoMJdSmrnr of the 
obi <sb$M> #age is poor 


6-5 


a. rn 



sguare, the Householder procedure was used tc determine the 
solutions for the first twc aircraft components and also to 
solve the diagonal su hmatrices during the iteration 
procedure.. Ihe input for this case is shown in Fig- 8.7-1 
and the output is shown in Fig. 8.7-2. 
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APPENDIX A 


fl — STAMD AB DIZE D FI LE HANDLING P RO CEDURES FOB EOTFAN PRO G RAM S 


Standardized POSTS A N procedures and subroutines for 
opening and closing files have been developed to facilitate 
using and coding PGTFAN programs and the conversion of these 
codes to different computer systems. 


A. 1 FILE CREATION 


This section describes actions tahen before and after 
any POTFAN program attempts to write a POTFAN file. 

Prior to writing any permanent file onto a unit, all 
EOTFAN programs call a system dependent subroutine as 
fellows : 


CAXl OPENH (NT, IFTYP, ID, IB) 


If IS is not zero, then NT and ID are considered subroutine 
inputs. NT is the logical unit number on which the file 
will be written and ID is the file creation identifier, 
which should also be the primary file identification number. 
If IK is zero, then ID is net considered a subroutine input 
and NT is- only the default unit number. In this case the 
program reads in ID and NT from a card via 215 format. If 
the value of NT on the card is zero, the subroutine replaces 
NT with the default value. 

If the value of ID determined in either case is then 
still zero and if it is possible cn the computer system 
being used, the program will replace ID with the current 
number on the identification number file and also update the 
identification number file- 

In addition to NT, ID, and IS, IFTYP is also input to 
the program. IFTYP defines the type of file being created 
according to the following table: 


IFTY P TYPE OF FILE 

1 Geometry 

• 2 Boundary condition 
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3 Influence matrix 

4 Velocity matrix 

5 Solutions 

6 Velocity at force sensing location of N 1 segments 

7 Velocity at force sensing location of N2 segments 

8 Constraint function txansformatior matrix 

9 Zeta plot file 

10 Constrained influence matrix 

12 Preset solution 

1i> Inverse or decomposition of influence matrix 

,10 External velocity 

17 Surface pressures 

18 Surface velocity 


Once ID and NT have teen determined, tne program opens 
(if possible on the system being used) the file for writing 
using a file name determined from ID and IFIYP. On IB H 
systems, opening a file consists of issuing a DDEP to the 
operating system. On the INFONET UNIVAC 1108 system, an 
LQUATE command is involved. This feature eliminates the 
need for jet control cards tc handle files on those systems 
for which FCHTBAN programs can open files. 

The program then rewinds the file and writes a message 
indicating which unit has teen opened ana the value of IE 
ana IFTYF. 

After the file has teen opened and written upon, it is 
released ty calling another system dependent subroutine as 
fellows ; 


CALI E&DFIL (NT) 


This subroutine writes an end-of-file nsarx on the unit 
and (if required ty the system being used) , releases the 
unit. The subroutine also writes a message indicating that 
unit NT has been closed. 


A. 2 FILE ACCESSING 


This section describes actions tahen before and after 
any POTFAN program attempts tc read any PCTFAN file. 

Pricr to reading any permanent file from a unit all 
FCTFAN programs call a system dependent subroutine as 
follows : 
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CALL OPENS (NT , IFTYP, ID , IE) 


If IS is not zero, then NT and ID are considered 
subroutine inputs. NT is the logical unit number from which 
the file is read and ID is the file access identifier, which 
should also be the primary file identification number- If 
IB is zero, then ID is not considered a subroutine input and 
NT is cnly the default unit number. In this case, the 
program reads in ID and NT from a card via 2l5 format. If 
the value of NT on the card is zero, the subroutine replaces 
NT with the default value. 

In addition to NT, ID, and IS, IFTYP is also input tc 
the program. IFTYP defines the type of file being read 
according to the table in -the previous section. 

Once ID and NT have been determined, the program 
attempts to open the file using a file name determined from 
ID and IFTYP. The capability to open a file from a FOETRAN 
program depends on the system being used- As explained in 
the previous section, this may involve a DDEF or EQUATE 
command and can eliminate the need for job ccntrol cards to 
handle files. 

The program rewinds the file and writes a message 
indicating which unit has been opened and the value of ID 
and IFTYP. 

After control is. returned to the calling program and 
the first record of the file has been read, all POTFAN 
programs check to see if the access identifier is equal to 
the actual primary file identification number existing on 
the first record. If net equal, the program writes an 
informational diagnostic message and proceeds. This feature 
is meant to be a helpful filekeeping technique tor those 
systems that do not permit automated file ccntrol. 

After the file has been read and there is no further 
use for it, it is released by calling another system 
dependent subroutine as fellows: 


CALI FILEND (NT) 


This subroutine rewinds unit NT and (if required by the 
system being used) releases the unit. 
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APPENDIX B 


E STA NDARDI ZED . FORM AT O F PC1 EAN F ILES 


A standard format has teen developed fcr POTfAN files. 
This format is applicable' to all files except scratch files 
and plot files. This standard has teen developed, for the 
iollowing reasons; 


1. to minimize the effects of changes in one POT FAN 
segment on ether POTFAN segments; 

2. to allow a program to be developed which can list 
and/or edit the contents of any PuTFAN file; ana 

3. to promote consistency among PCTFAN programs. 


Briefly, the standardized P01FAM rile consist of one or 
mere records. The first record is called the introductory 
lecord and contains miscellaneous data including tue primary 
identification number, a title, and real, inteyer, and 
logical parameters reflecting hew the data on the remaining 
records was calculated ard/cr how it is tc be used. The 
second and subsequent records generally contain the bulk of 
the data and are called data records. The latter records 
contain one or more arrays which are always either integer 
cr floating point numbers (i.e. integer and floating point 
numbers are not mixed cn a single record). A detailed 
description is given below. 


Firs t Rec or d (Introductory Record) 


This record is created by an . unformatted write 
statement such as the following: 

WRITE (NT) NCTIEi, (CTIME (N ) , K= 1 , NCTI EE) , NTITL, 

# (TIT! (N) ,N=1,NTITi) , NRECS, (IFORM (N) , N=1,NKECS) , 

#N1D, (ID (8) , 11=1, HID) , N10G, (LOG ( R) , R = 1 , NLOG) , 

# N I NT , (INT (N) ,K=1 ,NINT) , NFIT , (FIX (N) , K= 1 , KELT) 



The values of NCTIBE, NBECS, HID, NLCG , HINT, and NFL1 
are all at least one and can vary from file tc file even for 
files of the same' type (e.g. HINT may he different on two 
different- geometry files). An explanation of these 
variables is given below; 


KCTIME, Number of words in (CTIME) 


(CTIHE) 


NTITL 

(TITL) 

NBECS 

(IFORM) 


NID 

(ID) 


KLOG 

(LOG) 

RINT 

(INT) 

NFLT 

(FIT) 


Creation time in A4 alphanumeric format. Whether 
or not this array can be filled out depends on the 
availability of a system dependent subroutine to 
compute it. This array is used' only as a 
filekeeping aid. It is printed cut whenever a 
file is created or read.? 

The number. of words in (IITL) . Generally NTITL is 
a multiple of 20. 

Alphanumeric titling information (e.g. "Delta wing 
with flaps"). This array is to be written under a 
format such as (lx,20A4/). 

The number of records (including the first) 
comprising the file. NBECS is also the number of 
words in (IFOBH) . 

An integer array indicating the kind of numbers on, 
each record. A value c£ zero implies an integer 
and a value cf one implies a floating point 
number. IFORM (1) has nc significance. 

The number cf words in (ID) 

Identification number array. IE (NID) is the 
primary file identification number. In order tc 
keep track of files IE (NID) should be unique for 
each file. This number is printed out whenever 
the file is created or read. 

The number of words in (LGG) 

An array of logical parameters 

Number of words in (INT) 

An array of integer parameters 

Number of words in (FLT) 

An array of floating point parameters. If the 
remaining data on the file is dependent on Each 
number, then FLT (1) is the Hach number. 
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Seco nd an d Sub se quen t R ec or ds {Data Records) 


The remaining records o± PCTFan tiles contain one or 
irore arrays. If the data record contains more than one 
array, then all arrays on the record must fce of the same 
type (i.e. either integers or real numbers, but not both) 
ana all arrays must have the same number cf words. The 
records also contain array dimensions (Jl, J2, and J3) and' 
the total number cf words in all arrays cn the record (Ntf) . 
hollowing are some examples of code used to create data 
records: 


N W = J1*J2*J3 

WRITE (NT) J1,J2,J3,NW, ( { ( A (I , J , K) , 1= 1 , J 1 } ,1=1,12) ,K= 1,13) 


13 = 2 

Ntf = J 1*J2*J3 

WRITE (NT) J1,J2,J3,NW, { (A (1,1) ,1=1 ,J1) ,1=1 ,12) , 
M(B (1,1) ,1=1, Jl) ,1=1,12) 


J2 = 1 
J 3 = 1 
N V» — 11 

WRITE (NT) J1,J2,J3,NW, (A (I) ,I=1,NW) 


12 = 3 
J3 = 1 
NK = 3* J 1 

WRITE (NT) J1,J2,J3,NW, ( A (I) , 1= 1 , 1 1 ) , {0 (I) ,I=1,J1) , 
#{C(1) ,1=1,11) 


Note that in the above examples all dimensions with 
multiple arrays were written with the leftmost indice 
varying most rapidly. This practice is always followe 
unless it is strictly necessary tc do otherwise. 

No matter how a data record was created, it can be read 
in by either of the following; • 

READ (NT) J1,J2, J3,NW, (A (I) ,1=1 ,NW) 

READ (NT) J1,J2,J3,NW , ( ( (A (I, J, K) ,1=1,11) , J = 1 , J2) , K= 1 , 13} 


In the termer case, the data is packed solidly into colo. 
In the latter case, some a priori knowledge of Jl, 12, and 
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J3 of their maximum allowable values must have . been 
available in order to properly dimension (A), such a priori 
knowledge is generally contained as elements cf {INT} . 

Different data records may contain data of different 
types and may have differing valued of J1 , J2, J3, and NW. 
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AFPEKDIX C 


C_ ARRA Y _ NOTATION 


A shorthand notation fcx referring tc arrays m the 
internal and external documentation or POTJFAH programs has 
teen developed. “This notation should he made clear oy the 
rollowing examples: 


(A) 

4* (H) ) 

* 

A{N) 

(A (L,J) ) 


A (I, J) 


This implies that i is an array. 

This refers to all the words in (A) from 1 
through N. 

This refers only tc the Nth word of (A). 

This refers to all the words in the doubly 
dimensioned array A for which the first index 
varies fretr 1 to I and the second from 1 tc 

J. 

This refers to the element in (A) for which 
the first index is I and the second is J. 


(A (I, J) , a- 3 , K) This refers tc the words of (A) for which the 
first index is I and the second index varies 
from 3 to K. 

(A (I,*)) This refers to these elements of (A) for 

which the first index varies from 1 to 1 and 
■ the second index varies frem 1 to some value 
which for seme reason cannot fce defined. 


S^BODVCffilLOTOFraE 


C-1 



UROGRAM 


CALLEE BY 


CALLS 


ATR 

■ INV ES , ITEM, 

RIVNS 

None 

BLKINV 

INVSS 


U IN V 

FALL 

INVRS, ITBB , 
MIN V 

MATHLT, 

None 

HOUSE 

INVSS 


Ncne 

HEEC 

INV BS 


None 

INVOUT 

INVRS 


OEENtf, FILEND, TIMES! 

INVRS 

OPRN, ITEM 


OPENW, ELKIN V , FALL, 
HOUSE, flREC, 1NVOUT, 
MATFRC, MINV, SOUT, 
TRIP, SCLVE, ATR 

ITEM 

OPRN 


ATR, FALL, .INVRS, HATCC8, 
MATPRF , HATPfiO, OPENW, . 
PART, FILEND, TIMEST, 
TRIP, GPENR, RELFIL 


TABLE 4.1-1. Subprogram Calling Structure 
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PROGRAM 


CALLED hi 


CALLS 


HATCCM 

SCLNIO, USE 

Ncne • 

BAT ML 'I 

SCLNIO 

El AT PRC, SOOT, FA^L 

1*1 ATP PE 

ITRK 

None 

MAT PEC 

IN 7 RS , ITEM, MAT ML I 

Ncne 

MINV 

ELK IN V , ISVES 

FALL, TEIP 

OPENS 

ITEM, BBCS , RCIVNS, 
RIVNS, R V NS 

Machine dependent sub- 
routines only. 

OPENW 

XNV GUT , ITEM, SCUT 

Machine dependent sub- 
rGiitines only. 

OPRN 

SCLNIO 

INVE5, ITRM, PART 

PART 

CPRN, ITRM 

Ncne 

RBCS 

SCLNIO 

RELEII, RHSCOd, OPENS 

RCINVS 

RVNS 

OPENS, SELFIr 

KELFIL 

ITEM, RBCS, RCIVNS, 
RI7NS, RVNS 

N one 


TAEL! 4.1-1 .Subprogram Calling Structure (Cont'd) 
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skHsw* 



PROGRAM 


CALLED BY 


CALLS 


EH S COM 

8BCS 

None 

El V NS 

SGLNIO 

AIR,- OE-ENR, RELFIL 

RVNS 

SCLNIO 

CPENH, RCIVNS, 
RELFIL, RVNS2 

RVNS2 

RVNS 

None 

SOLN 

None 

SCINIO , TIKEST 

SOLNIC 

SOLN 

MAICOM, MATMLT , OPEN, 
SECS, EIVNS, RVNS 

SOLVE 

INVRS 

ETIME, ETIMEF 

SOUT 

HATMLT, INVRS 

GEENW, FILEND, TIMES! 

TIMES! 

INVOOT, ITEM, SOUT, 
SOLN 

System dependent sub- 
routines only. 

THIP 

MIN V, INVRS, ITRH 

None 

TABLE 4. 1-1 

Subprogram Calling 

Structure {Concluded) 
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common 

BLOCK 
NAB E 

SIZE 

USING SUBPROGRAMS 

CONST 

10 

Nearly all 

ESC OH 

1 10 

ITEM, EECS, SCUT 

INVCOK 

70 

IN V CUT , INVES , ITEM, MATCUM, 
PA1KLT, OPEN, PART, KBCS, RClVNS, 
HHSCOH, RIVNS, RVNS, RVNS2, SOLE, 
SCLNIC , SGUI 

1NVFLG 

50 

ITEM, UCIVNS, RVHS 

NIDCOM 

864 

INVCUT, INVES, ITEM, BATflLT, RbCS 
RClVNS, 81 V NS , EVNS, EVNS2, SUj.N, 
SCLNIC, SOUT 

NOCOB 

9 

INVCUT, INVES, ITEM, tiATCOM, 
EATH1T, CEBU, PANT, RBCS, RCiVNS, 
•RIVNS, EVNS, R V NS 2 , SOLN, SULNIO 

NTCOH 

200 

ITEM, RECS, RCIVNS, SCLNIC, SOUT 


TABLE 4.3-1 Ccirrn.cn Block Usage 
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fOfirRAM 

VARIABLE 


LOGICAL UNIT DESCRIPTION 


CO fig ENT 
VALUE 


IHL 

Tempcrary storage device 

11 

KST 

Temporary storage device 

13 

INCBD 

Temporary storage device 

14 

INLT 

Temporary storage device 

12 

INB 

Temporary storage device 

15 

INBT 

Temporary storage device 

17 

INVR 

Temporary storage device cr logical 
unit number for output of inverse or 
decomposition matrix file 

18 

NTIR 

Temporary storage device 

16 

NTCP 

Printed output device in 
conversational mode 

6 

NTP 

Printed output device in batch mode 

6 

NTCR 

Conversational mode input device 

5 

NTR 

Batch mode input device 

5 

RIN 

Batch mode input device 

5 

ROUT 

Printed output device in batch mode 

6 

OUT A ( ) 

Array of logical units for 
solution output files 

User 

Input 

NBCR ( ) 

Array of logical units for 
constant vector input files 

User 
Inp ut 

ISPS ( ) 

Array of logical units for preset 
singularity input files 

User 

Input 


TABLE 4.4-1 Summary of Logical Oaits Used by SOLN Code 
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FiGUBii 1.0-1 Potential Plow Analysis System Structure 
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FIGURE 4.1-1 Equation Solver .Program Subroutine Structure 
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FIGURE .4.2-1 Equation Solver Program Logic Flow Chart (Cont’d) 
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FIGURE 4.2-1 Equation Solver Program Logic Flow Cxiart (Coat'd) 


KESPRODUCmiTY OF 

ORIGINAL PAGS IS VO 


FIGURES- 11 






N 


FIGOBE .4. 




-1- Equation Solver Program Logic Plow chart (Cont*d) 

FIGURES- 12 









frlGUBE 4.2-1 Equation Sclver Program Logic Flow Chart (Concluded) 


FIGURES- 13 





FIGUBES— 14 


AX = B 



AIRCRAFT COMI 


AIRCRAFT COMI 


AIRCRAFT COMI 


SUB-MATRIX OR SUB-VECTOR BOUNDARY 

SUB-MATRIX FRAGMENT BOUNDARY 


FIGURE 8.0-1 Illustration of the System of Equations Solved 
for Test Cases 1 „ 2 , 3, and 5; a 10x10 Square Matrix Composed 
of Three Aircraft Components and 22 Matrix Fragments. 
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SUB-MATRIX FRAGMENT BOUNDARY 


FIGU0E 8*0-2. Illustration of the System of lijuations Solved 
for Test Cases 2, 6, and 7; a 17x10 Rectangular Matrix 
Composed of Three Aircraft Components and 22 Matrix Fragments. 
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DIMENSION AftSO, 150),6f l50),TlTLf5) ,IF0RMV(10) f IDVfS) ,L0GV(12) ,INT 
iVn2),FLTV<6),lP0RMBC9!,I08f?),L0GB f 9),lNT8(3),FLTB(t),F8Crin 

nlM£tfWl^-WR^-f3rnNCPCf3)VNfFCfS,33,MpPFr3f 3,31 , NTVffWT3,‘3)> NTBCW f3“ 
1) 

LOGICAL LOGv.LOGR 
- • LOGICAL CONV, BATCH 

fOMMON /CONST/ CONV, BATCH, DllMY(tf) ,NTCP,NTP,NTcR#NTR 
DATA TITl/«HTEfiT,4H CAS.4HE FO.OHp S0.4HLN / 

--GATA— If"0tfWV/o,9*r/, tFOR M S/0,6*!/ 

DATA L0GV/F,T,«*F,T,.F/,INTV/f),5, 2,a*ft, 10#10#1» l0#2l/,FLTV/0 t ,l,,0 t 
i.O.,0,.0,/ 

DATA LOGB/«*F,T,F,T,F,F/,iNTR/0, 10. l/,FLTR/i ,0/ 

DATA Fsc/l ,0, 10 * 0 , 0 / 

DATA NCMP /1/ ,NFPC/t/,MULT/l/,MExP/3/ 

“Na^EHST'/MDaTA/NCMP, NRPC,NCPC> NFPc# NPPF,NTVNW,NTBCW,MULT,MFXP,’NS0L' 
NTR«5 
wTCR»5 

NTP#|j 

NTCP«6 

lo CONTINUE 

- -- -PE-*{)-tfT*TP*OMPT*+€NTER Nj <) m 

IFfN'.LE.O) STOP 777 
NSQ«N*N 
N80L«i 
NRPC f 1 3 «N 
NCPCfO«N 

~ NpPPt‘i‘Tll'1 3*N 

PEADfS, MDaTaI)’ 
do 2 1*1, N 
DO 2 J*1 , N 

KK*( (N.IABS(I.J) )*MULTi**MEXP 
A ( I , vl ) «KK 

2- - -CON-TtMIE 

DO 6 I»1,N 

ct.o: 

Dtl U K»i,N 

■U CT*CT + A(I,K)*FLOATfK) 

B(I)«CT 

<r -CON^rWf- : 


FIGUBE 8.0-3 Program Used tc Generate the Matrix and Constant 
Vectors for Test Cases 1,3,4, and 5. 


1 


FIGDBES- 16 


reproducibility of the 

An^/xfxT IT "DA /Itn TS POfTR. 



30 


>5 000 


M T T T t, * *5 
fiinvsj 

n(_ogv*1£ 

M r M T V a 1 2 

MPLTVSft ♦ 

TpVC 0*1 
_r^TVf 

‘f KiTVc I ?)»2i 
mJDBb? 

NtdGRSQ 

NI>>!T«e? 

MFLTP.»! 

TD«f! ) a 1 
W p U M T B ♦ 

PU^JaJ 
m v F S s 3 
"F-SChI 1 
m 3 a 1 

MPpCsBsl +3 *ns0i 
w R c 1 

nD 2o TC = 1>C MP 
x'Csl 

TMTVfg^eViRPCf If) 

. T MTVr<?)sMRPCf 1C) - 

1 k*KR+NPPC f jf)*l 
fin SO JCai.NCMP 

t F f MFPf ( JC» t Q, ' 1 NjRPF ( l f JC. TC5 BNRPC (TC) 

LCsKC+^CPC MCH 
NPFCSVal+NiFPCfJC.TC) 

MTV Mb 1 

CALL 6P£Mw r NTVW, ?,NTVN.« f Jc. IC) , 1 3 
inv ( 3 ) bntvnw r JC , IC 3 
ruTVf s J8NCPC <\jc) 
mfpciCbnfpc f nc, m 

T M T V f J 1 ) B 0 

on 30 lPsi,NFPric 

tf (Jmrpf ( TP. jr, in.GT, iMTvr t n 3 iMvr H >snrpf t if , jc, icj 

wRT Tr ( NT V*v3 NnUHT.MlMJHI.NTTrUTlTL.NPCCSV.flFO^HVfTJ.tBl.NRECSVJ.N 
1 T nv , T 1 ; V , N|_OGV, LO(?V , N T NT V , I N| T V , NFL T V , FL T V 
wo T TP ( 6 , 00 00 3 MT!Tf„. 1 , ITl.,NPf--CSV,flF'nFMVfI)*I = l # NRfcCSV) 

F0PMaT(2X, » nTVN \i ,15 ,5 M/2V ,2015) 


FIGUBE 8.0-3 


Program Used tc Generate tte Matrix ana Constant. 
Vectors icr Test cases 1 # 3 # 4, and 5 (Cont’d) 


x&BMvmam of the 

ORIGINAL PAGE IS POOR 
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WRITE (6, 6001) NlOViIGV,NLOGV,LOGV,NINTV,lNTV,NFLTV,FLTV 

6001 FORMATf2X,UIS/?X,!5,i2L5/2X, 15. 1 2I5/2X , 15 , 2X , 6E 1 2 , 5 ) 

k#kr 

DO 40 IF«i,NFPCIC 
WI»NRPP( IF, JC,IC) 

-WfrNrPCTiJC 5 ' 

L«K+N1»1 
MW*N 1 *N2 

WrITE(NTVN) N2,N1,N3,NW, f ( A ( I , J) , I*K , L) , J«KC , LC) 

WRITF(6,6002) N2,Nl f N3,NW #( (A(I,J )# IiiK,L), Jj»KC,UC) 

6002 F0RMAT(2X, »NTVN 2 » '2X , 4JS, (/2X , | OE 12, 5) ) 

‘tort'*! “ ' ■ 

4o CONTINUE 

CALL RFLFILtNTVN) 

*C»LC+1 

50 CONTINUE 

NTBC»2 

CALt--€P£NWfNTpC,2>NTRtW{IC),rr • 

TDB(2)*NTBCW(jC) 

INTBf2)»NRPC(IC) 

INTe'f S)*N80L 

WRITE (NTPC) N0UMT r N0UMI,NTITL,TlTU,NRECS8, f I FORMS ( I ) , I * 1 , NRECSB ) , N 
lT0B,lD8 # NL0GR,L0GB,NlNTe,INT8,NFLT9,FLTB 
— WRITE(6;6003) NTITL.TITL, NRECSB, (IF0RMB(l),lFl,NREeS8-) 

6003 FORMAT (2X , * nT8C 1 • /2X, I5,5A4/2X,20I5) 

WRITE (6, 6004) NIOg, iDB,NL0G8,L0GB,NINTR, INTB,NFLTB,FLTR 

6004 FORMAT f 2X, 3I5/2X, I5,U5/2X, 15, 3I5/2X, 15, 2X,11E1 0,4) 

' N t «NRPC(IC) 

DO 100 isol«i,nsol 

- &&-200- IbKRtLR - - - 

200 

WRITE (NTBC) Nl ,NDUNI #NnUMl,Nl, (B(I) , I bKR»LR) 

WR-ITf (6,6005) N1,NDUMJ,N0UMI,N1, (-8(1) , IsKRpLR) 

6005 fORMAT( 2X, »NTBC 2i/2X,4tS,(/2x,lOEl2,5)) 

NW*Nl*NVFS 

--W4E<-N-He>" NF8C*NDt tMI,N&iiMI-,NF9C7FSC- - 

too CONTINUE 

vRiLR * 1 

CALC REiFIL(NTRC) 

20 CONTINUE 

STOP 

- F*f>~ - 


FIGURE 8-0-3 Program Used tc Generate the Matrix and Constant 
Vectors for Test Cases 1,3,4, and 5 {Concluded) 
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01 MENS I OK 4r 1 5fl» 15A) #Brl50) ,TITLf5) , T FORM v f 1 0 ) , I Dv f 3 ) , LOOV f 1 2) , INT 
IV f 1 2) f ^LTyefe) r f9) # ,LOr,Hf<9) , ( 1 ) ,F8C 'U ) 

nl»6NS!0M M^Prrs) ,Mc p Cf3» 3»3> # H«PF r 3,3, 3 J , ntvmw rl, 3 ) # ktrCw^s 

1) ‘ 

0 I « tNS 1 ON. t 'M5'1).vT f S),yC(3),IT(fe) 

LOGICAL tOGV.LOGQ 
LOGICAL CONV. HaTCh 

fOMMON /CONST/ C 0«V V , 8 A TCH * DUM Y ( U ) , N Tf P , N T P # NT £<? , N T R 
DATA T ITL/'JHTtST,/!* CaS, 4HF FO,4jHR SO,'JhLN / 
nfiTA IFO»MV/0, 9* 1 /, TFORmB/O, «* i / 

RATA L06V / F,T,«*F,T,ff/,J«gTV/0#5 < 2.A*fl,io»lO.|#10»2t/#FlTV/0,,l,,O t 

Ij0 9 o0 a> 0 e / 

RATA i OS0/<l*F # T,F, T,F # F/ f jNtP/O, t 0, J / ,F t Tg/J ,0/ 
n*TA F $ C / J , 0 , 10*0,0/ 

KAMEl. IST/MOaTa/NCMP, MPPC»NpPc,N 4 FPc,NRPF,NrvNw,N7 0CW,MnLT,MFXP,MSnL 

NaMEI IST/STZE/NR.wr 

MTHSS 

wTC.Ra'S 

V T P a a 

MTCPsft 

<0 rONTTNUE 

*?EADrS,$I Z fc) 
mS0snN*NC 
m«50i st 

OF. ADfSjMDftT at 

.TACUMsf 

?ACUMsi‘ 

msq 

00 1 Isl.NCMP 
M hM « p 

IIoMe.! 5SIACUM 
II fM)sNCPCr I) + lACHM»t 
TACUmstACUM + NCPC f y 1 

VI ( i 1 sjACUM 

vy fM) B N»PCf r)+.I4CljM-l 

.JACUMs JACUM + NKPc ( I I 
3 rONTtvljE 

[ MAXs?*NrMP 

Mgn 


FIGURE 8.0-4 Program Used tc Generate the Matrix ana Constant 
Vectors fcr Test Cases 2,6, ana 7 
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UU. I net *4 

ki*k»i 

• 90-5 

M*N + ? 

0(M)aIICU**Kt * 

f COWTT^tJg- - - -• - — ' 

t continue 

NS*i 

C4LU H0U9E'tLMAY,5,0,NS, Yl,YCl 
00 2 J*i,NC 

Y*YCf n+YCr 2>*,r+Ycr3 j*j* j 

00 - 2 I«i v N*~- 

Af I,J)s{ INR-ABScl.Y) ) *MULT , s **MEVP 
3 ' CONTINUE 

00 fc 1*1 , NR 
CT*0l 

. 00 <* K*i,NC 
- CT*CTYftfl-,lO*FLO-ATrK5 

6 CONTINUE 

WTTTLS5 
WIDV*J 
NL0GV*ig 

— xriNfv*ig - 

WFlT v *ft 

I0V { 2 ) *2 
lNTVf33*j 
f N7V f 12)*2l 

. 

wtOORSR 
NINTR81 
- ■ ■ MFLT*#i 

NDUMlal 

-lYUM-Irt" 

MVF5S3 

NFSColl 

■ us»r 

MREC8B*jV2*NSOL. 

KR*J 

— pcrro Tcw t rN cnR- 



riGHEES-20 - 



Kf!Ot 

TNTvr2)BMF>Pcan 

TNTV f t ) ) sNRPC ( ic ) 

LRoKR + fJRPC ( 10-1 
00 50 JCal,^CM° 

TFfNFPC ( JC, IC) ,E0, 1) Mp PFf 1 » JC,IC)aN p PCf IC) 

LCsKC + MCPC'f JC)-t 
MRECSVsUMPPCf .JC, TC) 

WTVM* 1 

C ALL OPEMwrNT V M,S,NTvN^fJC»lC)#n 
TOVf ?)fiNTvNwr JC f TO 
tWTVfSls*C°Cf JC) 

KlFPCTCaNFPrr JC»IC) 

T mt v r 1 1 )*o 

no 3n TFst.NFPCjr 

To IP r NPPff ( Tf» JC» IC ) ,GT . pt Vf 1 n ) INTV( 1 n sNRPf f IF» JC, jC) 

. WR I TF(MTVN) MntiMl f MnllMT , NT ITL, TIT|„^ NR£oSV, f IFDRHVf T ) , Isl , MP&CSV1 

tTnv,TDV # NtOGv/,LOnv,NTNTv,INTv,MpLTv,FLTv 
WRlTP(6,feonol W TTTL,TITL,NI>ECSV, f iformv ( I) , I* \ , n«ECS vj 
noOO FORM T ( ?Xr * NT^m 3 • /2 X , 1 5 , 5 A4/<?y , 20 T 5 T 

wRITp 60 01 ) NJtnv,mv f mlo.gv,logv/ ) I n, tv, iwtv.npltv.fitv 
*0 0 1 POO.MaTt2X,4I5/2/,!*5, 1 2L5 , T5 , 1 215/2* , J5 , 2x , for 1 2, 5 ) 

Ks*R 

n(] HO fFsl^MPPCjC 
NlsNRPFf TF, JC, IC) 

M2»NCPf f JCl 

( S i< 4 M 1 » 1 

wRITp(NTVNT N2 ( Nl,N3,NW,((A(J,J),IsK,LT,JcKCfLCT 
wrITF (6,6002) m 2* Nt>»3fNHlf r(A(T f J) f T«K t ( > ),JskC*i.CT 
*002 format (2*, »ntvn 2 , ygx.ais, r/2x, iOEta.5) j 
■< c L + 1 

no oONTtNijE 

CALL REIFIL f M T V iv T 
KC=lC+1 

Sf CONTINUE 

vTBC*<? 

CALL OPENrt f NTRC,2 # MTBCFfIC)»n 
TD8(?)*NTBCwf IC) 

THT0f2) B N«PCf IC) 

INTBf 3)sNsni v 

WfilTPUMToC) NO'JH T , NpUM r, NT TTL,T I TL,NRFC SR, ( IFQRMfi ( I ) , I a 1 , M»fcCSR) 


FIGUdE b - 0-4 Prograni Used tc Generate the Matrix ana Constant 
Vectors for Test Cases 2, fa, and 7 (Cont’d) 
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1 IDB,IDB,NLOGB,LOSB.W!NTB, INTB,NFLTB ( FLTB 
WRTTF ( 6 , 6003) NTlTL»TtTt»NRECSB, f jFoRMB f I ) » !■ 1 , MRECsB ) 
6005 FORMAT (2X, »NTBC 1 I /2X , 15, 5 A4/2X , 20 155 

WRITE (6,600a) N iDe# XDB* NLOGB , UOGB, N INTB , INTB , NFLTR „ F-LTP 
60 na BORMAT(2X,3I5/2x,i5,Pl5/2x.l5,3l5/2x, l5,2x,ilElO,4) 
Nl«MRPCtIC) 
no iflO ISQL«l f MSOL 
DO 200 IbKR, LR 
.200 B<n«ISOL*B(!) 

WR!TEtMTBC5 Nl f NDUMj,NpUMI,Nl,(B(!)#I»KR,LR) 

WRITE (6, 6005) M1,NDUMI,NDUMI,N1, (B(I),IoKR,UR) 

-60^5- *MTBC 2 ♦ /** , 4 l5yf/fX> iHl 2*5 ) ) 

NW«N)*NVFs 

WRITE (NTbC) NF8C,N0UMI,NDUMI,NF8C,F8C 

i oo CONTINUE 

vRcLPi 1 

CALI RELFIL(NTBC) 

-■^0" tOWjWUf- -■ 

STOP 

End 


FIGURE 8.0-4 Program Used to Generate the Matrix and Constant 
Vectors for lest Cases 2,6, and 7 (Concluded) 
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SEPBODUCIBILTry OF TBS 

original page is poos 



1 fl 

SWD * T A *C «Rs * , fafiPceS , 1 '■C PC*b , 3 , PCs S , i , 3 , ? , 2 f 5, 2 , 2 , «? , 
NRPFs? , a , 1 , i , 1 , 1 , 1 , 1 , -5 , 2 , 1 , 0 , 1 , ? , 0 , i , 1 , 1 , 1 , 1 , 0 , 1 , 1 , 0 , ! , 1 , 0 , 
vTVfv*!2 7n«?9 , 7nt>u, 7 Oh 1 , l{)bi, 70b J, 7*>h.'J , 70bS, 7l>he>, 70b 7, 

NTaf w: 7 Of 8 , 7 ft 9 , 7 07 0 » V1 'L Ts 1 , »F XP= 7, i*S0| s i? * E *n 


FIGURE 8.0-5 Input Used to Run SIEST Eroyram 


f IGURES-23 



$ S I ZF vRa 1 7 j n C = 1 0 “0 

$MD4TA NCMPs^,NRPC = ’^6* !j / N C p C = 5,3,5,t\'F p C = 3,-3»3,?,2|3,?»oc / 

NRPFa4,2, 1 t S,l, 3#2,3#2#S?# ( >, 1# 5, 0,3, 2, 1 # 3, 1 1 r***' 

nTVNW*70?7 , 7o 78 #7079, 7080, 7081, 7082,7083, 70 8-4,7085, 

NTBCW=7 0B6, 70 87.#7088#MULTS1 ,MEXP = 7,^SOL»2 SE 1 ^ 


1 FIGURE 8.0-6 Input Used to Run SIESIH Program 


’ sS®® ss ‘ 



TRUf 

SAMPt .-- f AST 1 FOP THf S01.N r.Or>t‘. 
SOll'TfO^- ttY (* I K£f. T 
tU Otr O^PfiST TION HEThOn, 

HATCH M * 1 1 j ft > 

THREE A/RC^aFT cnupUNfcNf S, 

two C. nsTani vk; Tors', 

CNTl 

1 

Pf AO INHiiENCE SUB-MATRIX FILES' 

VNWP 

?0*S<) 

7 060 
7(>6l 

7064 
706u 

7065 

7066 

7067 

SJNtt T*E EMTlPfc *«TRlx FITS IN 
rnRt, Cr-MH COMMAND IS NfiT 
NEEDfp, 

READ Uip.STAMT SUB*VtCTnK FILES. 

BCPt 

7068 

- 70e<J • 

7070 

SOLVE MY IU DECOMPOSITION, 

AND O^FiNf the SOLUTION 
SUB^VtCtOH OUTPUT FILES' 

SOL V 

o r- 

707! 

7r>7<? 

7073 

0 

ALL Dun*- , 

STff* 


FlQUiiE 0.1-1 Input for SCIN Cede for lest Case No. 1 
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POTFAN EQUATION SOLVING PROGRAM fSOLN). VERSION i.i 
DYNAMIC MEMORY * toooo 


We » — ^>*5/76 htotm 


+i Sample case i for the soln code' 
smrtrmm- B Y di r ec t 

♦ I LU. DFp OHPOSJ T'lON M£TH0p 

+i batch mode; * 

'*T A E AI R CRAFT com p on e nts; 

♦I TMO CONSTANT vectors; 

♦I CNTL 


Nf+tP* l- — mPA-T^ 0- 

♦ » RfAD INFLUENCE SUB. MATRIX FIiES' 

+1 VNRfc u * 


- FILE 7fl59 VNaPNe 


• ^AS been OPENE D for r eading ON UN !T- 


CREATION f IMfc.,07/07/76 00*0«*5a 

UNIT i REWOUND AND RELEASED 

£!c*rInt°;!(L J r PNC ~ ~ WA-S-ftEtN -OPTNED- FOR-REaD I N fi- im~ UNIT 
CREATION TIME, 07/07/76 00*0*150 

UNIT 1 REWOUND AND RELEASED 
-F-ILf 7o6l»VN»PNC H as B£g - N - 




- OPENEB FO R R Ea P INS ON UNIT — ^ 


CREATION T i ME ,07/07/76 00*08*5*1 

UNIT 1 REWOUND AND RELEASED 
“PTLS 7 062i 

CREATION TIME, 07/07/76 00*08*50 

UNIT i REWOUND ANO RELEASED 

—PILE TObl, VN*PWc~ H y S BE E N 

CREATION TIME, 07/07/76 00*08*50 

UNIT l REWOUND and released 
-pile 7oto t vn«.pnc 
CREATION TjME»07/07/76 00*08150 

UNIT i REWOUND AND RELEASED 

—FILL 7o65,VN»FNC H aS- 8EEN - 

CREATiON TjmE*07/07/76 00*08*50 

UNIT J REWOUND AND RElEa8E0 

-Pitt- “7066* VN*pNC — M a 3 BEEN OPENED -P6ft~ Rt~at> FNS ON UNIT 

CREATION TjmE*07/0T/76 00*08*50 


“1 


- o pe ned f -or - r eading on u mt — r 


~SN UNIT - r 


O P ENED FO R- R tg- gpING QN -UN pT y 


FIGURE 8.1-2 Output for SCLN Code for Test Case No. 1 
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♦ I 

♦ I 

♦ I 
+ J 

♦ : 


upmi i Ktwriywo ANp RELEASED 

rpt a rtn^t VN*RNt MAS BEEN- OPENtO FOR *EA0lttG ON UN J T 

CREATION T In Eeo7/07/7t> 04*0S»S« 

UNIT j RhWrUJND AND RELEASED 

siNrf the entire matrix fits in 

CURf , COMB COMMAND IS NOT 
NE E l “EO, 

REaU tOWSTiNT sUH»VECfOR Hlfs' 

OCRfc 

HaS BE ’ EN 0P tNED FOR READING ON UNIT 
Creation f i Mf, 0 7/07/76 oai08l*5a 

UNIT ? REWOUND and RELEASED 

rlulr 7 ^ 9 l BCmPHC H * S BEEN OPENED FOR -READING ON UNIT 

CREATjOf TTf*Es07/O7/7M 04t08»Sa 

UNIT ? Wf wi'UND AND RELEASED 

rBt*tfSi°: Pr r PNC HaS HtEN ° PENED F (]H RE aD f.Nfi ON UNIT 

CREATION I I^F*G7/t>7/7<? 04:08l*S4 

UNIT ? REwmUmD ANp released 
►: SOI VE BY Lu DECOMPOSITION, 

♦ J A *'D prF'JNt THE SOLUTION 

♦t sun. VECTOR OUTPUT FILES 

+ j S0| V 

methods 

TPS*. o 

TPR I u r - 

TOUT A, OUT A of , js-f , NfMPj « 

7 0 71 r 

7072 ,-i 

70 7 A 

ENTER SUB, wArT «FRAr,S2 22 i NV* „2 

MATRIX LOADED- FROM UNIT ji USING ?? RFaDSand Mb SKIPS 
NO PaRT I T I ON'! ng RtftUIRfO 5 

EXIT PART 


ENTER su«' TRIP* Na 
FXIT TRIP* 


10 


output 

SOLUTION 1 IMF rSFr' ) 


,019 


SOLUTION V F r T 0 R NO* 1 


FIGURE 8.1 2 Output for £CLN Code for Test Case No. 1 (Cont'd) 
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file 707i f so.PNc/ii8$ has been opened for writing on unit 5 

CREATION time s 07/22/7 6 UtQ6l23 


UNIT NO, 3 

component NO, 1 SOIUTtONSi 

,iooooooE+Ci t eoo oo ooE+oi ^- Soooo o oriot ;ao fr ooooF*o t~ - ■ ^oodoodf+o 

unit 5 ENDFILEO AND RELEASED 

FILE 7072,SO-PNC/LI8* HaS BEEN OPENED FOR WRITING ON UNIT J 
— ORtA-T TON- T IME--* o7 722/76 1 t~T 06 1 2U 


UNIT NO, 

component 


3 


-2 — SOLUTIONS, 


UNIT 


^OOOOOOF+O! 
ENDFILEO AND 


a 8000000E + 0 J 


,7OOOOO0F*Ot 
RELEASED 

HhAS BegN OPE^ED f UR- W*m NS ON~tm f- 3 

CREATION TIME s 07/22/76 Ul0«l24 


UNIT -N O, 3 

COMPONENT NO', 3 SOLUTjONSf 

,R000001E + 0! , 1 OOOOOOE + 02 


— tHtf-T"— ; 


OUTPUT 

SOLUT I ON TIME (SEC T .? .TCTt 




FILE 7071 ,50»PNC/t t 8$ HaS BEEN OP g Wgtr-p-O R WRI T ING ON ' UNI T- ~%~ 

UNIT NO, 3 

COMPONENT-Nt'r. \ SOLUT I ONS, 

, 2000000E +01 ,«OOOOOOFF<H ISOOOOOOEiOt ’SOOOOOmOl 9999999F+0 

UNIT 3 ENDFI|ED AND RELEASED 

FILE Tp- 72 t SU.PN C /L IB s HaS BEEN OPEN E D F O R — WRIT ING W~UNtT 3 

UNIT NO, 3 

C O MPONE N T N o, i — SO LUt I tW37~ — 

• 1200000E+02 ,1«OOOOOFa02 , 1 600OOQE+O2 


FXGURE 8,1-2 Output for SCI 8 Code for Test Case 80 . 1 (Cant’d) 
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UNIT 3 ENDUED AND pleased 

FJL? ft) f% ,-80-PNC/t TB* MAS SEEN OPfNEO tO« WHITING ON UN 1 1 3 

■UN'i T NO, 3 

.COMPONENT mu' ? ' SOLUl TONSt 

,.l«'0.00.00f. + 02 * iOOOOOOF +0? 

UNIT . 3 r^oT- IUO A Nf) MELEaSEO 

♦ t ' ALL ’fTONf , . 

♦ t STnP 
STfiP 777 


t luUHE 8.1-2 Output for 8CLN Code for lest Case No. 1 (Concluded) 
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FIGOHS 8. 


TRUfr 


sample case 2 for the soln code. 
R ECU H SUU R M A TR lX -a-m-&-? 


SOLUTION R y DIRECT 
HOUSEHOLDER method' 

frft TCH - MCDg-; 

three aircraft components. 

cntl 

s 

REaD INFLUENCE SUB.MaTRIX FILES'. 

vnre 

TOT? 


7078 

7079 
■w r 

7081 

7082 


708a 

7085 

AS S U M E THa T T HE M A TR IX Wt " 

NOT FIT IN CORp ALL AT ONCf, 


comb 


R EAD C O NS T ANT — SttB^ Vt C TU R- f ILFS, 

BCRE 

708fe 

7 0 6 7 

7088 

SOLVE BY HOUSEHOLDER METHOD# 


SoLV 


8U»»VECT0R OUTPUT FILES. 


r 

707a 
7075 
- 7 V ti r 


3 TOP 


AIL DUNE’. 


2-1 Input for SOLN Code for Test Case No- 2 


FIGUHES-30 


REPRODUCIBILITY OF THE 

HPIfitWAT. PAflW T<2 TWfo 


PQTPaH tOiiATION SOL V I Np PROGRAM (SOlN), VERSION 1,1 

miesisisdiiKttiitsssiiitttefiBiastssBsmssmit 

DYNAMIC Ml moRy a 10000 

Ti«r s o?/?2/7s iijo7i?i 

+ ? Samei.f fASf ? f<3» the sol n r or>E. 

tf RECTANGULAR MATRIX f 17 x10) 

♦ : Si‘L '*! } ON RY OIR6-CT 

fj H'U)''E HOLDER METHOD, 

*: H.* Tt H MODE , 

♦ I I»-REE AlPr.RafT COMPONENTS, 

+ 5 f.NH 

NCMPa n HERAT* 0 

+ 1 Rf AD INFLUENCE SIJU-MaTRI* FILES. 

♦ t VNR1 

FILE ffi?7«VM»PNt HAS HIE « OPENED FOR R£ AD I N(- ON tjN-J T j 

CREATjON t «07/07/7h 04109101 

UNIT i REWOUND AND RELEASED 

FILE 70?fi. Vh.PNc w A ij been opened FOR REaDINO On UNIT i 

CREATION Ti^E 8 07/07/7t> 0«»09|0l 

UNIT t REWOUND ANn RELEASED 

FILE 70/9,Vm#pnc has rEEn qPenEO FOR REaDInO Ott UN i T l 

CREATION TI *E *07/07/76 G<J*ORlOl 

UNIT 1 REWOUND AND RELEASED 

FILE 70b0,VN.PNC HaS BEEN OPENED FOR READING ON UNIT j 

TREATS T J,.F*07/O7/76 0«:09»Ol 

UNIT } Rf WOUND AND RELEASED 

FUF 70*1 , VN*»pNC HAS HfcEN OPENED FOR REaDIWO ON UNIT f 

CREATION T T oFs07/07/7h 0U*09l0l 

UNIT 1 REWOUND AND RELEaSEU 

FILE 7fl«2,Vf«*PNC w 4S HtEN OPfcNED FOR REaDING ON UNI f j 

CREATION 1 ToE*07/07/76 0 a X 0 0 1 0 I 

UNIT j REWOUND AND RELEASED 

FILE 7a*^,VM*PNr h A S rtfeEN OPENED FOR READING ON UNIT { 

CREATION J J»t:*07/07/76 - OGlOPJOl 

UNIT 1 REWOUND AND RELEASED 

FILE 7fiS«,VN-PNf HAS SEEN OPENED FOR Rfc ADiNO ON tjN | T j 

CREATION T TMf 8 07/07>7S 0fl|09|0l 


FIGURE 8.2-2 Output tor SCLN Code foe Test Case No. 2 


FIGURES-31 



unit 

-HHr fr 


t REWOUND and RELEASED 

"H y S- -&fc g N- O PgN€<f -f^R~F? £ - atrt - N6 W 'ttNl-T— f 


CREATION TIME, 07 / 07 /To 0 <l* 0 Rl 0 l 

UNIT l REWOUND A NO REUaSED 
"t| ASSUME THA-T-THE MaTRIXWTU. 

+1 nut FIT in cor f all at once'. 

+ 1 COMB 

■n r h h? ■c ons t ant sug*VEC TaR - file s , 

+* HCRE 

FlLg 7.06b ,8C«PNC 


UNIT 

FILE 


2 REWOUND AND 
70#7,RC*PNC 


HAS BEEN OPfNEn FOR READING UN 

fmtmfrt 

released 

has been opened for reading on 

IWTORtOT 

UNIT t REWOUND A«D RELEASED 

FILE 70»8iBC-PNC Has SEEN opened FOR READING ON 

"CREA-tt-gN fTHEr07'^tr76r ORt&RlOl 

UNIT B REWOUND AND RELEASED 
♦1 SOLVE BY HOUSEHOLDER METHOD, 

*1 

♦ I 

M SOLV 

TPS* 

IPRINT* 


UNIT g 
UNIT ? 
UNIT 2 


A*0 D E FI N E " T H E SD tt l T t tm - 
8IIR.V F CT0R output files; 


0 


C OUT A , ou T AN T > T g y j N C HP )~ m 

707« ,-i 

7o7B <\ 

707 h r. 


enter sub, p art #frags* h inv» .t 

MATRIX LOaDFD FROM unit 5 3 USING 22 READS AND {76 SKIPS 

Wfr-P*» Tm OKtws REftUtftEtr " 

EXIT PART 


EXIT HOUSr 


-v? fm*8 


^solutions 


tt utpu t 

SOLUTION TIME f sec; > , 03 « 


FIGCfEE a. 2-2 Output for SOLH Code for Test Case No, 2 (Cont'd) 

ftEJPEODOCIBlUTY OF THB 
otu«nal PASS IB POOR 


fIGURES-32 


SOLUTION vf'r TOR NO'. 1 


FILE 707/( ,SO.pN*c/LlBS H A S BEEN OPENED ¥ OH WRITING ON UNIT 3 

CREATION T I HR a 07/22/76 U?07|<lfc ' 

UNIT NO, 3 

COMMENT NO, 1 SOUtTTONSj 

.1O0O000F+01 , 2000000F.+01 ,3oOOOOOE*OI f «oooooie+oi 

UNIT 3 EMTHLED A W D RELEASED 

FILE’ 70 7 V SD*PNC/L!B1? *aS BEEN OPENED FOR WRITING ON UN I T 3 

CRFATION ?P"L s 0 7/22/70 Ul07|47 

ONM NO, 3 

component NO, i SOU 1 TTON 5 1 

,f>000000M0j , 7000000? * 0 1 ,«ooooooe+oi 

UNIT 3 E fc -D» ILEO AND RELEASED 

FILE 70 7#>,SO-pNc/LIB3 HAS BEEN OPENED FOR WRITING ON UNIT 3 

CREATION ? 1 1 >t = 0 7/22/76 l|l07t«R 

UNIT NO, 3 

COMPONENT NO, 3 SOLDI TONS* 

. , Q !>00000F + 01 «lQo0000f+02 

UNIT 3 fc'pf TIED ANo RELEASED 


SOLUTION vf f TnR NO'. 2 


FILE 70 7« , SO-PNC /L 1 8 3 HaS BEEN OPENED FUR WRITING ON UNIT 5 

UNIT *m, 3 

COMPONENT NO, 1 SOLUTIONS* 

• ?(’0<>0OtF+0J , ROOOOOOF +01 ,6000000? +01 ,8000001? +01 

UNIT 3 E^DMIED ANO RELEASED 

FILE 707‘i > S"*PNC/LIH$ HaS BEEN OPENED FOR WRITING ON UNIT 3 

tJNT T NO, 3 --- 

r opponent NO, 2 SOLUTIONS* 

, t^OOOoOE + 02 • 1 u QOOOOF + 02 ,lbOOOQOE+02 

UNIT 3 fc*nHt ED A Np PeLEASEp 

FILE 70 76, S! -PNC/LIBT H A S BEEN OPENED FOR WRITING ON UNIT % 


flGUHE 8.2-2 Output tor SCLN Code for Test Case No. 2 (Cont* 


, 4999999? +01 


,9999999F+0 I 


FIGUIiES-33 



UNIT NO, J 

co mp onent no; i — s otw t m si 

»i^oooooe 40 ? .aoooooor+o? 

UNIT l EndHlEO and released 

41 Af t DONE , - ■ * “ 

♦| STOP 
STOP 777 


FIGURE 8.2-2 Output tor SCLN Code for lest Case No. 2 (Concluded) 


\fIGUBES-34 



Sa^'U L A&E i FOR SOLN cODfc, 
tSOLUTJO-v «Y CRPAtING Mat*!* 
DecUMPOJilTlON FItE, REAOING IT 
BACK IN, f Hr N PERFORMING FOR^aRr 
fet 1 MI f IrtM AND BACKWARD 
SUBSTITUTION flY MUlT COMMAND' 
batch m,>df, 

ThRfcF Aff?Cf?AFf COMPONENTS, 

CNTI 


* 

hM« INFLUENCE MATRIX FILES, 

VN«P 
70 SB 
?06ft 
70ht 
706? 

7063 

7()6« 

70 bS 
7 06b 
7067 

SINCfc EmT|Pf MaTRIy FITS JM CORE , 
COMB C 0 ^ M 4 N [> IS NOT NEEDED'. 

perform matrix decomposition, and 

OFF INF fht DECOMPOSITION OUTPUT 

F ItF r.hv ITS UTt+'l ' 

I NVF 

0 70RM 1 

10 

DECUMPO^TTION MATRIX F Or $A«PU CASE 
USER fOHCO STOP HERE AND PERFORM 
RtMAlN!«fi COMMANDS IaTE/H, OR HE 
may Continue’. 

beau (O-.'STANT SUB*VfcCTOR FILES; 

hcrE 

?nt>H 

7069 

loin 

RFAO t«Ff OMP0SITI0N m a T R Ix FILE 
Hack jn' 

I VNR 
7098 

PfPFOPM FORWARD FLTMINaTION AND 

SUBSTITUTION, a NO DEFINE 


FIGURE 8.3-1 Input for SGIN Cede for lest Case Ko. J 


E1GDRES-35 


REPRODUCIBILITY OF THE 
08I4KNAL PAGE IS POOR 



SOLUTION SUb*VECTOP OUTPUT FILPS'. 


MULT 

-TW 

7078 

7079 


— Ail- OWE , 

STOP 


FIGURE 8.3-1 Input for SGLN Code for Test Case No. 


3 (Concluded) 


FIGDRES-36 



POTMN EQUATION SOLVING PROGRAM f sot n i , • version 1,1 

*tcaa**a*e«#3r**a*i*s*8**«as**saa*ej***s5t*ssss»**B*S3 

dynamic memory s ioooo 

Tiff - 0?/2?/7h 1 1 I 0W j OS 

4* SaMPU CA3F i FOR soln code, 

4t SOLUTION &Y CREATING MATRIX 

4} Of COMPOSITION file, READING IT 

4j 8-'-C“ IN, THEN PERFORMING FORWARD 

4-j fcll^tWAtlON ANf) BACKWARD 

4* SUBSTITUTION BY MULT COMMAND. 

4! B>. TCP MODE . 

♦ 5 Tr-REE AIRCRAFT COMPONENTS. 

+ 1 C N T L 

NcMPs \ TIERATb 0 

4 ; Rfci*' iNFLUFNCE MaT^T* FILES, 

4 » V M R t 

FILE TOBR.VN-PN C H A S been OPENfO FOR READING ON UNIT \ 

CRfeATfON I tmE S07/07/7P 0ajO«lS4 

UNIT t REwnUND AND RELEASED 

FILE 7060 , VN-pNf. HaS BEEN OPENED FOR READING ON UNIT \ 

creation HME>or/^ f/ft> o«st>8ts« 

UNIT 1 Rf WntiND aNO RELEASED 

Ell E 7oPi,VN.PNC Has BEEN opened FOR READING on UNIT t 

CREATION TlMEsQ7/07/7f> , 0<J;0«tS<J 

UNIT 1 REWOUND AND RHEaSEO 

F,||F 7oop, VH-PNC H A S BEEN OPENED FOR. READING ON UNJT j 

CREATION ■n=€*'07/«7/7fe 04:0«IS# 

UNTT 1 RLWilUND AND RELEASED 

FILE 7QPS.VN.PNC HaS BEEN OPfNED FUR READING ON UNIT j 

CRFATjOh T jmF * 07 /07/76 oa*ofti5« 

UNIT t RFWr.IjNO AND RELEASED 

FILE 706ii.VN.PNC H ft S been OPENED FOR READING on UNIT 1 

CREATION TTMErtE7707A7P imoatStt 

UNIT 1 REWOUND AND RELEASED 

Fill 706S.VN.pNC Has BEEN OPENED for READING UN UNIT t 

CRE A I j ON f f mF*Q 7/07/76 O<i:08t*jiJ 
UNIT t REWOUND ANp RELEASED . 


FIGURE 8.3-2 Output for SC.LN Code for Test Case No. 3 

FIGUKBS-3 7 gg^^BW® 



FILE 7066.VN-PNC H A S BEEN OPENED FOR READING ON UNIT 

C R EATION TlMCBO T V O y/Tfr 6 9 * 08 * 94 

UNIT l REWOUND AND RELEASED 

file 7op7.vn.pnc has been opened for reading ON unit 



A /I m . A .ft . „ 

94100*54 


1 


i 


UNIT 1 rewound and released 
m since entire matrix fits in core, 

M CO H 8 CO M MAND IS NOT NEEDED'. * - 

tt perform matrix decomposition, and 

♦I DEFINE THE DECOMPOSITION OUTPUT 

Fftr A - Np iT ST^tr; - — •• - 

+ i tNve 

METHODS fi NTINVWs 7oR5 IDC» \ 

— ENTER sod, P ART s FRAgSb Z2 r WV - »~« ry : : y 

MATRIX LOADED from UNIT li USING Z2 READS AND \ 7b SKIPS 

no partitioning required 


EXIT -F-*R-T 


**+ matrix inverse demanded 


ENTER 8UB' TRIP, N* 
EXIT TRIP 


10 


FILE 709% IN-PNC/LIBS 
Creation time s 07 / 22 / 7 b 
N- TTTL b Hi 


-t i-Mr-rsEr. > " .or# — 

has been opened for writing on unit ts 

11I08I1B 


<TITL(I),T»,,NTJTL|« DECOMPOSITION MATRIX FOR SAMPLE CASE 
UNIT JB ENDFILEp AND_RET-EASeo 

_l_l USE^R pOUtO 8 'T'GP HE R E"~ A NO "PE'ftF ORM : 

♦i remaining commands later, or he. 

*1 May CONTINUE, 

-n RtAO C O N STA N T SU B . VECTO R FlLE 8 r t — — 


+ i bcre 

file 70b8,BC*PNC HAS BEgN OPeNEd FOR ReAqING ON UNIT 

UNIT 2 REWOUND AND RELEaSEO 

FILE 7069,BC*PNC H a $ BEEN OPENED FOR REaDtNG-ON UNIT 

- C ^E- A - TfO N- Tf --ME « 077t ? 77 7 P fr R -T-Ofrf^-R 

UNIT 2 REWOUND AND RELEASED 


FILE 7070,BC*PNC 


has been opened for reading on unit ^ 


UNIT ? REWOUND AND R£lE*8EO 


FIGURE 8.3-2 Output for SOLN Code for Test Case No. 3 (Cont'd) 


FIGORES -38 



rfao 


DECOMPOSITION 

1*. 


MAtWIX FILE 


Has HEEN oPp.NEO 
MtOAMB 


M 
*1 

FI }VNK 

FItfc 70RS.lN.pNr 
CREATION fj*E - 07/22/76 
MACH NO, s ,J23« 

UNIT n Rt WOUND AND RELEASED 

PERFfmfT FO*W&R0 elimination and 
backward substitution, and define 

solution 5Up-VpCT0R OUTPUT FTLfS'. 


F OR READING on unit H 


♦ r 

♦ I 

.♦ I 

Ft MUl 1 
(0»T At I ) , 
7077 
?+>7« 

7079 


■in AuTf ! ) , 1st f NCMPjs 


SOLUTION '/EfToR NO'. 1 


FILE 70 77.Su.pNC/LlBS HaS BEEN OPENED FOR WRITING ON UNIT S 
CREATIn« Tff-E s 07/22/76 11lfiB-|23 


UNIT NO, % 

C OMPONFNT NO, \ SOLUTIONS t 

, !tjflOOOOt>©1 .^OOfrOOOF+Ol : ,3000000 Cf01 'AOOOOOOEfOI .SOOOOOOr + Ot 
unit t endfiiEo and released * e 

FILE 70 78.sn.PNC/lr8ii HAS- BEEN OPENED FOR WRITING ON UNJT y 
CREATION TIME « 07/22/76 

UNIT NO, 1 

- component no,- 2 •&tttuTrflNSi 

*60000001 + 01 ,7000000f + 01 , 8000000F *0 1 

UNIT y ENDEILFD AND RELEASED 

FILf . 7079 4 SO-pNc/UBt HaS BEEN OPENED FOR WRITING ON UNJT t 

CREATION Tjrt « 07/22/76 11106(20 

UNIT NO,- •=§ - - - 

COMPONENT nO, J SOLUTIONS! 

,90000017+01 . I 000000E ♦ 02 

»nit * e^df it ed and released 


FIGURE 8.3-2 Output for SC1N Code for Test Case No. J (Cont*d) • 


FIGURES— 39 



SOLUTION VErTOR NO'. 2 

■wnewmmmmmmm mikmwmm-mmmmmm 


FILE 7.077, Sn»pNc/L|8S HaS BEEN OPENED FOR WRITING ON UNIT i 


n 

Q 

a 

w 

« 

w 

i 

•tr 

O 



UNIT NO, 3 

component no; 1 SOLUTIONS* 

r g -sft o o ooF+o i e + w o do f + oi- 

UNIT % ENOFItEO AND RELEASED 
FILE 7076,SO«PNC/LIB$ H^S BEEN OPENED FOR WRITING ON UNIT 


UNIT NO, 3 

component no; g SOLUTIONS* 

, I2OOOO0I- +0-2 r H THH HW r' y frg rttUTODO O g + o * - 

" UNIT 3 ENDFILEO and released 

EIJJL. 70-7.5,, 80oPNC/LlB$ Has BEEN OPENED FOR WRITING ON UNIT .1 


UNIT NO, 3 

COMPONENT NO, 3 SOLUTIONS* 

itBooooofc+og »Eoooooor>-p? 

UNIT 3 EN D MLED an D RELEASED 
♦I AIL OONE, 

— rr — STtjP - - • 

STOP 777 


EIGURE 8.3-2 output for SG1N code for Test Case No. 3 (Concluded) 



FIGURE 


» n vc 

■ sampu rase a for soi.n code, 

S-Ot^Tte^ RV I TER A T TON PROCEDURE 
WITH t U DprpMPoSI T ION OF DIAGONAL 
BLOCKS, R F a 0 1 N G 1 K* OF ONE DErOMR. 

OSfD i.'JrtGflNAt BLnCK AND A pRE- 
DPTE RAINED GUESS AT ONE SOLUTION 
SUB* VP C TOR , AND SLIGHT OVER* 

RELAXATION, 

BATCH Mj.DE', 

three aircraft components, 

cnu 

i 

read DIAGONAL INFLUENCE SUB.MaTRT* 
ASSUeiAfEp HUH SECOND aIRCRaET COMPONENT, 

VNRF 

7063 

PERFORM DEroMPOSI HON Of THIS SUH. MATRIX 
AND STORE' 

INVE 

. tT J 

17 

DET UMRO-G I T T ON OF DIAGONAL SUB.M A T R T y 
FOR second AIRCRAFT rO M PnMfNT', 

USER rOULo STOP H£«fc aNo PERFORM 
SUBSEUONT COMMANDS LATER, 

M&w - Dt- TERM INF SOLUTION FOR Alt 

three aircraft components, 

r.NTi 

5 i 

. read ineluence sum. Matrix hce«j, 

AND T>£ P«E. COMPUTED A NO P PE* 

BfORfn DfrOMPOSlMOIv OF ONp Of 

the diagonal siiB* M atrices, 

vnrf 

705? 

70 60 
7061 
70tf? 

70Rh .1 

7064 

7065 

7066 

7067 

THE- CDM»- c«**ano is not u-sfo with 


8.4-1 Input for SCLN Cede for Test Case No. 4 


FIGUBBS-41 



THE ITERATION PROCEDURE t 

read constant sub-vector files, 

— * 

7068 

7069 

— • - - - 

SOLVE DIAGONAL BLOCKS By LU 
DECOMPOSITION, define the 
solut i on sur«vcctor files, -t he— 

PRE-OETERMINEO APPROXIMATE 
SOLUTION SUB-VECTOR files, 

AMOUNT OF OVE R .RELa X-A^TON, AND- 

AMOUNT of INTERMEDIATE PRINTOUT', 

SOLV 

„ 1 * 

7080 

7081 

ffrfr? - - 

O 

* 707 n 




*ttr- DtWF^ — 

STOP 


FIGURE 8.4-1 Input tor SOLE Code for Test Case No. 4 (Concluded) 


i FIGURES— 42 


POHaN £00* I Ion SOWING PROGRAM fSOLN)'. VERSION 1,1 

fSsz;**t*afiftS(Cscmit3SsX3stessarietszsB8iii3isxs 

DYNAMIC MfMORy » i 0000 

riMf. X 07/2,3/76 ) ! J 0 B X 3» 

♦ i SAMf'tF CASE 4 FOR SOfN coot". 

4t SOL u T t ON BY ITfRaTIOM PROCEDURE 

4j WIT^ III DFCOBPOSTT TOW OP OUGONaI 

♦ t HlOf-kg, Rt«nI N G IN OF one decamp. 

*t 06fn DTaOONaI, BLOCK and A PPE. 

M Of Tt KM I NED GUESS AT ONE SOLUTION 

m sub-vector, anq slight over. 

+ : RE L a * A 1 1 ON . 

+ S R/TC H Moot. 

+ : T.Pff A jPCP AF T COMPONENTS, 

4J fNU 

N{ Mp» , jTFRATs o 

4 ; BEAU PIaGOMaI. INFLUENCE SUB.Ma Iff t X 

+ } ASSOCIATED WITH SECOND AIRCRAFT COMPONENT, 

4 { VfjRF 

HIF 7<U>3, VN.pWf- Has BEEN OPENED FOP READING ON tfNJt I 

CRT A T I ON T I»*f 807/07/76 QHtQHt'iU 

unit i rewound and rEi Eased 

4* PERFORM OErOMPoSi T xON OF THjS SUtt.MATRfX 

4 : A v«D STORE , 

4 t I N V4 

METhOox < NTINVR3I 7096 Ir)fs 1 
ENTER SUB, PART tfFRAGS* 2 TNVs .2 

matrix loaded prom unit i3 using ? reads and n skips 
no partition TnG reoutrep 

EXIT PART 

*•* matrix inverse demanded 
Enter sub’ trip. n« 3 


FIGUHE 8.4-2 Output tor SOLK Code for Test Case No. 4 


FIGUHES-43 


EEPRODUCIbILITY OF THE 
OBIGKNAL PAGE IS POOR 



■■■ “ — — 

INVERSION OR DECOMPOSITION TIME <$EC.) ,007 

has been opened for writing on unit \& 

NTITL* j 7 itTO^m 

( TITL(I)» Is) #NTlTtV« DECOMPOSITION OF DIAGONAL SU6.M A T«IX 

UNIT-ta- ENqFlLED ANp REtg ' A ' S gp 

♦I FOR second AIRCRAFT COMPONENT' 

♦I USER COULD STOP HERE AND PERFORM 


♦ l NOW DETERMINE SOLUTION FOR all 

♦i three aircraft components, 

-n — im : 

NCMPs H iTERATt \ 

♦I READ INFLUENCE SUB. matrix files. 


♦I STORED DECOMPOSITION of one of 


♦ I THE DIAGONAL SUB-MATRICES, 



FILL 70SR,v N *PNC 
CREATION TIMEs07/G7/76 

UNIT 1 R EWOUND A N D RELEi 

FILE 7o*»O.VN*pNc 
CREATION T jMEa07/07/76 

llit T i . Hi unitiiK < ifA ■ tr . 


FILE 7o61.VN.PNC 
CREATION T JMF.507/07/76 
- U N I T — J- -REFUND- A-NO -Rttei 
FILE 7062 r VN.pNc 
CREATION T JmEs07/O7/76 

- U N I T 1 R E W OUND -aWO - R ELE) 

FILE 7o 9 6.IN-PNC 
CREATION TIME « 07/22/76 

FILE 706«;vn.pNc 
CREATION TIMEs07/07/76 

> i 11.1 * » . M r- i.i I « . ak. . k. h w . • r . 


FILE 7065.VN.PNC 
CREATION TimE.07/07/76 


has been opened for reading on unit i 
o«ioa*5ii 

tee 

has been opened for reading on unit \ 

O4:08|54 

• **• is 


has been opened for reading on unit j 
oaiosisa 


has been opened for reading on unit t 

06I08|5« 

j£ o 

has been OPENED FOR READING ON UNIT 4 

Hl 08 *ai 

m- — : 

has been opened for reading on unit j 

04106154 


has been opened for reading on unit j 

04106(54 


fIGURE 8.4-2 Output for SOLN Code for Test Case No. 4 (Cont’d) 


FIGORES-44 


UNIT 1 REWOUND AND RELEASED 

File 70D6, Vw-pNC has been OPENED FOR READING ON UNjr I 

CREATION T IHF.ff07/07/7#» 04t0flt3d 

UNIT t REWOUND and RELEASED 

FITf 70*7, Vh). PNC Has BEEN OPENED FfJR R£ AO 1 NO ON UNIT J 

■ CREATION T T MEs07/O7/7t> 0«S08|5« 

UNIT 1 ftpWnUttB AND RELEASED 

+j The COMB command is not USED with 

♦ I THf. ITERATION PROCEDURE, 

♦ I RcaD constant SUBSECTOR FTI.es. 

+ J nc pfe 

FILE 70CH,Rr.»pNC H A S BEEN OPENED for READING ON UNIT e> 

f Rp AT ION tJMf -07/07/76 O«T08J5« 

UNIT ? REWOUND aND RELEASED 

F II- 1 70*R,Br--PNC H a$ BEEN OPENED FOR READING ON UNIT ? 

CREATION f !«E*O7/07/7fc 0U*n«*S« 

.UNIT ? RE wnuM> AND RELEASED 

FILE 7D?0 # BC-PNC m A S BEEN OPENED FOR READING ON UNIT ? 

C RE A T'TTTN- T i f»F -"0 7 / 0 7 / 7 ft D«JDH|Sfl 

UNIT ? Rf WOUND AND RELEaSED 

♦ t SnLVE DIAGONAL Bl OCKS BY l.U 

♦ : m r.Hwposn ton, define thf 

+ j SHUT l ON SUB-VF.CTUp FILES, THf 

+r PkE. DETERMINED APPROXIMATE 

♦ r ■ ^LtrTrON SUB- VECTOR -FILES, 

+ | AMpi'NT OF OVFR.Rf-LAXAT ION, AND 

♦ I amount OF INTERMEDIATE PRINTOUT. 

4 1 sot v 

METHODS 
TPSs 1 

T PR T*T 5 

( OUT A, OUT AN 1 , 1st ,NCMP)s 
70 60 
70S J 
708? 

f INPS, tNPbNT , I* J ,NrMPl- 
0 t 

n ii 

7073 


FIGURE 8.4-2 Output for SCLN Code for Test Case No. 4 (Coji.t’d) 


FIGURES- 45 



HA* J I J 7 T 

OMEGAe | . 0^6 
ITERPb «$ 


ENTER SUg' TRfP, N« 5 
EXIT TRIP* 

FILf 70 9 fe t fNj.pN e H y ft " 8 ggN Q P gNED FOR REA D ING ON tmtT— a 

CREATION TIME.07/22/76 1*108141 

UNIT 4 REWOUND and released 


ENTER SUB' TRrP'. Nb ? 
EXIT TRIP 




H ER aT I- U N S OL t j TlO- N -PR OcEDu R g INTE R MEDIATE P R I N T Ou T 

initial guess 


FILE 7068.BC.PNC HaS BEEN OPENED FOR READING ON UNIT 4 

UNIT a rewound and released 

FILE ToPR^geaPNC — H A S BEEN O P ENED FOR READING ON UNIT g- ' 

UNIT 4| REWOUND AND RElEaSED 

FILE 7070.PC-PNC H A 3 BEEN OPENED FOR READING ON UNIT 4 

HttH 4 REWOUN D A N D RELEASED 1 - - ■ 


INFLUENCE OF UNIFORM A ND NON.UNJFORM PREESTRE A N . 

and r otation rate vector on a i rc r a f T~~ co yrP pNf ttt 1 — : 

,3115Jf; 08 15«017E+08 *78976Ef08 '.10474F + 09 T .i3062F*0R 


C^tPONEN^ NuMfUER 1 — " 

NUMBER OF ELEMENTSb S 
SOLUTIONS 

— ,i t87 ?t;oo — .mo R E+ o i — .gsg ^gg- ^ oi — 'ihrtorrtt 


Tl0^7tF+0t 


INFLUENCE OF AIRCRAFT COMPONENT | ON AIRCRAFT COMPONENT 


FIGURE 8.4-2 Output for S01N Code for Test Case N‘o. 4 (Cont'd) 


REPROP0CSBlUTy OF THE 
ORKKNAL PAGE IS POOR 


f IGORES-46 




.Mcofeje* .ssom + oft ,99ao9f + o? 

INFLUENCE Of UNIFORM AMO NOW. UNIFORM FREESTRfAM 
and ROTATION RaTF VECTOR ON AIRCRAFT COMPONENT p 
t 1^5BRf«0R .17B8OF+0V ,195597+09 

COMPONENT Mi iMrrER 7 

NUMBER (iF Ei E RENTS* 3 

SOI UT TONS 

,5188<»F } Q1 ' 6 7 366F f fl l 'J4672F + 02 

INFLUENCE Of aIRCRaFt COmPONFN t 1 UN AIRCRAFT COMPONENT S 

,3?«90F 4 C)f ,89 0B8F+0B 

T'NF L 1 *F NCf Of AIRCRAFT COhPONc NT ? ON a!RC«AFT COMPONFNT l 

,Bp9?Rf.0M ,S7<?t0F + O« 

INFLUENCE Of UNIFORM aND NON. UNIFORM FREE STRf’ A M 
AND PflTATfON- RATE VECTOR ON AIRCRAFT COMPONENT 3 
, l<m8F lfl T , I*773F + f)9 

Fill 70 7 5, Su-PNf HAS BEEN OPENED fQff READING ON-UNJT « 

CRFaTIOM 1 TMExn7/?P/7o litOhtPU 
UNIT /J REWOUND and released 

fOMPONg-MT Nt/Mpgw ■% 

NUMBER of El FMENTSe ? 

SOLUTIONS 

,900007,01 , tOOOOF+O^ 

ITERATION , LARGEST pEL&TlvE Errors .626337E+00 

ITERaTTON s LARGEST RELATIVE ERRORS ,97«581F*0t 

ITFRaHON LARGEST RELATIVE ERROR* ,U?70#F«O! 

iteration a largest relative epndr s .-?^7379e-02 

ITFRaTIUN SOLUTION PROCEDURE INTERMEDIATE PRINTOUT 


FIGURE 8.4-2 Output for SCLN Code for Test Case No. 4 (Coat'd) 
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tTE RATION m« Hi g tLOWABl E — 

INFLUENCE of AIRCRAFT COMPONENT 3 ON AIRCRAFT COMPONENT i 

---.-frW Sr;Ofe — ,? 5 5 TtEfO 7 -,- T 5 24 e F + 0 7 ,20585 r* 0 8 — ,tmfeOt+tl« ~~ ’ ~ 

INFLUENCE OF AIRCRAFT COMPONENT g ON AIRCRAFT COMPONENT j 

— ;\ i«-it + 0« — f .50R«*«E+05 — ^r^SSr+TT*" ' Sfr^OEiTT* — ’,33008't + OT — 

P.. 0 *fP* g NT yu^o£R i 

SOLUTIONS 

,10000Ri01 ,2OGO0FfOl ,30000EfOI e «OOO0F*0l *S0005F+OI 


INFLUENCE OF AIRCRAFT COMPONENT j ON AIRCRAFT COMPONENT 
t 35«l<»E;0fc ;i«7P3Ea08 'S5o72!U07 


INFLUENCE Of AIRCRAFT COMPONENT g ON AIRCRAFT COMPONENT ? 
,I0212 f + 0« ,27111E + o8 ,6«020F + OS 


COMPONFNT NUMBER 2 

NUMBER OF ElEMENTSs 3 

— SOttittO N S 

,SRRR&e, 01 .70000E+OI •79999Eao1 

INFLUENCE OF -y l RC RAF - T - COM P ONE N T f 'ON- sT R CW T COMPONEN T y 

■ |7660E *07 ;«6302F + 06 

— INFLuENCE OF a IrCraFt COmPONEN t 2 — ON .aIrChaFT COMPONEN T 3 

,5?a8«E + 0« ,2«221 Et08 

— COMP O N E N T nu m b e r 5 : 

NUMBER OF ELEMENTS# 2 

SOLUTIONS 

.ROOlHfi 01 71000&E+02 : : - 

Iteration 5 largest relative error# »«7or22e«.03 


ITERATION t, LARGEST RELATIVE ERRnR* 9 82U8PE*0« 

FILE 7O8O,SO«PNC/L!0» H*S BEEN OPINED FOR WRITING ON UNIT | 


FIGURE 8.4-2 Output for SOLN Code for Test Case No. 4 (Coat'd) 


RilPR^UCliilLlTY OP THb 
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CREATION T?Pf? s 07/22/76 1110*105 

UNIT ^ ENDFIlFO AND RELEASED 

UNI T NO, 3 

COMPONENT HO' 1 SOLUTIONS! 

,9999999F*00 ,i999999F + 0l ,?999998E + 0t .3999996E+01 ,50000BIF+0| 

Fit? ?OeH , S^’-pNC/t ” HaS BEEN OPENED FOR WRITING ON UNIT 3 

CREATION T T * E * fil/W/lb llto^lnfc 
UNTT ? ENDFuEO AND REt EASED 

UNIT NO 3 

COMPONENT VO' 2 SOLUTIONS* 

,^99970FF0t { 7ft00O021N 0 1 ' 79999ft J f + 0 I 

FIL,E 70*«? t S»«PNC/LlH* ‘ Has BEEN OPENED F OR WRITING ON UNIT J 
CREATION TTMf, at 07/22/7 b 11*09(06 
UNIT 3 FNDMlFU AND RELEASED 

UN'TT NO, 3 

COMPONENT ho' 3 SOiUTtONSj 

,9('0D007F + 01 .IOOOOOOF+0? 

NO. HE JTfR^TiONS TO CONVERGES ft LARGEST RFLATlyE ERROR* '.«?! 1 86E«0<> 


SOLUTION VFrTfiP NO' 2 


ITERATION S 0 L 1 1 T 1 0 N PROcEDiiRE INTERMEDIATE PRTNTOnT 


TNJTI 

At <S»»ERS 









FILE 

7<K*R,Hf -PNC 

was 

BEEN 

OPENED 

FOR 

reading 

ON 

UNIT 

n 

UNIT 

a rewound and 

RELEASED 








FILE 

7oo q »Bf:-PRC 

Has 

BEEN 

OPENED 

FOR 

reading 

ON 

UNIT 

4 

UNIT 

a rewound and 

«E l EaSED 








FJU 

7O70.RC-RNC - ‘ 


BEEN 

opened 

FOR- 

reading 

ON 

UNIT 

a 

UNIT 

41 pE WOUND AND 

pElfASED 









FIGURE 8.4-2 Output for SC-LN Code for Test Case No. 4 (Cont'd) 
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— influen c e - op uni f o r m and NON^u y r r r tm ^ tr ee- s tre a m 

AND ROTATION RATE VECTOR ON AIRCRAFT COMPONENT 1 

.62 3o 7 E;0» .JOeOJE + O 9 ,l5T95E;<jR ,2(>947E + o9 .2M23E*o9 

component NUMBER 1 

NUMBER OF ft EMENTSa S 

»O tUTrON 9 

,lRR 7 ae;oi ; 3 R 0 iRE + 0 i *. 565 a 0 E + 01 r . 7 ia 60 F + 01 .', 2 U« 2 F 402 . 

I Nr Li lc.NCt Dr A, I “C* Sr T C“““ UNEN^ p uN TfTreC“A F j CQMruNtNj 2 

.12201F+09 ,52t4«F + 0« ;iRR22F+08 


I N FtUE N CE - OfM ? N l POR M aNO N ON. UNIFO R M F REEST R EaM 

AND ROTATION RATE VECTOR ON AIRCRAFT COMPONENT 2 
,3H76p;0R .35775E + 0R .39U8E + 09 


COMPONENT NUMrER ? 

number of elements* 3 

SOLUT I O N S 

>3767Ef01 r .t3«73F*02 ;29544F*02 


# 657SlF + 07 # t78l6E + 07 


f -ON-AtPCRA^T ODMPONENy- — * 


I N FLU EN C - E-OF - xI R C R aFT -C O MP O N E N T ? ON AIRC R a PI 1 COMPO NENT T 

.ITiBPF + O 1 * ,7<Hl2OF + 08 


“A-NO-NON •UNtFORM’-pPEf SETTLE A"M — 

and rotation r a te vector on aircraft component 3 
# 3R«96e;09 ;335<|6E + 0R 

FlL r TOF^P wP N C - — Ha3 Be e n OP e N gn FO R READING B N -tfW" "T 

CREATION TlME«07/22/76 Ut06l2« 

UNIT U rewound and released 


COMPONENT NUMrER 3 
NUMBER OF Ei EmENTSb 2 

S OLUT I O NS 

, 1 &OOOE *02 * 20000E+02 


FIGURE 8.4-2 Output for SOLN Code for Test Case No. 4 (Coat'd) 
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ITERATION 

1 

LARGEST 

relative 

ERRORa 

.626S37F+00 

ITERATION 

? 

LARGEST 

RELATIVE 

ERRORS 

.974581E-01 

iteration 

% 

largest 

RELATIVE 

ERRORS 

, l5?70UE,0t 

ITERATION 

u 

t AftHgST 

RELATIVE 

ERRORS 

,?R757R.r,oa 


IHRaTJON SOLUTION PROCEDURE INTERMEDIATE PRINTOUT 


ITERATION WI T H Jo ALLOWABLE 




INFLUENCE f)t- AIRCRAFT COMPONENT ;» 

,1?01/E 4 07 ,47Jq5F407 .150S0E+08 

ON AIRCRAFT 
* OB 

COMPONENT 

.99921F+0B 

t 

INFLUENCE 01 aTRCRaET COMPONENT ^ 

# 2^?alE*0R # 4l9?8F*0b t 33flfe7E f 06 

ON aircraft 
, 1 7340 F +07 

COMPONENT 
• 660 1 6F + 0 7 

1 

COMPONENT NMNrrpT? 1 

NUHRFP OF ft E M F N T S a S 

solutions 

,20000f*01 ’40000F + 01 'boooof+oi 

f . B0000F + 01 

I0001F + 02 


INFLUENCE OF AIRCRAFT C0HP0NFN t | 

,?^N8 7F+t)8 f l t0t4r + OB 

ON aTRCRaFt 

COMPONENT 

fi 

INFLUENCE OF aIRCRaFT COMPONENT ^ 

,20423?, 0* 08 r J2804F + 0R 

ON AIRCRAFT 

COMPONENT 

? 

COMPONENT Ni iM i* E 4 ? 

NUMRfR UF ft fWENTS* 3 

SOt. UT TONS 

, 1 <?000f ,02 .14000F+02 .16000F+02 


* 


INFLUENCE OF aTRCRaFT COMPONENT { 

.3S3iRf t 07 ;<?2fcft3F408 

ON AIRCRAFT 

COMPONENT 

3 

INFLUENCE of AIRCRAFT COMPONENT P 

,llS7?F40O ,48443F+0« 

ON AIRCRAFT 

COMPONENT 

t 


FIGUBE 8.4-2 Output for SC1N Code for Test Case No. 4 (.cont'd) 


{jjr Tuj 

OJUaNSL PAGE IB POOR 
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COMPONEN T NUMpgR 3 

NUMBER OF ELEMENTSb a 
SOLUTIONS 

rt fl ooOf»oa *, 2000 0F » o^ 


ITERATION * LARGEST RELATIVE Errors . ,470922E«.03 


ITFRATION h LARGEST RELATIVE ERROR* ,82U86F.04 

FILE 7080,SO-PNC/LJB$ H*S BEEN OPENED FOR WRITING ON UNIT 5 


UNIT NO, 3 

C O M P ONENT NO r . 1 Sfo UT IO NS I 

.20000O0E+01 .S999999F+QJ ,5999997E;0i ,7999992F*0I ,1000018F+02 

FILE 7081 ,80-PNC/LtB* H*$ BEEN OPENED FOR WRITING ON UNIT 3 

Wf — r- e - NPHteo-ANP relea sed 


UNIT NO, 3 

C OM PO NENT noI — -t — 8 P riiT t ON S-|— 

,U9999«E*02 B 1 400000F+02 ,tS99996E;02 

FILE 7082.SO.PNC/LIB* H*9 BEEN OPENED FOR WRITING ON UNIT 3 
turn 3 E N PFILED AND R E LEASED - r - 


UNIT NO, 3 

■ CO MP O N E N T NO* ^ SOLUTIO N S; 

,J M 00001E402 r ,2flOOOOOF*02 



♦ 1 A| l Pone, 
♦t stop 

sto p 777 


-~^t1S6E.0« 


PIGUBE 8.4-2 Output for SCiN Code for Test Case No. 4 (Concluded) 
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TRUf- 

s a m p l t case b for soln code 

HV A two STe? PROrloURff. 

E1RS7, 1W(3 AIRCRAFT tOMpo^EfvjTQ ARE SOLVED 
ny OlRFfT tu decomposition, THEN THE 
INFLUENCE Of A THIRD aIRC«aFT COMPONENT 
ON f Hf. FIRST Two IS COMPUTED THE 

ttfRatton procedure, the direct solutions 
for tmf first two aircraft components a«e- Used 
as Ini rut guesses to the three component 

ITERATION SOLUTION' 

CNTL 

? 

RE A 0 t NF L Uf Nf F SUB. MATRIX FILES 
" pur thc first two aircraft 
compone nts', 

VNRE 

705R 

7060 

7062 

7063 

SINCE ENTIRE Matrix PITs IN corf, 

COmH Command Is NOT needed', 

REaO constant suh. vector files 

bcre 

7068 

• 7ft** - - 

SOLVE pv [U DECOMPOSITION, 

and define the solution 
SU f(«VT f. T OR OUTPUT FILES 
SOL V ■ 

0 p 

?n#3 

7084 

. 0 

U$EP rOttLO STOP H£ p € ANp P*?rFOrM 

pfMAtMNf; commands Later, or he 
may continue, 

~i nf-Tt-RMiNf !NFLt;ENre nF THIRD 

AIRCRAFT COMPONENT VIA ITERATION 
PROCEDURE’ 

CNTL 

* 1 

REREAD aLL INFLUENCE SUB_M A Tf?tx 
f It FV 


FIGUfi E 8.5-1 Input for SOLN Cede for Test Case No. 5 


*> 
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ftEPEODUCfelLITy OF THE 
0 EI 9 CSTAL PAGE IS POOR 




7067 


BCRE 


THE COMH COMMAND MaY NOT BE USED 
REREAD aU CONSTANT SUB-VECTOR FILES, 



Jviiuirl 


DECOMPOSITION, DEFINE THE 
SOLUTION SUBSECTOR fILES, THE 


> 1 ^ ■ “i-* 


il ul l^l<4 * M 


SOLUTION SUB-VECTOR FILES, 
MAXIMUM number of ITERATIONS, 
- a m ou n t -O F OVE R »RELA *- mON - , 

AND AMOUNT OF INTERMEDIATE 
PRINTOUT'. 


0 

TOSS 
— mt- 
7087 

0 


2 708a 

10 



FIGURE 8.5-1 Input for SCLN Code for Test Case No. 5 (Concluded) 
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ftEPRODU C5B1LITY OF THE 
ORIGINAL PAGE IS POOR 



POTfAN EQUATION SOLVING PROGRAM fSOtN), VERSION 1,1 

titisciliitttttsieieiBttisxicstctcistasssaiiBKesssss 

DYNAMIC Mt MORy x 10000 


TI«e * 


1 1 


Sample CASE «S for SOLN CODE'. 

SOLUTION HY A 7*0 STEP PROCEDURE'. 

FIRST TWO AIRCRAFT COMPONENTS a r E SOLVED 
Mr I'IWEcT LU DECOMPOSITION, THEN, THE 
INrtUpNCE Of A Third AIRCRAFT COMPONENT 
Of; THE FIRST Two IS COMPUTED hY the 
ITERaTTON PROCEDURE, THE DIRECT SOLUTIONS 

for Tkf first two aircraft components are used 

AS INITIAL GUESSES TO THE THREE COMPONENT 
ITERATION SOLUTION^ 

? TjERAis 0 

REAP INFLUENCE SUB. MATRIX FILES 
FUR The first TWO AIRCRAFT 
COMPONENTS^ 


+ * 

♦ X 
+ 1 

♦ X 

♦7 ’ • 

♦ * 

♦ X 

♦ 1 
+ * 

♦ X 

♦ i rNTi 
NCMPs 

♦ I 

♦ ! 

+ 1 

+j vNRF 

ph t . 7 -ov? , v*-pn£ has been opened for reading on unit j 

TREAT ION T IHE»07/07/7t> 0<TTOM<5<4 

UNIT 1 REWOUND AND RELEASED 

FILE 70BO.VM.PNf HAS BEEN OPENED FOR READING ON UNIT t 

CREATION T j ME »07/07/7P 0«*OB|5« 

UNIT 1 REWOUND AND RELEASED 

Fltf-ttn^^. vm*pnc . mas been opened for Reading -on unit \ 

CREATION T TME*07/O7/7ft Q4:0*|5« 

UNIT l REWOUND ANO RE|.EaSED 

F ILF 70B3,Vn.PNc H A S BEEN OPENED FOR READING ON UNIT \ 

CREATION r TMFb07/07/76 0<1I0BI5<< 

UNIT 1 REWOUND AND RELEASED 

*7 S-i Nf-fr -f-NT IRE Mjffr IX FITS IN CORfc , 

+j comb COMMANo IS NOT NEEDED. 

+1 RI ad constant sub.vector files. 

♦t bcre 

FILE 706R,Rr*PNc HaS RggN QPeNEo FOR REAolNfl ON UNIT 2 


FIGURE 8.5-2 Output for SCLN Code for Test Case No. 5 
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CREATION T!M£*07/07/7<> 04»08l5<l 

W 2 R EW O UN D -* N O - R ELEa8ED 

FILE 70fe9,8C«PNC H*S BEEN OPENED FOR REaOING- ON UNIT £ 

CREATION TiME»07/07/7b Oa»08|5« 

turn- g r e w oun d and r eleas ed- r— 

♦* SOLVE BY LU DECOMPOSITION, 

♦I AND define the solution 

ti SUgwV E CTO R O UT P U T FILES. 

♦ I SOlV 
METHODS A 

IPS* 0 r-— - 

IPRINT# a 

(OUTA,OUTANT # |ei,NCMPj* 


7 O B 3 flf 

708 a n 

ENTER SUB! part «FrAG 8 « 10 INV* »2 



no partitioning required 
exit part 


“TNO w -SKIPS 

\ 


ENTER SUr^ TRIP. Hm (i 
EXIT TRIP 


OUTPUT 

SOLUTION TIME (SEC r J ’.Oil 


SOLUTION vErToR NO r . 1 


FILE 7085,SO*PN C /LJ8S HaS BEEN OPENED FOR WRITING ON UNIT T 

cr eati on t ime -«• - vnttrr * 1 i m i s s 

UNIT NO, 3 

C OM P ONE N T NO" i SOLUTI ONS ! 

,1000S5Sr + 0| ,200S3a2F;oi t 3003<>88E;01 t 399038iE+0l t a9«J5J2E*0i 

UNIT 3 ENDFILEO AND RELEASED 

- F - It -r^oS a.SO^PNC/LlBl M as B EEN ope ned for w riting on unit — 5 

CREATION time , 07/22/76 UtOR|56 

UN I T NO, 5 

COMPONENT NO'. 2 SOLUTIONS* 


FIGURE 8.5-2 Output for SOLS Code for lest Case No. 5 (Cont^) 
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.SMSfeflSE + Oj ,6fe386^ft(- +0 j , ! «0 j Ef + 02 

omit i ewitfo-Awo released 
output 

SOLUTION TI*E rstC’.) .001 


SOtOTlOW Vf r Tfjff tiO\ ? 


f!U 70«i.Sn.^w C /t- IH ^ Hi A S BEEN OPENED F OR WRITING ON UNIT 3 

UNTT NO, 3 

COMPONENT fif'. 1 SOU) r TOWS f - • 

»?O01l()7F+0| ,«006fe«lF*0I .6007377L +0 1 t 7<?80 762h 4 0 1 , 900662 3f ♦ 0 1 

.UNIT 3 E^DFIiED AND RELEAS'D 

f ILf- 70#«' Sf>«PNC/tlB$ HaS BEEN OPENED FOP WRITING ONUNJT 3 

uNTr NO, 1 

COf+PONE-M Nfi', 2 SnLtif?ON$t 

* l 1 67 1 ? 1 F a 02 .n277«0F + 0a .2R22801E+Q2 

UNIT 3 EMSMlFD AND RELEASED 
T« U S E n c^Ot-D STOP Hfcft£ AND PERFORM 

+ « RfM^INING r 0MHaN 0 S LaTE«, OR he 

M M ' y CONTINUE , 

rt ■ - - Tri; rr*MlWf -TNft.tjENcf- Of fKt*D 
•f ! AIRCRAFT COMPONENT VIA ITERATION 

+ 1 PwfVEDlJRE, 

+ | f N T t 

NCMPs •* TTFRATfi j 
♦ I REREAD ALL INFLUENCE SIJB.M A TR I x 

>1 ftLFS. 

+ 1 VNRf 

F i Lp 70SR.vE.RNc has been opened for Reading on unit i 

CRFATjnN T I ^Es07/0 7/7o 0*UQ«|S<J 

UNIT f RfcWfiUND and RELEASED 

Fit E 70h0.VN.PNc H A S BEEN OPENED FOR READING ON UNIT , 

•eRFATtON-Tji-Fx^TAfrt/ft 0*70*7^* 

UNIT t RFwiiUnD AND released 

FILE 7oM .VN.pNc H A S BEEN OPENED FOR READING ON UNIT * j 

CREATOR T t hE*07/O7/7W 0«T0S|S« 

unit \ rewound and released 


FIGURE 0.5-2 output for SCLN Code for Test Case Ko. s (Cont’d) 


a^RUUUCIBlUTY OF'THif. 
PIGU RES-5 7 , jpjomAL PAGE IS POOP 




itiBkiiMiaar ivifh fli .1 H HP! 


♦I DECOMPOSITION, define the 

it SOLUTION SU8«VECT0R PILES, 7 Hf. 

Vt SOLUTION SUS. VECTOR FILES, 

♦i maximum number of iterations, 

V| AHOUN ir OF OVgRwRgLAXATlON ; 

ij AND amount of intermediate 

it PRINTOUT., 


METHOD* fi 


FIGURE 8*5-2 Output for SOIN Code for Test Case No. 5 (Cont'd) 
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IPS* ? 

JPHfMhr- - • • 

( OUT A# OUT ANT , I*| , NCHPj- 
708*4 n 

708 b a 

7087 

flNPS#rNPSNT,!#l#NCMp)* 
-7083 - a- 

7084 r, 

0 't 


MAXJT1* t 

OMFGA* 1 t 1 (* o 
TTERP* 0 

/* 



FNTEP SUB", 
EXIT TRIP 

TRIP. 

N* 

h 

fNT£R SUp' 

fxit trip 

tptp. 

N* 

S 

fNT£R SUh' 
MIT TRIP* 

TRIP. 

W R 

2 


S-Ot-t+TtfTN-^rtrtCrP NO 1 ,; t " 


p 1U 

70b8 # BC-PNC 

HAS 

BEEN 

OPENED 

F OR 

reading 

ON 

UNIT 

n 

UNIT 

(1 Rf' WOUND AND 

RELEASED 








FI IF 

70O9,BC«*PNC 

has 

BEEN 

OPENED 

FOR 

READING 

ON 

UNIT 

a 

UNIT 

a rewound and 

RELEASED 








PILr nr, mpar. ■ 

has- been 

OftNED 

FOR 

-PE AD I NO 

ON 

UNIT 

a 

UNIT 

a RfWiHTNO aNd 

released 








FIU 

7085.sd.pnc 

Has 

BEEN 

OPENED 

•FOR 

reading 

ON 

UNIT 

a 

r Rfc A T ] QN n*F*07/g?/76 IHORIS^ 







UNIT 

« REWOUND AND 

released 








FILE 

70B«,SU.PNC 

h*s 

BEEN 

OPENED 

FOR 

reading 

ON 

UNI T 

4 


CTTfrAT jfTN--TyiKt-*lTT7P?'/7<i tlTO^TSft '■ 

UNIT a RfcwnUNP and RElEaSED 

I TE« A T t 0N 5 t'APGEST R£tATl v E ERROR* , 476246f>00 

flGURE 8.5-2 Output for SOLN Code for Test Case No. 5 (Cont'd) 
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ITERATION g LARGEST RELATIVE ERRORb e U5208E+00 


ITERATION 


5 LARGEST RELATIVE ERROR* 


,25S57«Eo0.t 
V303b3Tr«0^ — 


ITERATION 5 LARGEST RELATIVE ERRORb .122R35E»03 

ITERATION h LARGEST RELATIVE ERROR* ,S7S02OE*O« 

FILE 70«S,SO*RNC/UBs HaS BEEN OPINED FOR WRITING ON UNIT 3 

C R EATION Tt*g b 0 7/22/ 7 6 n > T fi ?«b 

UNIT 3 ENDURED and RELEASED 


uru i Ptu, 3 

COMPONENT NO' l SOLUTIONS! 

, 10000001 + 01 ;2oOO0OOF+Ot JOOOOOOE+Oi 

FItE TO O fe t SO*PNC/LlSS H &S BEEN OPENED F O R WR T 

CREATION TIME s 07/22/76 llUot«7 
UNIT 1 endfiled AND RELEASED 


;3R99<mg + oi 

JG“ ON _ tTNli' — I 


UNIT NO, 3 

COMPONENT NO r . 2 SOfUT IONS | 

— — — i-ym997E;ot — , 7 o o o o o o r »ot — 7 8 o d o o iirm - — 

FILE 70®7,SO«PNC/L!Bs HaS BEEN OPENED FOR WRITING ON UNIT 3 
CREATION TIME s 07/22/76 11*10848 


UNIT NO, 3 


t 899999*E+0i , 1 000000E *02 


# 5000026? +0 1 


-^Twaor.tra- 


SOLUTION vErTOR NO' 2 


UNIT a REWOUND AND RELEASED 

FILE 7069.BC-PNC HaS BEEN OPENED FOR READING ON UNIT a 


FILE 7070, BC. PNC 


has been opened for reading on unit a 


PIGURE 8.5-2 Output for SOLS Code for Test Case Ho. 5 (Cont^) 
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UNIT a RtwuUND AN(? RELEASED 

F-itf 70H3,3+>»P*C' Hj\S BEE* OPENED FOR READING ON UNIT u 

CREATION T I M E*0 7/22/76 1110*155 

unit a rewound and released 

FILE 7oS« t §o«fPNc Has BEEN OPENED FOR READING ON UNIT a 

CREATION T T ME*f)7/22/7f> 1110*156 

UNIT a Rt'WOUND AND RELEASED 


ITERATION 

I 

largest 

relative 

error* 

,«762«6E + 00 

ITERATION 

? 

t ARGfST 

relative 

ERROR* 

, US208F + O0 

ITERATION 

% 

LARGEST 

RELATIVE 

ERROR* 

,23557UE,01 

ITERATION 

0 

largest 

relative 

ERRORs 

,303637g*0? 

ITERATION 

<*> 

LARGEST 

RELATIVE 

ERROR* 

, 122435F*03 

iteration 


largest 

relative 

ERRORb 

,573020?. 0« 


F-iif ■fOHs t $v»pNC“/ , t.-if?i5' ha-s been opened for writing on unit -5 
UNIT 3 ENDHLED and released 

UNIT NO, % 

COMPONENT NO' I SOLUTIONS! 

,?000000t + 01 ."OOOOOOFaOI ( 5R99999E + 0 1 ,7999999E + 01 . I000005f +02 

FILE- 7DS#T-,-ST}wRNe/i-JR7: «*$ BEEN OPENED FOR WRITING ON UNIT ? 

UNIT ^ endfiled and released 

UNTf NO,, 3 

COMPONENT no' 2 SOI UTlONSj 

,Jl999R9F+02 , 1 900000F+02 ,16000021*02 

F^tf-tOM | -Stj»RNcytI»r' “--HsS BEEN OPfNED FOR vrtmiNG ON UNIT 1 
UNIT 5 endfilfd and released 

UNTT NO, 3 

COMPONENT 3 SOi UT JONS j 

, 1 /999R9F + 02 < .2000000T +02 

NO, OF ITERATIONS TO CONVERGE* h LARGEST RELATIVE ERROR* r .973C20E-04 

♦ i ai i done, 

♦ t S T ftP 
STOP 777 


ilGUSE 8.5-2 output for SCLH Code for Test Case No. o (Concluded) 


FIGURES-61 


ebproducibeuity of the 

OSKEMAL PAGE IS POOR 



true 


SaMPU CASE 6 FOR SO(,N CODE, 


that a rectangular matrix c 1 7 x i o t 

IS SOLVED INSTEAD OF A SQUARE 


PROCEDURE WITH INVERSE OF DIAGONAL 
blocks, reading in of one inverse of a 


IfJfV I'i 1 !! » ■«— -I— ■ I W V It f 1 


CNTL 


GUESS at 0Nr SOLUTION SU9«VECT0R, and 
SLIGHT DVERwRELAXATION, 


three AIRCRAFT COMPONENTS, 


read DIAGONAL INFLUENCE SUB-MaTpIX 

associated with second aircraft component, 

, — „ ... — - — 

7081 

perform inverse of this sub*naTrix 


INVF 


1 7097 

INVERSE OF DIAGONAL SUB. MATRIX 
FOR SECOND AIRCRAFT COMPONENT’, 




■ ■ 9 1 .1 4 ^ 4 * 


SUBSEQUENT COMMANDS later 
now DETERMINE solution foA 



.1 «f» rTMf.» Bk 


and the previously computed and 
STORE ft INVERSE OF ONE OF THE DIAGONAL 



FIGDBE 8.6-1 Input for SOLN Code for Jest Case No. 6 


FIGOBES-S 2 



FIGURE Q 


708S 

the ct'Mh comma wo is mot used with 
■ THE IfJrHATIttM- PPfleffDupt ’ 

RCR£ M£aD Cn^STAM SUB*v£CTnR FILES. 

708fe 

7087 

7088 

‘ %?***[* lnV ?*S* Of *>I Attest StOCKS 
HOIjEEHOlOER PROrEDuHE, DEFINE THE 
SOLUTION SUh»V.ECTOR FIlES. THE 
m-&fcTfRHlPJED APPROXIMATE SOLUfTON 
SUS-yfCtOp FILES, AMOUNT of 
OV fK.kE, AXATTGN, And Amount of 
PRINTOUT 

SOL V * 

1 1 
'7088 ■ - 

7089 

7090 


8Y 


1 707 & 
10 

1,85 

S 

ALL DONe 
—9ft)P 


6-1 Input for SCLN Code for Test Case No. 6 (Concluded) 


IIGUEES-63 


■iCESPBODUCIBlLiry OF THE 
OfilSfNAL PAGE IS POOS 



POTfAN EQUATION SOLVING program (SOLN), VERSION 1,1 

tiftmiMitfstssatiKsittiiiisissaaKEstssiBSKaaKig 


DYNAMIC memory c 10000 


SAMPLE CaSe, 6 FOR SOLN CODE'. 


V) ' ■ * ** * 1 C"’ ’ ’"'Of' ■» I v. AUfc' I 

that a RECTANGULAR MATRIX UTXlo) 
IS SOLVEO INSTEAD of a square 




PROrEDURE WITH INVERSE OF DIAGONAL 
BLOCKS, READING IN OF ONE INVERSE OF A 


I'L’IIM-IBHI 


♦ I 

♦ I CNTL 


GUESS AT ONE SOLUTION 8UB»VECTOR, AND 

slight over-relaxation; 


wn rivut t 

Three aircraft COMPONENTS, 


read DIAGONAL INFLUENCE SUB.MaTRJX 
ASSOCIATED WITH SECOND A|R C PafT COMPONENT’, 


file 70*1' VN»pNc 

CREATION T!M£ B 07/Q7/76 


HAS BEEN' OPENED FOR READING ON UNIT j 
0«I09I01 


♦| PERFORM INVERSE OF THIS SUB-M^TRIX 

♦I AND STORE, 


METHOD* , NTINVNo 70*7 IDC# 0 
ENTER SUB. PART fiFRAGS* ? TNV* 0 < 

no partitioning required 
exit part 


and (r — Skips 


*** ' matrix inverse demanded 


EXIT HRfcc: INVERSE ON UNIT ft 


FXGOBE 8.6-2 Output for SOLN Code for Test Case No. 6 


FIGURES- 64 




INVERSION OR DECOMPOSITION TIME.fSfcc,) .145 

FILE 70^7, «AS SEEN OPENED FOR NR I TINS ON UN|7 ft} 

CREATION TJmf « 07 / 22/76 l!ti 5 H 2 
NT IT! « \f 

(TITlf I)# t*I*NTjTU>« INVERSE OF DIAGONAL RUB.M*TRIX 
UNIT IS ENOf'UEO AND RELEASED 
+ | F'iR second AIRCRAFT COMPONENT. 

♦ I W»f*- r wtt> STOP HERE iND PERFORM 

+ 1 S!'RS fc 'olj E NT f OMm aNqS LaTEw’. 

+ | NOW DETERMINE SOLUTION for 

+ 1 ALL THREE AIRCRAFT COMPONENTS*. 

♦ I CNU 

NCMP* \ HERAT, 1 

♦ * READ ?*rLt)ENCt &U8.MATRTX FILES 

+i and ThF previously computed and 

M STOKED INVERSE OF One OF the diagonal 

+ 1 S 1! 8*Ma TRICES, 

+| VnRE 

FILE 7 Q 7 7 , V^-PNC Mas reen opened FOR READING ON UNIT j 

CREATION f Hf *07/07/76 03*0^10! 

UNIT 1 Ft WOUND AND RELEASED 

FILE 707ft, VN.RNc HAS BEEN OPENED FOR READING ON UNIT t 

CREATION I tMf *07/07/76 ORlORinl 

UNIT i REWOUND AND REl E aSED 

FILE 70 79, VN.pNC W AS BEEN OPENED FOR READING ON UNIT j 

CRtAfiftR fTHF*O 7 / 07 / 7 ft OUjOPrOI 

UNIT t REWOUND AND RELEASED 

FILF 7O60,VN-PNC h aS BEFN OPENED FO« READING ON UNIT [ 

CREATjON TiME*07/07/76 04|0<>ial 

UNIT i REWOUND A NO RELEASED 

FILF 7fl97»I**PNC HAS BEEN opfn&d FOR READING ON UNIT 0 

CREATjON fywfr- *- 67 /? 2/76 1 UHH 2 ’ 

U*N I T u REWOUND ANn RELEASED 

FILE 7062.VN.PNC H*3 BEEN OPENED FOR READING ON UNIT \ 

CREATION TIMt*07/07/76 oa * O^t o 1 

UNIT 1 REWOUND ANO RELEASED 

FILE 7c>ft3,VN-PNC H AS BEEN OPENED FOR READING ON UNIT j 

■watjon fiMg » c 7 / 07/76 <mo*tot 

UNIT ! REWOUND AND RELEASED 

FILE 70»«,VN.PNC HAS SEEN OPENED FOR READING ON UNIT \ 

CREATION I J*E*07/07/7e 0«IO9f0l 

UNIT i REWOUND aND RELEASED 


FIGURE 8.6-2 Output for SCLN Code for lest Case No. o (Cont'd) 
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FILE 7085.VN.PNC 
C R EATIO N TIH£»07/07/76 


HaS been opened FOR READING ON UNIT 

twr^rrl 


UNIT 1 REWOUND AND REL£a8ED 
♦I THE COMB COMMAND IS NOT USED WITH 

READ constant sub.vector files'. 


+ 1 

I I pf bit 

FILE 70 8 fa.SfP N c 

CREATION UME»07/07/76 0«*09JOl 

UNIT a REWOUND and released 

"m-E-7 0 8 7.BC. pNc H^S-&e gN - 0 P E N&0 

CREATION TImE« 07/07/7S 0<ll09|<H 

UNIT E REWOUND aNO RELEaSEO 

■ HIE 7 0 8 8.BC» PN C NgS BEEN OPENED 


h as been o p ened f o r-rf apI n s on unit — ? 


PGR - RTAOtNO-GN -tTN I T g 


FO R R EADI NG- O N - U N IT — ? 


CREATION TimE*07VO7/78 04 S 09 1 0 1 

UNIT 2 REWOUND AND RELEaSEO 

+1 C OMPUTE INVE RSE- -Of-Pi-AGGN A L" BL OC KS~gt 

♦I HOUSEHOLDER PROCEDURE, DEFINE The 

+1 SOLUTION SUB. VECTOR FILES, THE 

PR g.P g TgR H iN g ft- A P P R OXIMAT E SOLUT I ON — 


■tr 
+ i 
♦ i 

♦I SOL V 
METHOD* 
-IPS * 


SUB.VECTOR FILES, AMOUNT 0? 
OVER-RELAXATION, and AMOUNT of 
I - NTE RM e pT-m- PR-I NT-OUT 1 , 


— 1 

IPRINT. n 
(OUTA.OUTANT, j«j ,NCMP)« 
7068 


7089 a 

7090 ft 

(INPS y |N F s NT » i m»NCHP): 

0 ft 

0 ft 

7073 a- 

MAXITI* t n 
OMEGA* I . OSO 

• ITERP* % 


*** matrix inverse demanoed 


ENTER 8U8' HREC; 


7 ROWS 


s columns; matrix on unjt ia 


FIGURE 8.6-2 Output for SOLN Code for lest Case No. 6 (Coat'd) 
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EXIT HREc; iNVpRSf, ON UNIT 12 


SUB. ATN tr a nsposed A 

FILE 7097,IN-pNC 

CREATION tI*fci07/2?/76 


7-bv» s matrix 
has been opened fur reading on unit 
llUStt? 


UNIT a Rf wnUND ANn RELEASED 


0 


S-Ufi-, At+f- WTaNSFOSM -r • O.BY-* 1 WAtfriX 

« 

*** MAlRTX INVERSF DEMANDED 

enter $ur' h rec. « rows ? columns' matrix on unit n 

EXIT HRfce' INVERSF ON UNIT ^ 

SUB, aTR TRANSPOSED a «-8V* ? MATRIX 


SOLUTION VE r TOR NO'. 1 

— • ■ ' ... 

ITERATION SOLUTION- PROCEDURE INTERMEDIATE PRINTOUT 
INITUL GUEsS 

FILf: 7-Dfrftr.ftf-pNC OMrNfrD 'FAR REiDf^ tW UNIT a 

UNIT t\ REWOUND AND RELEASED 

FILE 70 6 7,BC-PNC HaS BEEN OPENED FOR READING ON UNIT « 

UNTT a REWOUND and released 

FILE 70«B,B(.-PNC Has been opened for reading ON UNIT U 

UNIT H REWOUND and RELEASED 

INFLUENCE OF UNIFORM and NON. UNIFORM FReESTRFaM 

and rotation rate vector on aircraft component i 
, i6#3?F*tO ,21S«6 Ea 10 .29693E+10 # 5*?<l9P*tO ,«2939f + 10 ,956«0F + 10 ,SU53i 

COMPONENT number j 

NtfNfrE# UF-ft f -H tNTg« -- - — - 

SOLUTIONS 

,10189**01 ,l*522F+0i ,J2019 e + 0I .2475iF+0i ,Ufli7F + 02 

INFLUENCE of AIRCRAFT COMPONENT j ON AIRCRAFT COMPONENT ? 


FIGURE 8.6-2 Output tor SCLN Code for Test Case No. 6 (Cont'd) 


FIGURES-67 


REPRODUCIBILITY OP THE 
ORIGINAL PAGE IS POOR 



.3I995E+10 .19929EM0 ,i20HE+l0 8 69707E+09 ,38721E + .09 ,20«30E,*09 


INFLUENCE OF UN | FORM aNO NON. UNIFORM FREESTREaM 

^ RA 5^6Q gCT ^ W ° N * lRCR ^ FT COMPONgN T ? 


COMPONENT NUMrER ? 


HWimun wr tPJtPMO* 3 

SOLUTIONS 

,25062E;0t ;S8«77E + 01 , ,!«'»9|E*02 


INFLUENCE OF aIRCRaFt COMPONEN t j ON AIRCRAFT COMPONENT 

g 1 0 1 3^E + 09 ,«67i2E + o8 ,19627**08 ,75323l-+07‘ 


INFLUENCE OF aIRCRaFt COMPONENT s ON aIRCRaFt COMPONENT 

,358a9Fi|0 g?2H0E+ 10 # 13208 f+10 .757UF+09 

INFLUENCE OF UNJFORM AND NON. UNIFORM FREESTREaM 
and ROTATION rate VECTOR on aircraft component s 


file 7073 ,S 0 «PN C 
CREATION T IME* 07 / 22/78 


0 33’32Wt * 

has bee w opened for reading on unit a 
U 10612* 


component NUMBER 3 


SOLUTIONS 

.90000E+OI ;10000E+02 


iteration 


\ LARGEST relative Errors 


,678b90E400 





ITERATION $ with If) ALLOWABLE 


INFLUENCE OF AIRCRAFT COMPONENT ? ON AIRCRAFT COMPONENT \ 


FIGURE 8.6-2 Output for SOLN Code for Test Case No. 6 (Cont'd) 


FIGURES— R8 



,542*>6EiO* ,i>«S4E*09 ,226aiE + 09 , <i239<JE + 09 f 7S521F + 09 .12927E + 10 .21347 

INFLUENCE OF a TRC.WaFT COMPONENT 3 ON AIRCRAFT COMPONENT \ 

,5238»F f 0'j ’.3i323F + 0P '.I3j7u>07 ,43967*401 '.12442 f;o 8 *,3t062F + o8 .70246 

COMPONENT NUMttEr 1 

NUMBER of el ements* s 

sot-oT-jo^a — - 

,1000C>F*01 .19998F + 0I ,30009ei0i .J9979F+01 t 50l0aF + 01 

INFLUENCE OF aIRC«aF7 COMPONENT , ON AIRCRAFT COMPONENT ? 

,l«S99F,i0 r ,!W90F + tO >6956F + 09 ,3ai80E + o9 ,20809F 4 09 ;|07*>0F4<)9 

|NFlHE*t£- OF' A'fRO^tFf COM^ONf^t ' 3 ON aIRCRaFT COMPONENT ? 

,14663P.*09 ;£8b«aF>09 'S2925F + 09 .93262F + 09 ,j5777E + iO r .2S756F + l0 

fOMPONFNT NMMpF-ft ? 

NUMBER OP f) FME NTS« 3 

SOLUTtONS 

,59RW,tTt' ;7«Ot7f40i .79989^ + 01 

INFLUENCE OF aIrCraFT COMPONENT t ON AIRCRAFT COMPONENT 3 

,239«8E4p8 ,959P7e + 07 ,54910*4.07 

INFLUENCE Of AIRCRAFT COMPONENT 2 ON AIRCRAFT COMPONENT * 

,23-M-9^jlft ,«52?5irvfr4- ,-«T>«?Rf*09 

COMPONENT number 3' 

NUMBER of El EHENT3* g ' 

SOI uttons 

,90009E<01 ;99f99|40l 

ITERATION ti LAfttCftT RELATIVE ErrOr* I S29720E.fl2 

ITERATION h L ATN1C6T RELATIVE ERROR* t l<lR 7 42F*02 

ITERATION i LARGEST RELATIVE ERROR* t 4l99l3E,03 

iteration r largest relative error* t ii7206p,03 


ITERATION f, tRRfflT RELATIVE ERROR* ,S26H7e.0« 

FILF 708«.SU.PNC/L!IIF Ha« BEEN OPENED FOR WRITING ON UNIT 3 


FIGURE 8.6-2 Output tor SCLN Code for Test Case No. b (Coat'd) 

REPRODUCIBILITY OF THE 
FIGURES - 69 ORKHMAL PAGE 18 - POOP 



CREATION TIME « 07/22/76 

tm*T — !~ W ItEO - A -N D RELEASED 


UNIT NO, 5 ^ 

• 1 UOOOOOf +0 1 ’.I999999E + 01 ;5000002E;Ot r ,3999987E*01 .S000063IT + Q 1 

FILE 7069.SO-PNC/LI8S H*S BEEN OPENED FOR WRITING ON UNIT 3 

-c reation t i »t • oy/ga/7 6 rntm? 

UNIT $ ENOFIlED AnO RELEASED 

-tunr-wv~T 

COMPONENT NO' 2 80lUTj0N8» 

,5999966E+01 ,7000009f* 01 # 7999996t;oi 

-ETtE 7090 ,SO*PNC/l IB? N AS BEEN OP E NE D' FOR WRI T IN G ON UN TT — 3 

CREATION TI«E • 07/22/76 UIJ«lS3 
UNIT 3 ENDFILED AND RELEASED 


UNIT NO, 3 

COMPONENT NO' 3 SOLUTIONS* 

^COOOOlE+Ol 79999 9 9 9 E #01 

NQ/ QF ITERATIONS TO CONVERGE* 9 LARGEST RELATIVE ERROR* ;326i«7E.0<i 


S O L UT I O N v g r To R NO' 1 



INITIAL GUESS 


FILE 7.0S6 .BC-PNC HaS BEEN OPENED FOR READING ON UNIT 9 
UNIT a REWOUND AND RELEASED 

UNIT <i REWOUND AND RELEASED ' ** 
FILE 7088,BC.PNC HaS BEEN OPENED FOR READING ON UNIT a 
UNI T 9 R E WO UN D AN D RELEASE D 


I n D^ R O TATiQN TJllT AN0 N0N,,UNIF0RM FREESTREaM 

* .33673EUO .97091E + 10 ,59387E;iO ,7taS8E*iO [s5877E;iO ,9l280E + iO , J0227K*! 


FIGURE 8,6-2 Output for SCLH Code for lest Case No. 6 (Cont’d) 


FIGURES-70 



COMPONENT NjiMpf R I 
WHM’# «F ELEMENTS* 5 
SOLUTIONS 

,20379f*oI .J7099E+01 .69038F+01 .99502E+01 .22833F+G2 


influence OF AIRCRAFT COMPONENT t 
,6399 JF,io ,39858EfI0 ,2«023E*t0 


ON AIRCRAFT COMPONENT ? 
.13991F+10 .77991F+09 ,90860F*09 


INFLUENCE OF 
AND ROTATtON 
.1081 IF ,11 


HNJFORM AND 

rate VECTOR 
.tt360E>t 1 


NON. UNIFORM 
ON AIRCRAFT 
, 1 21 ^ 3F+ 1 1 


FREESTREaM 

COMPONENT 

# !22<*9iUH 


? 

,129%0E + U 


,12T9ep + ll 


COMPONENT N!>Mt-Ff} p 
•NttHftFR-Of FtEftEnTSb • ? 

SOLUTIONS 

*50 1 pip « 0 ( .1J695F+02 *299fl2F*0;> 


INFLUENCE OF AIRCRAFT COMPONENT t DN AIRCRAFT COMPONENT J 

,2027lF.0 f ? ',93«23E + 0S ;3925JE + 08 .1U665F + O0 

INFLUENCE OF aIRCRaFt COMPONENT p ON aIRCRaFt COMPONENT s 

, 7 1 69 7f 4 i f) t ««279F + !0 .26915F4.10 *l?l«2F + tO 

INFLUENCE OF UNIFORM and non. uniform FREFSTRFaM 
and ROTATION rate VECTOR on AIRCRAFT COMPONENT i 

vmw+u 

FILE 707i.Sn.PNC uAS BEEN OPENED for READINo ON UNIT 9 

CRE ATION T I ,4 E*07/22/7to 11106*29 

UNIT n RfWfUJND AND RElEaSED 


COMPONENT NUMBER T 
NttMPftt -OF r, ff-WFNfSs - "2 • ' ■ 

SOLUTIONS 

,I 8000 F 4 02 ',20000F*02 

TTFRATION i LARGEST RELATIVE ERROR* 

-f tF"Rat ton — ^--t-Ai^tsr-REt'Artve-fwRoff* - 

iteration a largest relative errors 

ITERATION <, LARGEST RELATIVE ERROR* 


f 67fl6S0f>00 
, 195333F*0O 
,620199^.01 
, 189026E.01 


FIGURE 8.6-2 Output for SCLN Code for Test Case No. 6 (Cont'd) 


FIGURES-71 



ITERATION SOLUTION PROCEDURE intermediate printout 
ITERATION 5 WITH io ALLOWABLE 


# 1085Jg;o9 .22ROBE+OP t 95283E*09 .84&88E+09 


-component - - \ 
15l09E;t0 ,258S<Jg + io 


I N FLUE N CE O F aIEC R aFT C OM P ONENT 5 — O N AIRCRAFT CO MPO NENT T 

,I0«78f f 0fc >S6i|7E+08 '26J4 IE;o7 s 879S3E+07 r 4 24884fr;Q8 r # 62124E;<)8 


NUMBER OF elements* 

SOLUTIONS 

— r^ fr oo ft frTtn — ts w ^ e+o! — T .6 ooo Qg ;oi 


~ 4 ~TT T58g^01 VlOom + Og 

O n^ AI RCRA ET component ^ 2 

ON AIRCRAFT COMPONENT 2 

TtOt rS tEttO — T it S S 4 F» 1 0 — ^mTr?10 


TNFtUENCE OF aI&CRaFT COMPONENT 
»3p788g-^-tO — 722579E+ 10 r 


i 


INFLUENCE Op AIRCRAFT COMPONENT * 
lg 9 S* 5 Ei f><» — .B7288 E ^09 , 10 5 ftS m <T 


COMPONENT NUMBER 

1 1 1 1 jj n c ft n r 
NUHOCK wr 

SOLUTIONS 

,U988g;o2 1 ^OOSE + Oi ’.15998**02 


INFLUENCE OF AIRCRAFT COMPONENT i 
,109506*09 ,46895**08 ,19199**08 

INFLUENCE OF AIRCRAFT COMPONENT p 
,«68J9E + I 0 '.28369F + 10 ' 4 I664SE*10 


ON aircraft component 
,69820f*07 


ON AIRCRAFT COMPONENT 
.9J7S7E+09 


COMPONENT NuMrFr i 
NUMBER OF ELEMENTS* 2 

-SOtrUTl O N S — 

. 1800 IF +02 , 20000£*02 

ITE R ATIO N sr — LA R GEST R Et AT I ~VC~ £ RR OR « ^SgSTaOE^ O S 

'iteration h largest relative error, , 149742*. og 


ITERATION 7 LARGEST RELATIVE ERROR* .4I9913E.03 


FIGURE 8.6-2 Output for SOLN Code for Test Case No. 6 (Cant'd) 


,4269SE*10 
, 14049**04 
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w 

£71 

a 

SO 

KJ 

t/1 

t 

'-a 

ui 


9 


ITFRmON « LARGEST RELATIVE ERROR# 

iteration q largest relative error# 

FILE 70«8,$(i.pnc/L 16* Has been OPENED 
UNH - % EM>ML,Et> AND RELEASED 


, 1 1 7g06E*03 

,526i«7F*04 
FOR WRl T lN(i ON 


unit 


UNIT NO, S 

— rOXPONENT i SOttiljONSj 

pt t , rt if 00o0o0f * 01 ' 3RRR9R9F + 0I >00000^401 

r,i tos9,so-pnc/l !«.* has been opened for writing on unit 

UNIT ^ FAILED and released 


799997SE + 0: 


UNTT NO, ■* 

-C- O NP O NtNT- W^ - 2- SOtU-TTONS! 

» I 19R993I- + 0R , 1 U0OOO2F + 02 ,K99999EiO? 

Flu 7090,S(i-PNC/t.IB$ Has been OPENED FOR WRITING on unit 5 
UNIT % F>OHt fct) A no RELEASED 5 


.1000013*402 


UNIT NO, 3 

COHPONLNt-NO' i SOrUTtONSl 

• 1>'OOOAOF«02 • .2000000*402 

NO, OF ITERATIONS to CONVERGE# 9 LARGEST RELATIVE ERRORS 
♦ : Ail DONE f 

4 1 STOP 

— Mnf* 7fl -■ ■ ■ 


'.326l«7E-0a 


ilGUiJE 8.6 2 Output tor SCLN Code for Test Case No. 6 (Concluded) 



true 


SAMPLE CASE 7 FOR SQIN CODE, 


^ f iKf 1 ■ 


vTlUl'J 1 


THAT A ftECTANGULA# Matrix (17X1 0 ) ts 

!°K V ,I?-. IN !. TEAD ^°! a S0UARe matrix noxioi, 

FIRST, TWO AIRCRAFT COMPONENTS aR£ SOLVED 

by -direct householder procedure, then# the 


on THp FIRST Two IS COMPUTED BY THE 
ITERATION PROCEDURE, THE DIRECT SOLUTIONS 




VNRE 


AS INITIAL GUESSES TO THE THREE COMPONENT 
ITERATION SOLUTION, 


read influence $ub;hatrix files 

-P O R^~T-Hg -p y RgT TW O AI R CRAFT 

COMPONENTS'. 



since entire matrix fits in core, 

COMB COMMAND IS NOT NEEDED'. 



SOLVE bv HOUSEHOLDER PROCEDURE, 
AND define the SOLUTION 






j>-j jm *-* u 


remaining commands later, or he 
may CONTINUE. 




AIRCRAFT COMPONENT VH ITERATION 
PROCEDURE' 


FIGURE 8.7-1 Input for SOLN Cede for Test Case No. 7 
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3 1 

RERfcAU *LL INAlUENcE SUB.HaTRJ* 

-Fjt'fc-S* 

VNRE 

7077 

7078 

7079 
70*0 

■708? 

7083 
708U • 

708S 

THE COMP COMMAND May NOT BE USED 
WHH Fttf t TfRa T TON PROCEDURE. 

RP Rt AU A L l CONSTANT SUB. VECTOR riitS*. 

RCRE 

708#* 

7087 

7088 

' SOt-VE -^T frfrtwat 8 U>C 8 s MY 

HOUSEHOLDER PROCEDURE, DEFINE THE 
SOLUTION SUh.VeCTOR FILES, 1 Hf 
FRE^OFTFRMINEO APPROXIMATE 
SOLUTION SUr. VECTOR FILES, 

MAXIMUM mUMRER OF 1 1 F W A T ! OKiS , 

-■•A-M&UNT •• OP -fTV€R-R€iAXAUDN, 

AND AMOUNT OF INTERMEDIATE 
PRINT'UJT 4 

SOL V 

1 P 
7093 

W 

- 7095 

0 

t 7091 
£ 709 g 
10 

-+-.H7 - - 

0 

ALL D i i m f ' t 

STOP 


FIGUttE 8.7-1 Input for SCLN Code for Test Case No. 7 (Conducted) 
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POTFAN EQUATION SOLVING PROGRAM iSOLN), VERSION 1,1 

*8S«*a8S«*8t8«««**att«8**as888a888*888e8asa8B*a*8*8*tt 


dynamic MFMORy a lOQOO 


S/MPLE CASE 7 FOR SOtN cOOE'. 


THAT a RECTANGULAR MATRIX (17X10) is 
SOLVED INSTEAD OF A SQUARE MATRIX (10X10), 




FIRST, TWO AlRCRAfT COMPONENTS ARE SOLVED 
BY DIRECT HOUSEHOLDER PROCEDURE, then; THE 

ON THg first two is computed bt the 
ITERATION procedure; the DIRECT solutions 


♦I AS INITIAL GUESSES TO the THREE COMPONENT 

m iteration solution; 

-tl CNTt 

NCMP* P iTERATa 0 
♦ I RFiD INFLUENCE SUBwMaTRIX- FILES 


♦ I COMPONENTS’, 

♦ I .VNRE 

FILE 70 7- T t V^ PNC H*S“Bf£N OPENED FOR READING ON UN I T \ 

CREATION TIMFtOT/07/76 OOtOOJOl 
UNIT 1 REWOUND AND RELEASED 


CREATION TjMEtOT/07/76 04!O*8Ol 

UNIT 1 REWOUND AND RELEASED 

CREATION TlHEa07/07/76 0410*101 

UNIT 1 REWOUND AND RELEASED 


CREATION TIMEa07/07/76 0410*101 

UNIT i REWOUND AND RELEASED 


t | w A 'Yt u L n 1 1 n|» i rv* n r * t w 4. 

♦I C0M8 COMMAND IS NOT NcEDcOl 


DN-tWtT- I 


FIGURE 8.7-2 Output for SOLN' Code for Test Case No. 7 


FIGURES-76 



+ f RE AD CONSTANT SUP-VECTOR FILES'. 

• *t -'HCRfe “ ' - • - 

Fll.f: 708b t RC«PNC Has PEER OPgNEri FOR R£AdI ng ON UNIT 2 

CREATION ?1 M F*07/07/76 OUlQPlol 

UNIT 7 REWOUND AND RELEASED 

FILE 70«7,RC*PNC H*S BEEN OPENED FOR READING ON UNIT ? 

CRFATiON 1 ihEi, 07/07/76 0«t09l0i 

UNIT - * Ri-wnttND Atnr RELEASED 
+| SOLwE RY HOUSEHOLDER RROrEOURE, 

+ 1 A fJ D DEFINE THE SOLUTION 

1 1 SUr^vfCTOR OUTPUT FILES'. 

M SOlV 
METHODS 
IRS* 0 

IPRJNTb '• 

( OU T A # OUT AnT.Js.J.NCHP)* 

7091 r. 

709P 

ENTER SUB. PART *FRAGS» If) tnv* «1 

'HATfrtr tO, '0^0 F R OM "UNIT |3 USING |ft READS AND' 3rt SKTPS 
NO .PARTITIONING RFOUlRfD 
FXTT PART 

ENTER SUr' HOtlSF. 13 ROWS 8 r OtUMNS ? SOLUTIONS 

EXIT HOUS*’ 

OUTPUT 

SOLUTION tjKF fSEC'. ) .025 

I ■ 

*** SOLUTION TIME Is THAT FOR ALL SOLUTIONS 


SOLUTION VErTQR NO' 1 


FILE 70^1 ,S'i.pNC/l!B$ h aS BEEN OPENED FOR WRITING ON UNIT 5 

CREATION TT'-F « 07/22/79 UU 6 HI 


UNIT NO, 5 

COMPONENT NU. 1 SOLUTIONS! 

t 99H}*62E*00 ,200T92SE*01 *298tj8t|5E 4 0 1 ’«0?I&2JF*O1 4 «75 1 636F+01 

UNIT 3 E w OF It ED AND RELEASED 


FIGURE 8.7-2 Output for SCLN Code for Test Case No. 7 {Cont*d) 


FIGURES-77 


REPROD UCIBlLITy OR UiB 

rHRJ75S>IA l'_ T1 A 



FILE 7OV2 t SO«pNe/U0$ HaS BEEN OPENED FOR WRITING ON UNIT J 
creat io n time m o7/g2/7 t — i n | ft i aa 


UNIT NO, 5 


,feSn 6 iOE + 01 ;«802658Fi0l 

UNIT 3 ENOFILEO AND RELEASED 


l5l8873’E'f02 


SOLUTION V£r TOR NO; 2 


FILE 7091 .SO-PNC/LIBS HaS BEEN OPENED FOR WRITING ON UNIT 3 


UNIT NO, 3 

COMPONENT NO' 1 SOLUTIONS! 

; ~ f 99 t>7 ' 3g £ + 01 ,-4fl TS84Rrm — H*r*S**5t + 01 ‘ SnirttWf *t % 

unit 3 endfi,eo and released 

FILE 7092,Sn«PNC/LlB* HAS BEEN OPENED FOR WRITING ON UNIT 3 


UNIT NO, 3 

COMPONENT No', 2 . 80tUTl0N$j 

r HT* T g gE >Og . 9 fe0930«E + 01 ,30377a6g;02 

UNIT 3 ENOFILEO AND RELEASED 
tl USER COULD STOP HERE AND PERFORM 


— rr- 

♦ i 

♦ « 

Remaining commands laTeR, or hf 

mav continue'. 

DETERMINE INFLUENCE OF THIRo 



♦ i 

♦ i 

♦ i 

AIRCRAFT COMPQW|TIt — VTA'~ITeRaTiON 
PROrEDURE, 

CNTt 



flcnn 1 i I r K A 1 4 J 

♦i REREAD ALL INFLUENCE SUB.MaTRJX 

♦ 1 FILES’, 

— n vnre - - - — - 


FILE 7077.VN.PNC HaS been opened for reading ON UNIT t 

CREATION TIME*07/07/7S 0<ll09»0l 

-tmtr — r -*w t mp-gN t r re l e ased : 

FILE 7078, VN*RNc HaS BEEN OPENED FOR READING ON UNJT j 

CREATION TiME«07/07/76 OSIORgfll 

tmfT 1 REWOUND AN D R ELEASED 

FILE 7079.VNUPNC HaS BEEN OPENED FOR READING ON UNIT j 


FIGURE 8.7-2 Output for SCLN Code for Test Case No. 7 (Cont*d) 
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CREATION t IMf:s07V07/76 Oa*0R|Ol 

UNH- ANfl RfifcASt^ " 

FILE 7060, VN-PMC Has BEEN opened fop reading ON UNIT J 

CREATION 1 ! M F*07/07/76 oaiO^IOl 

UNIT 1 REWOUND AND RELEASED' 

FILE 7o»1,VN«PNc H a S BEEN O p f N£P FOR READING ON UNIT j 

CREATION T jmF *07/0 7/76 0«109|<U 

-Ufrff-- -t— RtWMTNfT-AND' RELEASED - 

FILf 7 q62,VN*PNc Has REEN OPENED FOP READING ON UNIT j 

CREATION I I *07/07/76 OdlORJOl • ■ 

UNIT \ Rewound -AND RELEASED 

FILf 706S,VN^PN(; H a S BEEN OPENED FOR REaDING ON UNIT j 

CREATION TT m F*G7/07/76 0 « * 0^ » 0 I 

ONi V t" ■RtwfHHrf)' - A NO* R E t E A SFD 

FILE 7()«<I,VN.PNC Has BEEN OPENED FOR reading ON UNIT t - 

CREATION T thE«07/07/76 0«J0R»M 

UNIT 1 REWOUND AND RELEASED 

FILf 7 o05.VN«PNC h a S BEEN OPENED FOR REaDING ON UNIT ( 

CREATION time * 07/07/76 0«*OR»Ol 

•UN-IT t TrtnmtJND"ANo RftfTSED 
f* TMf comb command MAY NOT be used 

+ t WITH THE ITERATION PROCEDURE*. 

+ t RFRFAD all constant SUB. VECTOR FILES'. 

m rcrf 

file 7066, RC. PNC HAS been opened for reading ON UNIT p 

CREAMDN ^i+tt»D-7/T>7/76 0«toRtol’ 

UNIT ? REWOUND aND RELEASED 

FILE 70^7 * PC .PNC h A S BEEN OPENED FOR REAplNG ON UNIT 2 

CREATjON Tt^f *07/07/76 O«»0RlOt 

UNIT ? RFwnuwD AND RElEaSEO 

FILE ?o6«*BC.PNC ' H A S BEEN OPENED FOR READING ON UNIT p 

CREATION T THE *07 /0 7 /T6 1 OatO**T>1 
UNIT p REWOUND aND RELEASED 
M SOLVE diagonal BLOCKS by. 

M HOUSEHOLDER PROCEDURE, DEFINE THf 

M SOLUTION SUB- VECTOR FILE'S, THf 

PRE*DETERMfN£D APPROXIMATE 
•-♦ r - - S^U+fTTN'SUn-rVfC'TBR FILES',’ ’ 

♦j maximum number of iterations' 

+1 amount OF OVeR-RpLAXATlON, 

M A#»o AMOUNT- OF TNTERMf DT A Tf 

M PRINTOUT’. 


FIGURE 8.7-2 Output for SGLN Code for Test Case No. 7 (Coat'd) 
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REPRODUCIBILITY OF THE 
QglSMAL PAGE IS POOS 



+ 1 S01.V 

M trT 'H Opa \ 

IPS# 2 
IPRINT# n 

— TOUT A » o U t-a-n Tvi-tr^yNC-MPf* — 

7093 n 

709<l n 

TOTS n 

(INPSrl^PsNTf !*1 iN0HP>, 

7091 o 

— mi j> 

0 (i 

MAX jt 1 a l n 

PM ECAal 1 10 0 

ITERP. o 



ENTER SUB' HRfC, 7 ROWS 

EXIT HREC^ flTV- g ft S g ON-UNIT n 

SUB, ATR TRANSPOSED A 7 » 8 V« 

*** MATRIX inverse demanded 

•EN TEr SUB 1 , n r ec ’. t r r OH S 

exit, hrec' inverse on unit 12 

SUB. aT R T-Ra NSPOS ED-A 

*** MATRIX inverse demanded 

ENTER SU B ' MRrc‘. « Rows 

EXIT HRfcc' INVERSE ON UNIT t2 

SUB, atr transposed A <UBY« 


s columns; matrix on unit ti’ 


<5 MATRfX 


s-columns'. ma t ri x on uni r ra — 

— 5 — HJ TRI X 


z columns; matrix on unjt 11 


? matrix 


solution v e?tor no; 1 

m&mt* * - rn - &mwmrn ' ii ' Wm m mmmmm 

FILE 7086 .BC-PNC HaS SEEN OPENED FQR READING ON UNIT a 


FIGURE 8.7-2 Output for SOLN Cods for Test Case No. 7 (Cont*d) 
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UNIT « REWOUND AND H£l.£ aSED 

FILft 7t>6?,KC*RNe H*9 BEEN OF£NeO"+'flf? READING ON UNIT a 

unit u rl wound and wEt E ased 

FRF 708B.BC.PNC HaS BEEN OPENED FOR READING ON UNIT a 

unit a REWOUND AND released 

FILE 70R1.SD-PNC H A S BEEN OPENED FOR READING ON UNIT a 

creation t i' j e* 07/22/76 Uii 6 iaj 

UNIT « REWtTttND" «Np RELEASED 

FILE 70R3.sn.pwC H 6 S BEEN OPENED FOR READING ON UNjT « 

creation i m-#o7/?j?/?6 imin# 

UNIT 'i REWOUND and RELEASED 

iteration , largest relative error* .oToiise+oo 

ITERATION ? LARGEST RELATIVE ERROR* ,1«7M!F*00 

ttfratiun x largest relative errors ,?S 886 of-<u 

ITERATION ti URGES! RELATIVE ERRORS t 7S0a95E.03 

ITERATION s LARGEST RELATIVE ERROR* .153576E-02 

iteration * largest relative error. .swsie-oi 

iteration 7 LARGEST RELATIVE ERROR* .aOSaOUF.OW 

flit 70^V-S rr **Nt/U&S HAS BEEN OP£NEO'~Ft>R WRITING ON UNIT ? 

CREATION Time a 07/22/76 ntlft|I 6 

UNIT 3 ENOrTLFD AND RELEASED 

UNT7 NO, 3 

COMPONENT NO" 1 SOLUTIONS, 

9979RR8£TD0 ’?00000IF+0I *?RRR9R«r + oi ’ 90 O0DO 7F ♦ 0 1 , ?£ * 0 1 

FILE 709<i, sOoPNC/LlBt * W AS BEEN OPENED FOR WRITING ON UnIT 3 
CRFATION T I mE * 07/22/76 HUBflB 
UNI! 3 ENOULFD AND RELEASED 

UNIT NO, 3 

-- COMPONENT ' NO' 2 SOjUTTDNSl 

,6000016E*01 .6999996f, + 0i .8OO0OO1F+O1 

FILE 709«5,SO*PNC/l:iB$ has BEEN OPENED FOR WRITING ON UNIT 5 

CREATION TJmE s 07/22/76 j tltlftli* 

UNIT 3 ENDULED AND RELEASED 


FIGURE 8.7-2 Output for SCLN Code for lesfc Case No. 7 (Cont’d) 
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SSHaSKS 1 



UNIT NO, i 

CO MP ONENT no; - 1 ' -fr OtUTrO N Sj 

,9000000E + Ot MoOOOOOF + 02 


-N6- r - 


~«052o«e*cr« 


SOLUTION vEfTOR NO' 2 


UNIT a REWOUND AND RELEASED 
FILE 70*7, BC- PNC Ha8 BEEN OPENED FOR READING ON UNIT 

U N IT ti -R E W OU N D A N p - RELEASED 


— »- 


FILE 70S8,BC-PnC 
UNIT a RF WOUND AND 

-F ILE 7p9l.SO.P Nt- 

CREATjoN T r ME*07/22/76 U»16|4S 

UNIT 4 REWOUND AND RELEASED 

■ FILE 7 p 9 ?,SO*PN C H fl g BfgN- 

CREATION TTME*07/22/76 Ull6*44 

UNIT a REWOUND AND RELEASED 


has been opened fqr reading on unit 4 

RELEASED 

u 1 f i n f. y . k t ft fir H r n f _ n r . n T UF rt A t t (Ut t 

HAir tt.N rOR KCAui NR ON uN IT <j 


TTPg NgD FOR READI N G ON UN f~T 4 


ITERATION 1 LARGEST RELATIVE ERROR* a a701t9E+00 

ITERATION r - LA R G E ST R gLATlVg E R RO R* ,14 7 611^ 0 

ITERATION * LARGEST RELATIVE ERROR* B 2S8860E.0t 

iteration u largest relative error 3 b 7so«9se«o2 

ITE R ATION 5 — LARGEST RELATIVE ER R O R* , 1 53S78E..02 

iteration h largest relative error b ,25985ie*03 

ITERATION t LARGEST RELATIVE ERROR* .405204F.04 

FILE 709S,SO-PNC/UB$ Has BEEN OPENEO for WRITING ON UNIT i 

■Utm 3 EN p FIlf P AND -Rf LEA -W * — — 

UNIT no, 5 

COMPONENT ftQ\ 1 — SO LU T IONS 1 — 

*2QOOOOOE*Oi .4000003E + 01 ,5999996£;0t .8000015F+01 1 9999P25t ♦ 0 1 


FIGUBE 8.7-2 Output for SOLN Code for Test Case No. 7 (Cont*d) 
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FILE 709q, sfWpNC/LIBS H*S BEEN OPENED FOR WRITING ON UNIT i 
ttf o ~ +KD- RCt'eA 8E0 

UNIT NO, 3 

• OOwPONEnI VO; ? 30| U T 1 0N8 1 

,1?0000 JF+O? , 1399999f* 0? , 1 bOOOO 0E+ 0? 

FILE 709 ^, SO.PNc/LIBt m*S BEEN OPENED FOR WRITING ON UNIT } 

-"■WfHhf — T“€^F'ft-F'0-*'A-NO"RFlfa«eD 

UNIT NO, 3 

tOwFOWfcN! f?t>' 1 SOiUTlONSj 

,1''0000.0F40? .2000000F+02 

~-NO'/^r-Ht^rfxO^^--TO-CTmveReF» 7 tARGFSt RELATIVE ERRORS «0520<1E.09 
♦ i aii done, 

*t STOP 
-STOF 777 


FIGURE 8.7-2 Output for SCLN Code for lest Case No. 7 (Concluded) 



